Electronic Transport in Metallic Systems and 
Generalized Kinetic Equations* 



This paper reviews some selected approaches to the description of transport properties, 
mainly electroconductivity, in crystalline and disordered metallic systems. A detailed quali- 
tative theoretical formulation of the electron transport processes in metallic systems within a 
model approach is given. Generalized kinetic equations which were derived by the method of 
the nonequilibrium statistical operator are used. Tight-binding picture and modified tight- 
binding approximation (MTBA) were used for describing the electron subsystem and the 
electron-lattice interaction correspondingly. The low- and high-temperature behavior of the 
resistivity was discussed in detail. The main objects of discussion are nonmagnetic (or param- 
agnetic) transition metals and their disordered alloys. The choice of topics and the emphasis 
on concepts and model approach makes it a good method for a better understanding of the 
electrical conductivity of the transition metals and their disordered binary substitutional al- 
loys, but the formalism developed can be applied (with suitable modification), in principle, to 
other systems. The approach we used and the results obtained complements the existent theo- 
ries of the electrical conductivity in metallic systems. The present study extends the standard 
theoretical format and calculation procedures in the theories of electron transport in solids. 



Keywords: Transport phenomena in solids; electrical conductivity in metals and alloys; 
transition metals and their disordered alloys; tight-binding and modified tight-binding approx- 
imation; method of the nonequilibrium statistical operator; generalized kinetic equations. 



* International Journal of Modern Physics B (IJMPB), Volume: 25, Issue: 23-24 (2011) p. 3071-3183 



A. L. Kuzemsky 
Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 
141980 Dubna, Moscow Region, Russia. 
E-mail:kuzemsky@theor. jinr.ru 




Abstract 



Contents 



1 Introduction 

2 Metals and Nonmetals. Band Structure 

3 Many- Particle Interacting Systems and Current operator 

4 Tight-Binding and Modified Tight-Binding Approximation |o 

4.1 Tight-binding approximation 11 

4.2 Interacting electrons on a lattice and the Hubbard model 

4.3 Current operator for the tight-binding electrons 

4.4 Electron-lattice interaction in metals 

4.5 Modified tight-binding approximation 

5 Charge and Heat Transport 

5.1 Electrical resistivity and Ohm law 

5.2 Drude-Lorentz model 

5.3 The low- and high-temperature dependence of conductivity 

5.4 Conductivity of alloys 

5.5 Magnetoresistance and the Hall effect 

5.6 Thermal conduction in solids 

6 Linear Macroscopic Transport Equations 

7 Statistical Mechanics and Transport Coefficients 

7.1 Variational principles for transport coefficients 

7.2 Transport theory and electrical conductivity 

8 The Method of Time Correlation Functions 

8.1 Linear response theory 

8.2 Green functions in the theory of irreversible processes 

8.3 The electrical conductivity tensor 

8.4 Linear response theory: pro et contra 

9 The Nonequilibrium Statistical Operator Method and Kinetic Equations 

10 Generalized Kinetic Equations and Electroconductivity 

10.1 Basic formulas 

10.2 Electrical conductivity 

11 Resistivity of Transition Metal with Non-spherical Fermi Surface 

11.1 Generalized kinetic equations 

11.2 Temperature dependence of R 

11.3 Equivalence of NSO approach and Kubo formalism 

11.4 High-temperature resistivity and MTBA 

12 Resistivity of Disordered Alloys 

13 Discussion 



1 



1 Introduction 



Transport properties of matter constitute the transport of charge, mass, spin, energy and momen- 
tum. 1 " 8. It has not been our aim to discuss all the aspects of the charge and thermal transport 
in metals. We are concerned in the present work mainly with some selected approaches to the 
problem of electric charge transport (mainly electroconductivity) in crystalline and disordered 
metallic systems. Only the fundamentals of the subject are treated. In the present work we aim 
to obtain a better understanding of the electrical conductivity of the transition metals and their 
disordered binary substitutional alloys both by themselves and in relationship to each other within 
the statistical mechanical approach. Thus our consideration will concentrate on the derivation of 
generalized kinetic equations suited for the relevant models of metallic systems. 
The problem of the electronic transport in solids is an interesting and actual part of the physics of 
condensed matterP^H includes the transport of charge and heat in crystalline and disordered 
metallic conductors of various nature. Transport of charge is connected with an electric current. 
Transport of heat has many aspects, main of which is the heat conduction. Other important as- 
pects are the thermoelectric effects. The effect, termed Seebeck effect, consists of the occurrence 
of a potential difference in a circuit composed of two distinct metals at different temperatures. 
Since the earlier seminal attempts to construct the quantum theory of the electrical, thermaP^^ 
and thermoelectric and thermomagnetic transport phenomena,^ there is a great interest in the 
calculation of transport coefficients in solids in order to explain the experimental results as well 
as to get information on the microscopic structure of materials F2H35] 

A number of physical effects enter the theory of quantum transport processes in solids at various 
density of carriers and temperature regions. A variety of theoretical models has been proposed 
to describe these effectsP^^IIISEaEilEsllMHIl Theory of the electrical and heat conductivities 
of crystalline and disordered metals and semiconductors have been developed by many authors 
during last decades QHsKqBMHH There exist a lot of theoretical methods for the calculation of 
transport coefficients QSiQMSHSDlSHlID as a ru j e having a fairly restricted range of validity and ap- 
plicability. In the present work the description of the electronic and some aspects of heat transport 
in metallic systems are briefly reviewed, and the theoretical approaches to the calculation of the 
resistance at low and high temperature are surveyed. As a basic tool we use the method of the 
nonequilibrium statistical operator3ZE2(NSO). it provides a useful and compact description of the 
transport processes. Calculation of transport coefficients within NSO approach^ was presented 
and discussed in the author's workp^ The present paper can be considered as the second part of 
the review articled The close related works on the study of electronic transport in metals are 
briefly summarized in the present work. It should be emphasized that the choice of generalized 
kinetic equations among all other methods of the theory of transport in metals is related with 
its efficiency and compact form. They are an alternative (or complementary) tool for studying of 
transport processes, which complement other existing methods. 

Due to the lack of space many interesting and actual topics must be omitted. An important and 
extensive problem of thermoelectricity was mentioned very briefly; thus it has not been possible 
to do justice to all the available theoretical and experimental results of great interest. The ther- 
moelectric and transport properties of the layered high-T c cuprates were reviewed by us already 
in the extended review article.^ 

Another interesting aspect of transport in solids which we did not touched is the spin transport!^' 
The spin degree of freedom of charged carriers in metals and semiconductors has attracted in last 
decades big attention and continues to play a key role in the development of many applications, 
establishing a field that is now known as spintronics. Spin transport and manipulation in not 
only ferromagnets but also nonmagnetic materials are currently being studied actively in a vari- 
ety of artificial structures and designed new materials. This enables the fabrication of spintronic 
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Tablel. Five Categories of Crystals 



Table 1: Five Categories of Crystals 



Type of Crystal 


substances 


ionic 


alkali halides, alkaline oxides, etc. 


homopolar bounded (covalent) 


diamond, silicon, etc. 


metallic 


various metals and alloys 


molecular 


Ar, He, 2 , H 2 , CH 4 , etc. 


hydrogen bonded 


ice, KH2PO4, fluorides, etc. 



properties on intention. A study on spintronic device structures was reported as early as in late 
sixties. Studies of spin-polarized internal field emission using the magnetic semiconductor EuS 
sandwiched between two metal electrodes opened a new epoch in electronics. Since then, many 
discoveries have been made using spintronic structures .^"^ Among them is giant magnetoresistance 
in magnetic multilayers. Giant magnetoresistance has enabled the realization of sensitive sensors 
for hard-disk drives, which has facilitated successful use of spintronic devices in everyday life. 
There is big literature on this subject and any reasonable discussion of the spin transport deserves 
a separate extended review. We should mention here that some aspects of the spin transport in 
solids were discussed by us in RefspSES 

In the present study a qualitative theory for conductivity in metallic systems is developed and 
applied to systems like transition metals and their disordered alloys. The nature of transition 
metals is discussed in details and the tight-binding approximation and method of model Hamil- 
tonians are described. For the interaction of the electron with the lattice vibrations we use the 
modified tight-binding approximation (MTBA). Thus this approach can not be considered as 
the first-principle method and has the same shortcomings and limitations as describing a transi- 
tion metal within the Hubbard model. In the following pages, we shall present a formulation of 
the theory of the electrical transport in the approach of the nonequilibrium statistical operator. 
Because several other sections in this review require a certain background in the use of statistical- 
mechanical methods, physics of metals, etc., it was felt that some space should be devoted to this 
background. Sections 2 to 8 serves as an extended introduction to the core sections 9-12 of the 
present paper. Thus those sections are intended as a brief summary and short survey of the most 
important notions and concepts of charge transport (mainly electroconductivity) for the sake of a 
self-contained formulation. We wish to describe those concepts which have proven to be of value, 
and those notions which will be of use in clarifying subtle points. 

First, in order to fix the domain of study, we must briefly consider the various formulation of the 
subject and introduce the basic notions of the physics of metals and alloys. 

2 Metals and Nonmetals. Band Structure 

The problem of the fundamental nature of the metallic state is of long standing It is well 
known that materials are conveniently divided into two broad classes: insulators (nonconducting) 
and metals (conducting) EMSE] More specific classification divided materials into three classes: 
metals, insulators, and semiconductors. The most characteristic property of a metal is its ability 
to conduct electricity. If we classify crystals in terms of the type of bonding between atoms, they 
may be divided into the following five categories (see Table 1). 

Ultimately we are interested in studying all of the properties of metals At the outset it is 
natural to approach this problem through studies of the electrical conductivity and closely related 
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Figure 1: Schematic form of the band structure of various metals 



problem of the energy band structure* 
The energy bands in solids^®^ represent the fundamental electronic structure of a crystal just 
as the atomic term values represent the fundamental electronic structure of the free atom. The 
behavior of an electron in one-dimensional periodic lattice is described by Schrodinger equation 

where V is periodic with the period of the lattice a. The variation of energy E(k) as a function 
of quasi-momentum within the Brillouin zones, and the variation of the density of states D(E)dE 
with energy, are of considerable importance for the understanding of real metals. The assumption 
that the potential V is small compared with the total kinetic energy of the electrons (approxima- 
tion of nearly free electrons) is not necessarily true for all metals. The theory may also be applied 
to cases where the atoms are well separated, so that the interaction between them is small. This 
treatment is usually known as the approximation of " tight binding" P21 In this approximation the 
behavior of an electron in the region of any one atom being only slightly influenced by the field 
of the other atomsP^^] Considering a simple cubic structure, it is found that the energy of an 
electron may be written as 

E(k) = E a — t a — 2t/s(cos(k x a) + cos^a) + cos(^a)), (2-2) 

where t a is an integral depending on the difference between the potentials in which the electron 
moves in the lattice and in the free atom, and tp has a similar significance^"^ (details will be 
given below). Thus in the tight-binding limit, when electrons remain to be tightly bound to their 
original atoms, the valence electron moves mainly about individual ion core, with rare hopping 
from ion to ion. This is the case for the cl-electrons of transition metals. In the typical transition 
metal the radius of the outermost d-shell is less than half the separation between the atoms. As 
a result, in the transition metals the d-bands are relatively narrow. In the nearly free-electron 
limit the bands are derived from the s- and p-shells which radii are significantly larger than half 
the separation between the atoms. Thus, according this simplified picture simple metals have 
nearly- free-electron energy bands (see Fig.l). Fortunately in the case of simple metals the com- 
bined results of the energy band calculation and experiment have indicated that the effects of the 
interaction between the electrons and ions which make up the metallic lattice is extremely weak. 
It is not the case for transition metals and their disordered alloys 

An obvious characterization of a metal is that it is a good electrical and thermal conductor JIHIIMIMIH] 
Without considering details it is possible to see how the simple Bloch picture outlined above ac- 
counts for the existence of metallic properties, insulators, and semiconductors. When an electric 
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current is carried, electrons are accelerated, that is promoted to higher energy levels. In order 
that this may occur, there must be vacant energy levels, above that occupied by the most en- 
ergetic electron in the absence of an electric field, into which the electron may be excited. At 
some conditions there exist many vacant levels within the first zone into which electrons may be 
excited. Conduction is therefore possible. This case corresponds with the noble metals. It may 
happen that the lowest energy in the second zone is lower than the highest energy in the first 
zone. It is then possible for electrons to begin to occupy energy contained within the second zone, 
as well as to continue to fill up the vacant levels in the first zone and a certain number of levels 
in the second zone will be occupied. In this case the metallic conduction is possible as well. The 
polyvalent metals are materials of this class. 

If, however, all the available energy levels within the first Brillouin zone are full and the lowest 
possible electronic energy at the bottom of the second zone is higher than the highest energy in 
the first zone by an amount AE, there exist no vacant levels into which electrons may be excited. 
Under these conditions no current can be carried by the material and an insulating crystal results. 
For another class of crystals, the zone structure is analogous to that of insulators but with a very 
small value of AE. In such cases, at low temperatures the material behaves as an insulator with 
a higher specific resistance. When the temperature increases a small number of electrons will be 
thermally excited across the small gap and enter the second zone, where they may produce metal- 
lic conduction. These substances are termed semiconductor s JHHSMIS] anc j their resistance decreases 
with rise in temperature in marked contrast to the behavior of real metals (for a detailed review 
of semiconductors see 

RefsPEH]). 

The differentiation between metal and insulator can be made by measurement of the low frequency 
electrical conductivity near T = K. For the substance which we can refer as an ideal insulator 
the electrical conductivity should be zero, and for metal it remains finite or even becomes infi- 
nite. Typical values for the conductivity of metals and insulators differ by a factor of the order 
10 10 — 10 15 . So big difference in the electrical conductivity is related directly to a basic difference 
in the structural and quantum chemical organization of the electron and ion subsystems of solids. 
In an insulator the position of all the electrons are highly connected with each other and with the 
crystal lattice and a weak direct current field cannot move them. In a metal this connection is not 
so effective and the electrons can be easily displaced by the applied electric field. Semiconductors 
occupy an intermediate position due to the presence of the gap in the electronic spectra. 
An attempt to give a comprehensive empirical classification of solids types was carried out by 
Zeitz^ and KittelP^ Zeitz reanalyzed the generally accepted classification of materials into three 
broad classes: insulators, metals and semiconductors and divided materials into five categories: 
metals, ionic crystals, valence or covalent crystals, molecular crystals, and semiconductors. Kittel 
added one more category: hydrogen-bonded crystals. Zeitz also divided metals further into two 
major classes, namely, monoatomic metals and alloys. 

Alloys constitute an important class of the metallic systems E3im[55][ 56l | 59H 6i] This class of sub- 
stances is very numerous] EE9H61] me £ a i a lloy is a mixed material that has metal properties and 
is made by melting at least one pure metal along with another pure chemical or metal. Examples 
of metal alloys Cu — Zn, Au — Cu and an alloy of carbon and iron, or copper, antimony and 
lead. Brass is an alloy of copper and zinc, and bronze is an alloy of copper and tin. Alloys of 
titanium, vanadium, chromium and other metals are used in many applications. The titanium 
alloys (interstitial solid solutions) form a big variety of equilibrium phases. Alloy metals are usu- 
ally formed to combine properties of metals and the exact proportion of metals in an alloy will 
change the characteristic properties of the alloy. We confine ourselves to those alloys which may 
be regarded essentially as very close to pure metal with the properties intermediate to those of 
the constituents. 

There are different types of monoatomic metals within the Bloch model for the electronic struc- 
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ture of a crystal: simple metals, alkali metals, noble metals, transition metals, rare-earth metals, 
divalent metals, trivalent metals, tetravalent metals, pentavalent semimetals, lantanides, actinides 
and their alloys. The classes of metals according to crude Bloch model provide us with a simple 
qualitative picture of variety of metals. This simplified classification takes into account the state 
of valence atomic electrons when we decrease the interatomic separation towards its bulk metallic 
value. Transition metals have narrow d-bands in addition to the nearly-free-electron energy bands 
of the simple metalsP^^ T n addition, the correlation of electrons plays an essential roleP^^^H 
The Fermi energy lies within the (i-band so that the d-band is only partially occupied. Moreover 
the Fermi surface have much more complicated form and topology. The concrete calculations 
of the band structure of many transition metals (Nb, V, W, Ta, Mo, etc.) can be found in 
Refg p[32H35]|53l[63H66] and in Landolt-Bornstein reference booksP 1221 

The noble metal atoms have one s-electron outside of a just completed d-shell. The <i-bands of 
the noble metals lie below the Fermi energy but not too deeply. Thus they influence many of 
the physical properties of these metals. It is, in principle, possible to test the predictions of the 
single-electron band structure picture by comparison with experiment. In semiconductors it has 
been performed with the measurements of the optical absorption, which gives the values of various 
energy differences within the semiconductor bands. In metals the most direct approach is related 
to the experiments which studied the shape and size of the Fermi surfaces. In spite of their value, 
these data represent only a rather limited scope in comparison to the many properties of metals 
which are not so directly related to the energy band structure. Moreover, in such a picture there 
are many weak points: there is no sharp boundary between insulator and semiconductor, the the- 
oretical values of AE have discrepancies with experiment, the metal-insulator transition^ cannot 
be described correctly, and the notion " simple" metal have no single meaningP^ The crude Bloch 
model even met more serious difficulties when it was applied to insulators. The improved theory 
of insulating state was developed by KohrPSI within a many-body approach. He proposed a new 
and more comprehensive characterization of the insulating state of matter. This line of reasoning 
was continued further in RefsP^^"^ on a more precise and firm theoretical and experimental 
basis. 

Anderson^ gave a critical analysis of the Zeitz and Kittel classification schemes. He concluded 
that "in every real sense the distinction between semiconductors and metals or valence crystals as 
to type of binding, and between semiconductor and any other type of insulator as to conductivity, 
is entirely artificial; semiconductors do not represent in any real sense a distinct class of crystal"^ 
(see, however Refs.^'^'^211^'). Anderson has pointed also the extent to which the standard 
classification falls. His conclusions were confirmed by further development of solid state physics. 
During the last decades a lot of new substances and materials were synthesized and tested. Their 
conduction properties and temperature behavior of the resistivity are differed substantially and 
constitute a difficult task for consistent classification^ ( se e Fig.l). BokipS carried out an interest- 
ing analysis of notions "metals" and "nonmetals" for chemical elements. According to him, there 
are typical metals (Cu, Au, Fe) and typical nonmetals (O, S, halogens), but the boundary be- 
tween them and properties determined by them are still an open question. The notion "metal" is 
defined by a number of specific properties of the corresponding elemental substances, e.g. by high 
electrical conductivity and thermal capacity, the ability to reflect light waves(luster), plasticity, 
and ductility. Bokij emphasizes,^ that when defining the notion of a metal, one has also to take 
into account the crystal structure. As a rule, the structure of metals under normal conditions are 
characterized by rather high symmetries and high coordination numbers (c.n.) of atoms equal to 
or higher than eight, whereas the structures of crystalline nonmetals under normal conditions are 
characterized by lower symmetries and coordination numbers of atoms (2-4). 
It is worth noting that such topics like studies of the strongly correlated electronic systems,^ 
high-T c superconductivity,^ colossal magnetoresistance^ and multiferroicity^ have led to a new 
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Figure 2: Resistivity of various conducting materials (from RefP2) 

development of solid state physics during the last decades. Many transition-metal oxides show 
very large ("colossal") magnitudes of the dielectric constant and thus have immense potential 
for applications in modern microelectronics and for the development of new capacitance-based 
energy-storage devices. These and other interesting phenomena to a large extend first have been 
revealed and intensely investigated in transition-metal oxides. The complexity of the ground states 
of these materials arises from strong electronic correlations, enhanced by the interplay of spin, or- 
bital, charge and lattice degrees of freedom.^ These phenomena are a challenge for basic research 
and also bear big potentials for future applications as the related ground states are often accom- 
panied by so-called "colossal" effects, which are possible building blocks for tomorrow's correlated 
electronics. The measurement of the response of transition- metal oxides to ac electric fields is 
one of the most powerful techniques to provide detailed insight into the underlying physics that 
may comprise very different phenomena, e.g., charge order, molecular or polaronic relaxations, 
magnetocapacitance, hopping charge transport, ferroelectricity or density-wave formation. In the 
recent workj^ authors thoroughly discussed the mechanisms that can lead to colossal values of the 
dielectric constant, especially emphasizing effects generated by external and internal interfaces, 
including electronic phase separation. The authors of the work^Sl studied the materials showing 
so-called colossal dielectric constants (CDC), i.e. values of the real part of the permittivity e' ex- 
ceeding 1000. Since long, materials with high dielectric constants are in the focus of interest, not 
only for purely academic reasons but also because new high-e' materials are urgently sought after 
for the further development of modern electronics. In addition, authors of the work^ provided a 
detailed overview and discussion of the dielectric properties of CaCu^Ti^O\2 and related systems, 
which is today's most investigated material with colossal dielectric constant. Also a variety of 
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Table 2. Metallic and Semimetallic Elements 



Table 2: Metallic and Semimetallic Elements 





tii i mhpr 


ci \ prn pnfc 
ClClilCllljO 


alkali metals 


5 


Li, Na, K, Rb, Cs 


noble metals 


3 


Cu, Ag, Au 


polyvalent simple metals 


11 


Be, Mg, Zn, Cd, Hg, Al, Ga, In, Tl, Sn, Pb 


alkali-earth metals 


4 


Ca, Sr, Ba, Ra 


semi-metals 


4 


As, Sb, Bi, graphite 


transition metals 


23 


Fe, Ni, Co, etc. 


rare earths 


14 




actinides 


4 





further transition-metal oxides with large dielectric constants were treated in detail, among them 
the system Lai- x Sr x NiOt± where electronic phase separation may play a role in the generation 
of a colossal dielectric constant. In general, for the miniaturization of capacitive electronic ele- 
ments materials with high-e' are prerequisite. This is true not only for the common silicon-based 
integrated-circuit technique but also for stand-alone capacitors. 

Nevertheless, as regards to metals, the workable practical definition of Kittel can be adopted: 
metals are characterized by high electrical conductivity, so that a portion of electrons in metal 
must be free to move about. The electrons available to participate in the conductivity are called 
conduction electrons. Our picture of a metal, therefore, must be that it contains electrons which 
are free to move, and which may, when under the influence of an electric field, carry a current 
through the material. 

In summary, the 68 naturally occurring metallic and semimetallic elements^ can b e classified as 
it is shown in Table 2. 



3 Many-Particle Interacting Systems and Current operator 

Let us now consider a general system of ./V interacting electrons in a volume 0, described by the 
Hamiltonian 

N -q N 

H = (E £ + E (*)) + 2 E «(* " 3) =Ho + H 1 . (3.1) 

1=1 i=l vfij 

Here U(r) is a one-body potential, e.g. an externally applied potential like that due to the field 
of the ions in a solid, and v(fi — fj) is a two-body potential like the Coulomb potential between 
electrons. It is essential that U(r) and v(f{ — fj) do not depend on the velocities of the particles. 
It is convenient to introduce a "quantization" in a continuous space^^ via the operators \P^(r) 
and \P(r) which create and destroy a particle at r. In terms of and \P we have 

H = f d^\f)(^ + U(f))m (3-2) 
+ \ j j d 3 rd 3 r'^(f)^(f)v{f-f)^(f)^(f). 
Studies of flow problems lead to the continuity equatiorP^ES 

^ + VJ=0. (3.3) 
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This equation based on the concept of conservation of certain extensive variable. In nonequi- 
librium thermodynamics 42 the fundamental flow equations are obtained using successively mass, 
momentum, and energy as the relevant extensive variables. The analogous equation are known 
from electromagnetism. The central role plays a global conservation law of charge, q(t) = 0, for 
it refers to the total charge in a system. Charge is also conserved locally.^ This is described by 
Eq. (|3.3p . where n(r,t) and j are the charge and current densities, respectively. 
In quantum mechanics there is the connection of the wavefunction ip{f, t) to the particle mass- 
probability current distribution J 

J(r,t) = - W), (3-4) 

2mi 

where ip(r, t) satisfy the time-dependent Schrddinger equatiorP^^ 

ih^(r,t) = Hi>(r,t). (3.5) 

Consider the motion of a particle under the action of a time-independent force determined by a 
real potential V(r). Equation (|3.5p becomes 

It can be shown that for the probability density n(r, t) =ij}*ij) we have 

£+V/=0. (3.7) 

This is the equation of continuity and it is quite general for real potentials. The equation of 
continuity mathematically states the local conservation of particle mass probability in space. 
A thorough consideration of a current carried by a quasi-particle for a uniform gas of fermions, 
containing N particles in a volume £1, which was assumed to be very large, was performed within a 
semi-phenomenological theory of Fermi liquid.^ This theory describes the macroscopic properties 
of a system at zero temperature and requires knowledge of the ground state and the low-lying 
excited states. The current carried by the quasi-particle k is the sum of two terms: the current 
which is equal to the velocity Vk of the quasi-particle and the backflow of the medium.^ The 
precise definition of the current J in an arbitrary state \(p) within the Fermi liquid theory is given 
by 

j = <HE;>>' ( 3 - 8 ) 



where pi is the momentum of the iih particle and m its bare mass. To measure J it is necessary 
to use a reference frame moving with respect to the system with the uniform velocity hq/m. The 
Hamiltonian in the rest frame can be written 



tf = yj|J + V. (3.9) 



2m 

i 

It was assumed that V depends only on the positions and the relative velocities of the particles; 
it is not modified by a translation. In the moving system only the kinetic energy changes; the 
apparent Hamiltonian becomes 

H y iPi-H 2 +v = H _ mi yP ljrN ip4f (310) 
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Taking the average value of H q in the state \<p), and let E q be the energy of the system as seen 
from the moving reference frame, one find in the lim q — > 

where a refers to one of the three coordinates. This expression gives the definition of current in 
the framework of the Fermi liquid theory. For the particular case of a translationally invariant 
system the total current is a constant of the motion, which commutes with the interaction V and 
which, as a consequence, does not change when V is switched on adiabatically. For the particular 
state containing one quasi-particle k the total current is the same as for the ideal system 

(hk) 

Jk = — • 3.12 

m 

This result is a direct consequence of Galilean invariance. 
Let us consider now the many-particle Hamiltonian (|3.2j) 

H = H l +H 2 . (3.13) 

It will also be convenient to consider density of the particles in the following formP2l 

n(r) = 6(r — fj). 

i 

The Fourier transform of the particle density operator becomes 

i(q) = J d 3 r exp(— iqr) 5{f — fj) = exp(— iqfj). (3-14) 



i 



The particle mass-probability current distribution J in this "lattice" representation will take the 
form 

J(r) = n(f)v = - y^{— S(r — fj) + 8(r — ri) — } = (3.15) 
2 ^-^ m m 

i 

o yZi— exp(-iqri) + exp(-iqri) — }, 
2 m m 

i 

[ri,Pk] = ifrbik- 

Here v is the velocity operator. The direct calculation shows that 

[n(q),H] = exp(-iqTi) + exp(-iqfi) — } =qJ{q). (3.16) 



Thus the equation of motion for the particle density operator becomes 



dt h ' h 

or in another form 



' / " ( ' /) - l -[H,n(Q)} = -UM (3-17) 



dt 

which is the continuity equation considered above. Note, that 

[n(g),fl-i]_ = [n($,H 2 ]- = 0. 



dn(f) divJ(r), (3.18) 
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These relations holds in general for any periodic potential and interaction potential of the electrons 
which depend only on the coordinates of the electrons. 
It is easy to check the validity of the following relation 

[{n(q),H],n\q)] = [qJ{r),rJ {q)\ = (3.19) 

This formulae is the known f-sum rule^ which is a consequence from the continuity equation (for 
a more general point of view see RefP). 

Now consider the second-quantized Hamiltonian (13, 2j) . The particle density operator has the 
forrrJSKHID 

n(r) = e\p\r)\l/(r), n(q) = / d 3 r exp(— iqr)n(r). (3.20) 



Then we define 

j(r) = —(^W _ Wxp\). (3.21) 

2mi 

Here j is the probability current density, i.e. the probability flow per unit time per unit area 
perpendicular to j. The continuity equation will persist for this case too. Let us consider the 
equation motion 

= -I[n(iO,fli] - \W),H 2 ] = (3.22) 
eh 



2mi 



(^(r)V 2 V(r) - V 2 ^ f (r)tF(r)). 



Note, that [n(f),H 2 ] = 0. 
We find 

= -Vl(r). (3.23) 
Thus the continuity equation have the same form in both the "particle" and "field" versions. 



4 Tight-Binding and Modified Tight-Binding Approximation 

Electrons and phonons are the basic elementary excitations of a metallic solid. Their mutual 
interactions^'Sll'SSH manifest themselves in such observations as the temperature dependent re- 
sistivity and low-temperature superconductivity. In the quasiparticle picture, at the basis of this 
interaction is the individual electron-phonon scattering event, in which an electron is deflected 
in the dynamically distorted lattice. We consider here the scheme which is called the modified 
tight-binding interaction (MTBA). But firstly, we remind shortly the essence of the tight-binding 
approximation. The main purpose in using the tight-binding method is to simplify the theory suf- 
ficiently to make workable. The tight-binding approximation considers solid as a giant molecule. 



4.1 Tight-binding approximation 

The main problem of the electron theory of solids is to calculate the energy level spectrum of 
electrons moving in an ion lattice. 123351 The tight binding methocP™!] for energy band calcula- 
tions has generally been regarded as suitable primarily for obtaining a simple first approximation 
to a complex band structure. It was shown that the method should also be quite powerful in 
quantitative calculations from first principles for a wide variety of materials. An approximate 
treatment requires to obtain energy levels and electron wave functions for some suitable chosen 
one-particle potential (or pseudopotential), which is usually local. The standard molecular orbital 
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theories of band structure are founded on an independent particle model. 

As atoms are brought together to form a crystal lattice the sharp atomic levels broaden into bands. 
Provided there is no overlap between the bands, one expects to describe the crystal state by a 
Bloch function of the type, 

W = E Rn), (4.1) 

n 

where (j){r) is a free atom single electron wave function, for example such as Is and R n is the 
position of the atom in a rigid lattice. If the bands overlap or approach each other one should use 
instead of <j)(r) a combination of the wave functions corresponding to the levels in question, e.g. 
(a0(ls) + b(j)(2p)), etc. In the other words, this approach, first introduced to crystal calculation 
by F. Bloch, expresses the eigenstates of an electron in a perfect crystal in a linear combination of 
atomic orbitals and termed LCAO methodP^E^H 

Atomic orbitals are not the most suitable basis set due to the nonorthogonality problem. It 
was shown by many 

authorPHSSHSZl t h a t the very efficient basis set for the expansion (|4.ip is 
the atomic-like Wannier functions {w{r — i? n )}.^^HS3 These are the Fourier transforms of the 
extended Bloch functions and are defined as 

w(f-R n ) =AA- 1 /2^ e -^^ (f) . (4 . 2) 
k 

Wannier functions w(r — R n ) form a complete set of mutually orthogonal functions localized 
around each lattice site R n within any band or group of bands. They permit one to formulate 
an effective Hamiltonian for electrons in periodic potentials and span the space of a singly energy 
band. However, the real computation of Wannier functions in terms of sums over Bloch states is 
a complicated taskPSBHZI 

To define the Wannier functions more precisely let us consider the eigenfunctions ipj:(r) belonging 
to a particular simple band in a lattice with the one type of atom at a center of inversion. Let it 
satisfy the following equations with one-electron Hamiltonian H 

H^(r) = E(k)^(r), ^(f+ R n ) = e'^^r), (4.3) 

and the orthonormality relation (ip^\^g) = Srg where the integration is performed over the N unit 
cells in the crystal. The property of periodicity together with the property of the orthonormality 
lead to the orthonormality condition of the Wannier functions 

J d 3 rw*(r- R n )w(r- R m ) = 5 nm . (4.4) 

The set of the Wannier functions is complete, i.e. 

^2w*(r -Ri)w(f-Ri) = 6(r'-r). (4.5) 

i 

Thus it is possible to find the inversion of the Eq. (|4.2p which has the form 

^(r) = AT 1 / 2 e iMn w(r- R n ). (4.6) 

k 

These conditions are not sufficient to define the functions uniquely since the Bloch states ip^(r) 
are determined only within a multiplicative phase factor tp{k) according to 

w {f) = A/--V2 Y e^u % {r), (4.7) 
k 
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where ip(k) is any real function of k, and uAf) are Bloch functions.^ The phases tp(k) are usually 
chosen so as to localize w(r) about the origin. The usual choice of phase makes real and 

positive. This lead to the maximum possible value in w(0) and w{r) decaying exponentially away 
from r = 0. In addition, function tpt(r) with this choice will satisfy the symmetry properties 



It follows from the above consideration that the Wannier functions are real and symmetric, 



Analytical, three dimensional Wannier functions have been constructed from Bloch states formed 
from a lattice gaussians. Thus, in the condensed matter theory, the Wannier functions play an 
important role in the theoretical description of transition metals, their compounds and disordered 
alloys, impurities and imperfections, surfaces, etc. 

4.2 Interacting electrons on a lattice and the Hubbard model 

There are big difficulties in description of the complicated problems of electronic and magnetic 
properties of a metal with the d band electrons which are really neither "local" nor "itinerant" 
in a full sense. A better understanding of the electronic correlation effects in metallic systems 
can be achieved by the formulating of the suitable flexible model that could be used to analyze 
major aspects of both the insulating and metallic states of solids in which electronic correlations 
are important. 

The Hamiltonian of the interacting electrons with pair interaction in the second-quantized form 
is given by Eq. (|3.2p . Consider this Hamiltonian in the Bloch representation. We have 



w{f) = (w(r))* = w(—f). 




(4.8) 



A.- 



k 



Here (p a (k) is the Bloch function satisfying the equation 




(4.9) 




(4.10) 



k 



We find 




(4.11) 



mn klmn 



a ka + 



ka 




k 4 k 3 k 2 ki a Pt lv 
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Since the method of second quantization is based on the choice of suitable complete set of orthog- 
onal normalized wave functions, we take now the set {w\(f— R n )} of the Wannier functions. Here 
A is the band index. The field operators in the Wannier-function representation are given by 

^(r) = J2^x(r-iL)a n Xa, = E w *^~ ^) a l\a- ( 4 - 12 ) 

n n 

Thus we have 

4a. = ^ £ e~^a\^ a nX(7 = N~W £ e^a^. (4.13) 

k k 

Many of treatment of the correlation effects are effectively restricted to a non-degenerate band. 
The Wannier functions basis set is the background of the widely used Hubbard model. The 
Hubbard modeP^^I j S; i n a certain sense, an intermediate model (the narrow-band model) and 
takes into account the specific features of transition metals and their compounds by assuming 
that the d electrons form a band, but are subject to a strong Coulomb repulsion at one lattice 
site. The single-band Hubbard Hamiltonian is of the forrrp2MD 



H = Uj4a a j° + U / 2 Yl n i° n i-v ( 4 ' 14 ) 
ija 



Here a\ a and are the second-quantized operators of the creation and annihilation of the elec- 
trons in the lattice state w(f — Ri) with spin a. The Hamiltonian includes the intra-atomic 
Coulomb repulsion U and the one-electron hopping energy tij. The corresponding parameters of 
the Hubbard Hamiltonian are given by 



t ij = J d 3 rw*(r- Ri)H x {r)w(f- Rj), (4.15) 



= 2 

U= I I d A rd A r'w*{r - Ri)w*{r' - Ri) J - w(r f - Ri)w(r - R { ). (4.16) 
J J \f — r'\ 

The electron correlation forces electrons to localize in the atomic-like orbitals which are modelled 
here by a complete and orthogonal set of the Wannier wave functions {w(f — Rj)}- On the other 
hand, the kinetic energy is increased when electrons are delocalized. The band energy of Bloch 
electrons E(k) is defined as follows: 

^ = N' 1 E ( k ) eMik(Ri ~ Rj], (4.17) 

it 

where N is the number of lattice sites. The Pauli exclusion principle which does not allow two 
electrons of common spin to be at the same site, fifu — fiicr, plays a crucial role. Note, that 
the standard derivation of the Hubbard model presumes the rigid ion lattice with the rigidly 
fixed ion positions. We note that s-electrons are not explicitly taken into account in our model 
Hamiltonian. They can be, however, implicitly taken into account by screening effects and effective 
d-band occupation. 



4.3 Current operator for the tight-binding electrons 

Let us consider again a many-particle interacting systems on a lattice with the Hamiltonian (|4.11 j) . 
At this point, it is important to realize the fundamental difference between many-particle system 
which is uniform in space and many-particle system on a lattice. For the many-particle systems 
on a lattice the proper definition of current operator is a subtle problem. It was shown above that 
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a physically satisfactory definition of the current operator in the quantum many-body theory is 
given based upon the continuity equation. However, this point should be reconsidered carefully 
for the lattice fermions which are described by the Wannier functions. 

Let us remind once again that the Bloch and Wannier wave functions are related to each other 
by the unitary transformation of the form 

y? fe (r) = N-W^wi?- R n )exp[ikR n ], (4.18) 

Rn 

w(f-R n ) = N- 1 / 2 ^2^p k (f)eM-ikRn}- 
k 

The number occupation representation for a single-band case lead to 

V a (r) = w(f- RnWa, ^(f) = E W *^- Rn)aL- (4-19) 

n n 

In this representation the particle density operator and current density take the form 

"(f) = EE^( f " RiMr- Rj)a\ a a ja , (4.20) 

ij a 

^ = ^ E Yyir- RiYtwir - Rj) - Vw*(f- R^f- ^-)kW- 

ij a 

The equation of the motion for the particle density operator will consists of two contributions 

*?£L = -l[n(f), #i] - l -[n(r),H 2 ]. (4.21) 

The first contribution is 

[n(f),Hi] = E E F nm{r){tmiano a i<j ~ h n a\ a a ma ). (4.22) 

mni <J 

Here the notation was introduced 

F nm {f) =w*(r- R n )w{r- R m ). (4.23) 
In the Bloch representation for the particle density operator one finds 

Un4<T a ma), (4.24) 

mni a 

where 

F nm {k) = j d 3 rexp[-ikr\F nm (r) = J d 3 r exp[-ikr\w* (r - R n )w(f - R m ). (4.25) 
For the second contribution [n(f),H 2 ] we find 



'^^EEE^-^- (426) 



I'M'', ^ 

mn fst a a 1 



(^{mf^ls^a^a^ata'asa - (fm\H 2 \st)a ] m(jl a ] fa a ta ta sa 
Hi " s\H 2 \tn)a ] fa al al a t(T a na t - (fs\H 2 \nt)a ] fa a\ al a tu ta na ^. 
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For the single-band Hubbard Hamiltonian the last equation will take the form 

[n(f),H 2 ] = ^ E E F nm(r) — n n _ CT ). (4.27) 



The direct calculations give for the case of electrons on a lattice (e is a charge of an electron) 

dn{r) 



dt 

(h J2J2l w *(r~ Ri)V 2 w(f- R 3 ) - VV(f- Ri)w{r- R 3 )]a\ a a ]a 



(4.28) 



2mi 

ij a 

ij a 

Taking into account that 

divj(r) = (4.29) 
^ EEK( f - £*)VVr - R-) - vV(f - R\)w{r - Rj)]ala ja , 

ij a 

we find 

= -div/(f) - ieUj2 E #i(04r o i»("i-<r ~ n i-°)- ( 4 - 30 ) 

ij a 

This unusual result was analyzed critically by many authors. The proper definition of the current 
operator for the Hubbard model has been the subject of intensive discussions IIOUEII To clarify 
the situation let us consider the "total position operator" for our system of the electrons on a 
lattice 

N 

R = Y,Rj- ( 4 -31) 

In the "quantized" picture it has the form 

R = Y^J d 3 r&\r)Rj&(r) (4.32) 

3 

= EEE/ d 3 rR jW *(f- R m )w(f- R^al, 



j mn 



j m /j, 

where we took into account the relation 

-;3 „,„,,* 



We find that 



d 6 rw*(f - R m )w{r- R n ) = 6 mn . (4.33) 

m 
m 

[R,al a a ia ]- = 0. 
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Let us consider the local particle density operator nj CT = a\ a ai a . 



l -[n ia ,H}- = ^2Uj{a\ a a ja - a^ a a ia ). (4.35) 



dt ft 

j 



It is clear that the current operator should be defined on the basis of the equation 



-i 



j = e(—)[R,H]„. (4.36) 
Defining the so-called polarization operato J^ 103 *^^^ 



m a 



we find the current operator in the form 

j = V = e (-^) ~~ Rn)t m n^ ma a na . (4.38) 

mn a 

This expression of the current operator is a suitable formulae for studying of the transport prop- 
erties of the systems of correlated electron on a lattice The consideration carried out in 
this section demonstrate explicitly the specific features of the many-particle interacting systems 
on a lattice. 



4.4 Electron-lattice interaction in metals 

In order to understand quantitatively the electrical, thermal and superconducting properties of 
metals and their alloys one needs a proper description an electron-lattice interaction.^ In the 
physics of molecules^^ the concept of an intermolecular force requires that an effective separation 
of the nuclear and electronic motion can be made. This separation is achieved in the Born- 
Oppenheimer approximationP^lH^ Closely related to the validity of the Born-Oppenheimer ap- 
proximation is the notion of adiabaticity. The adiabatic approximation is applicable if the nuclei 
is much slower than the electrons. The Born-Oppenheimer approximation consists of separating 
the nuclear motion and in computing only the electronic wave functions and energies for fixed 
position of the nuclei. In the mathematical formulation of this approximation, the total wave 
function is assumed in the form of a product both of whose factors can be computed as solutions 
of two separate Schrodinger equations. In most applications the separation is valid with suffi- 
cient accuracy, and the adiabatic approach is reasonable, especially if the electronic properties of 
molecules are concerned. 

The conventional physical picture of a metal adopts these ideas^"^ and assumes that the elec- 
trons and ions are essentially decoupled from one another with an error which involves the small 
parameter m/M, the ratio between the masses of the electron and the ion. The qualitative ar- 
guments for this statement are the following estimations. The maximum lattice frequency is of 
the order 10 13 sec _1 and is quite small compared with a typical atomic frequency. This latter 
frequency is of order of 10 15 sec -1 . If the electrons are able to respond in times of the order of 
atomic times then they will effectively be following the motion of the lattice instantaneously at 
all frequencies of vibration. In other words the electronic motion will be essentially adiabatic. 
This means that the wave functions of the electrons adjusting instantaneously to the motion of 
the ions. It is intuitively clear that the electrons would try to follow the motion of the ions in 
such a way as to keep the system locally electrically neutral. In other words, it is expected that 
the electrons will try to respond to the motion of the ions in such a way as to screen out the local 
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charge fluctuations. 

The construction of an electron-phonon interaction requires the separation of the Hamiltonian de- 
scribing mutually interacting electrons and ions into terms representing electronic quasiparticles, 
phonons, and a residual interaction!^ 33 ^ 52 ^ 86 ™ 88 ! For the simple metals the interaction between 
the electrons and the ions can be described within the pseudopotential method or the muffin-tin 
approximation. These methods could not handled well the d bands in the transition metals. They 
are too narrow to be approximated as free-electron-like bands but too broad to be described as core 
ion states. The electron-phonon interaction in solid is usually described by the Frohlich Hamilto- 
nianP*EJ3 We consider below the main ideas and approximations concerning to the derivation of 
the explicit form of the electron-phonon interaction operator. 

Consider the total Hamiltonian for the electrons with coordinates fl and the ions with coordi- 
nates R m , with the electron cores which can be regarded as tightly bound to the nuclei. The 
Hamiltonian of the N ions is 



fj2 N 2 ZN ZN 2 



m=l t=l *J=1 

N 

Vi{R m -R n )+Y^ Uie{n; Rm). 



n>m m=l 



Each ion is assumed to contribute Z conduction electrons with coordinates fl (i = 1, . . . , ZN). 
The first two terms in Eq. (j4.39p are the kinetic energies of the electrons and the ions. The third 
term is the direct electron-electron Coulomb interaction between the conduction electrons. The 
next two terms are short for the potential energy for direct ion-ion interaction and the potential 
energy of the ZN conduction electrons moving in the field from the nuclei and the ion core 
electrons, when the ions take instantaneous position R m (m = 1, . . . , N). The term Vi{R m — R n ) 
is the interaction potential of the ions with each other, while Ui e (fl; R m ) represents the interaction 
between an electron at fl and an ion at R m . Thus the total Hamiltonian of the system can be 
represented as the sum of an electronic and ionic part. 



where 



H = H e + H h (4.40) 

f-2 ZN i ZN 2 N 

H * = £ V 'l + 2 ^ T^—r\ + E U ie (r7,R m ), (4.41) 

i=l i,j'=l ' * m=l 

and 

h 2 N 

m=l n>m 

The Schrodinger equation for the electrons in the presence of fixed ions is 

H e V(K, R,f) = E(K, R)V(K, R, r), (4.43) 

in which K is the total wave vector of the system, R and f denote the set of all electronic and 
ionic coordinates. It is seen that the energy of the electronic system and the wave function of 
the electronic state depend on the ionic positions. The total wave function for the entire system 
of electrons plus ions <fr(Q,R, r) can be expanded, in principle, with respect to the as basis 
functions 

$(Q, R,f) = J2 L (Q, K, R)^(K, R,r). (4.44) 

K 
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We start with the approach which uses a fixed set of basis states. Let us suppose that the ions of 
the crystal lattice vibrate around their equilibrium positions with a small amplitude, namely 
R m = R^ + u m , where u m is the deviation from the equilibrium position R^. Let us consider 
an idealized system in which the ions are fixed in these positions. Suppose that the energy bands 
E n (k) and wave functions ip n (k,r) are known. As a result of the oscillations of the ions, the 
actual crystal potential differs from that of the rigid lattice. This difference is possible to treat as 
a perturbation. This is the Bloch formulation of the electron-phonon interaction. 
To proceed we must expand the potential energy V(f — R) of an electron at f in the field of an 
ion at R m in the atomic displacement u m 

V{f- R m ) ~ V(f- R°J - u m VV(r- R°J + . . . (4.45) 

The perturbation potential, including all atoms in the crystal, is 

V = -J2^mVV(f-R° m ). (4.46) 

m 

This perturbation will produce transitions between one-electron states with the corresponding 
matrix element of the form 



M mKm = J r m {Kr)V^Mf)d 6 r. (4.47) 

To describe properly the lattice subsystem let us remind that the normal coordinate is defined 
by the relationPMl 

(R m - R\) = u m = {h/2MM) 1 ' 2 Y, Q^M exp(ig7&), (4.48) 

q,u 

where N is the number of unit cells per unit volume and e u (q) is the polarization vector of the 

phonon. The Hamiltonian of the phonon subsystem in terms of normal coordinates is written 
af P[86] 

where ji denote polarization direction and the q summation is restricted to the Brillouin zone 
denoted as BZ. It is convenient to express u m in terms of the second quantized phonon operators 

u m = (h/^M) 1 / 2 Y.lU'H^r^M^Mm^u + expHgifoj^J, (4.50) 

in which v denotes a branch of the phonon spectrum, e u (q) is the eigenvector for a vibrational 
state of wave vector q and branch u, and bt u (b^ u ) is a phonon creation (annihilation) operator. 
The matrix element M m k )n q becomes 



M mKnq = -{h/2MMf/ 2 Y,{^{k - q)A mn (k, q)[co u (k - q)]' 1 ' 2 '(b^ + b^)) . (4.51) 
Here the quantity A mn is given by 



A mn (k, q)=M r m (k, ?)W(r)^ n {q, r)d\. (4.52) 



It is well knowrP^ISS] that there is the distinction between normal processes in which vector (k — q) 
is inside the Brillouin zone and Umklapp processes in which vector (k — q) must be brought back 
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into the zone by addition of a reciprocal lattice vector G. 

The standard simplification in the theory of metals consists of replacement of the Bloch functions 
i'niQ,^) by the plane waves 

in which V is the volume of the system. With this simplification we get 

A mn (k,q)=i(k-q)V((k-q)). (4.53) 



Introducing the field operators ip(r),2jj' (r) and the fermion second quantized creation and anni- 

;tron of wave vector / 



hilation operators ->, a n ^ for an electron of wave vector k in band n in the plane wave basis 



nk 

qn 



and the set of quantities 

/ . \l/2 . 

r mTM /(«, q) = - (h/2Muj u {k - q)j e v (k - q)A mn (k,q), 
we can write an interaction Hamiltonian for the electron-phonon system in the form 

H ei = A/" 1/2 E q) (^«^- ? > + «^U>) • ( 4 - 54 ) 

nlu 

This Hamiltonian describes the processes of phonon absorption or emission by an electron in the 
lattice, which were first considered by Bloch. Thus the electron-phonon interaction is essentially 
dynamical and affects the physical properties of metals in a characteristic way. 
It is possible to show^ that in the Bloch momentum representation the Hamiltonian of a system 
of conduction electrons in metal interacting with phonons will have the form 



where 



H = H e + H l + H ei , (4.55) 

H e = J2E(p)aiap, (4.56) 
v 

.. |g|<9m 

F - = E E r^P-P^Mkv + bl^)- (4.58) 

" p'=p+q+G 

The Frohlich model ignores the Umklapp processes (G ^ 0) and transverse phonons and takes the 
unperturbed electron and phonon energies as 

E(p) = - E F , uj(q) = v° s q, (q < q m ). 
Here v s is the sound velocity of the free phonon. The other notation are: 

\q\ =q, \p\= p, q m = (67r 2 n;) 1/3 - 
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Thus we obtain 

H ei = ^2v(q)al +q a p (b q + bt q ), (4.59) 
p,q 

where v(q) is the Fourier component of the interaction potential 

v(q) = g{uj{q)/2) l /\ g = [2E F /3Afn^ s 2 ] 1/2 • 

Here ?ij is the ionic density. The point we should like to emphasize in the present context is that 
the derivation of this Hamiltonian is based essentially on the plane wave representation for the 
electron wave function. 

4.5 Modified tight-binding approximation 

Particular properties of the transition metals, their alloys and compounds follow, to a great 
extent, from the dominant role of d-electrons. The natural approach to description of electron- 
lattice effects in such type of materials is the modified tight-binding approximation (MTBA). 
The electron-phonon matrix element in the Bloch picture is taken between electronic states of 
the undeformed lattice. For transition metals it is not easy task to estimate the electron-lattice 
interaction matrix element due to the anisotropy and other f actor s IHHHl2i] There is an alterna- 
tive description, introduced by Fr6hlichp22HI2H anc i w hich was termed the modified tight-binding 
approximation (MTBA). In this approach the electrons are moving adiabatically with the ions. 
Moreover, the coupling of the electron to the displacement of the ion nearest to it, is already 
included in zero order of approximation. This is the basis of modified tight-binding calculations 
of the electron-phonon interaction which purports to remove certain difficulties of the conven- 
tional Bloch tight-binding approximation for electrons in narrow band. The standard Hubbard 
Hamiltonian should be rederived in this approach in terms of the new basis wave functions for the 
vibrating lattice. This was carried out by Barisic, Labbe and FriedelP^S They derived a model 
Hamiltonian which is a generalization of the single-band Hubbard model^ including the lattice 
vibrations. The hopping integral tij of the single-band Hubbard model (|4.14j) is given by 

Uj = J d 3 rw*(f- Rj) + v 'f(f- w (?- Ri)- (4.60) 

Here we assumed that V s t is a short-range, self-consistent potential of the lattice suitable screened 
by outer electrons. Considering small vibrations of ions we replace in Eq. (|4.60p the ion position 
Ri by (R® + Ui) , i.e. its equilibrium position plus displacement. The unperturbed electronic wave 
functions must be written as a Bloch sum of displaced and suitable (approximately) orthonormal- 
ized atomic-like functions 

d 3 rw*(r - R® - Uj)w{r- R% - Hi) « kj- (4.61) 

As it follows from Eq. (|4.6ip . the creation and annihilation operators a^ a ,a\ trT may be introduced 
in the deformed lattice so as to take partly into account the adiabatic follow up of the electron 
upon the vibration of the lattice. The Hubbard Hamiltonian Eq. (|4.14|) can be rewritten in the 

H = t J2hi*+Yl + ^ ~ - ^kW + U / 2 Yl n ™ n *-°- ( 4 - 62 ) 

For small displacements u^, we may expand t(R) as 

t(8j + u 3 - Rf - ui) « t(i?° - Rj) + ^-l^.goiu, -*) + ■■■ ( 4 -63) 

oR 3 1 
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Table 3: Slater coefficients 

Table 3. Slater coefficients 





element 


element 


element 


element 


element 


element 


q = 0.93 


Ti 


V 


Cr 


Mn 


Fe 


Co 


q = 0.91 


Zr 


Nb 


Mo 


Tc 


Ru 


Rh 


q = 0.87 


Hf 


Ta 


W 


Re 


Os 


Ir 



Using the character of the exponential decrease of the Slater and Wannier functions the following 
approximation may be usecP^SHEH 

^5> = - 9 „4-*(«). (4.64) 
dR \R\ 

Here qo is the Slater coefficient!^] originated in the exponential decrease of the wave functions 
of (i-electrons; q^ 1 related to the range of the d function and is of the order of the interatomic 
distance. The Slater coefficients for various metals are tabulatedP^S The typical values are given 
in Table 3. 

It is of use to rewrite the total model Hamiltonian of transition metal H = H e + Hi + H e i in the 
quasi-momentum representation. We have 

H e = Y, E{k)a\ a a kc + U/2N £ 

4it afc 2t a l 3 ; a fc44- 5 (^i -h + h-h + G). (4.65) 

ka kxkik-ik^G 

For the tight-binding electrons in crystals we use E(k) = 2 ^ Q t(a a ) cos(k a a a ), where t(a) is the 
hopping integral between nearest neighbours, and a a (a = x,y.z) denotes the lattice vectors in a 
simple lattice with an inversion center. 
The electron-phonon interaction is rewritten as 

H ™ = E E E dlkA^Ubl + b- qu )5(h -k + q + G), (4.66) 

fcfei qG vo 



( 1 \ 1/2 

&1 = \(MMujy(k))) (467) 
a a e v {k\ 



where 



F kkl = 2iq t{d a ) Z\ ' (sin(a a fc) - sin((a J x )j , (4.68) 

where J\f is the number of unit cells in the crystal and M is the ion mass. The e„ (q) are the polariza- 
tion vectors of the phonon modes. Operators b\ u and b qu are the creation and annihilation phonon 
operators and co u (k) are the acoustical phonon frequencies. Thus we can 

describ e™[i2Zl[T29j[i30| the 

transition metal by the one-band model which takes into consideration the electron-electron and 
electron-lattice interaction in the framework of the MTBA. It is possible to rewrite (|4.66p in the 
following fornPSEa 

H * = E E V "& * + ^)Q^ a t+qa a ^, (4-69) 
va kq 

where 

V(k,k + q) = T^yi/aE^X® (^J-sma a (k-q)) . (4.70) 
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The one-electron hopping t(a a ) is the overlap integral between a given site R m and one of the two 
nearby sites lying on the lattice axis a a . For the ion subsystem we have 

H i = \H ( P £ P «» + QtvQi" = E^X^V + 1/2), (4.71) 

qv qu 

where P qv and Q qu are the normal coordinates. Thus, as in the Hubbard modelp^ the d- and 
s(p)-bands are replaced by one effective band in our model. However, the s-electrons give rise 
to screening effects and are taken into effects by choosing proper value of U and the acoustical 
phonon frequencies. It was shown by Ashkenazi, Dacorogna and Peter^^^l that the MTBA 
approach for calculating electron-phonon coupling constant based on wave functions moving with 
the vibrating atoms lead to same physical results as the Bloch approach within the harmonic 
approximation. For transition metals and narrow band compounds the MTBA approach seems 
to be yielding more accurate results, especially in predicting anisotropic properties. 

5 Charge and Heat Transport 

We now tackle the transport problem in a qualitative fashion. This crude picture has many obvious 
shortcomings. Nevertheless, the qualitative description of conductivity is instructive. Guided 
by this instruction the results of the more advanced and careful calculations of the transport 
coefficients will be reviewed below in the next sections. 

5.1 Electrical resistivity and Ohm law 

Ohm law is one of the equations used in the analysis of electrical circuits. When a steady current 
flow through a metallic wire, Ohm law tell us that an electric field exists in the circuit, that like 
the current this field is directed along the uniform wire, and that its magnitude is J/cr, where J 
is the current density and a the conductivity of the conducting material. Ohm law states that, 
in an electrical circuit, the current passing through most materials is directly proportional to 
the potential difference applied across them. A voltage source, V, drives an electric current, I , 
through resistor, R, the three quantities obeying Ohm law: V = IR. 

In other terms, this is written often as: I = V/R, where I is the current, V is the potential 
difference, and R is a proportionality constant called the resistance. The potential difference is 
also known as the voltage drop, and is sometimes denoted by E or U instead of V. The SI unit 
of current is the ampere; that of potential difference is the volt; and that of resistance is the ohm, 
equal to one volt per ampere. The law is named after the physicist Georg Ohm, who formulated 
it in 1826 . The continuum form of Ohm's law is often of use 

J = a-E, (5.1) 

where J is the current density (current per unit area) , a is the conductivity (which can be a tensor 
in anisotropic materials) and E is the electric field . The common form V = I ■ R used in circuit 
design is the macroscopic, averaged-out version. The continuum form of the equation is only valid 
in the reference frame of the conducting material. 

A conductor may be defined as a material within which there are free charges, that is, charges that 
are free to move when a force is exerted on them by an electric field. Many conducting materials, 
mainly the metals, show a linear dependence of I on V. The essence of Ohm law is this linear 
relationship. The important problem is the applicability of Ohm law. The relation R ■ I = W is 
the generalized form of Ohm law for the current flowing through the system from terminal A to 
terminal B. Here I is a steady dc current, which is zero if the work W done per unit charge is zero, 
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while I 7^ or W 7^ 0. If the current in not too large, the current I must be simply proportional 
to W. Hence one can write R - 1 = W, where the proportionality constant is called the resistance 
of the two-terminal system. The basic equations are: 

V x E = Ann, (5.2) 

Gauss law, and 

^ + VxJ = 0, (5.3) 

charge conservation law. Here n is the number density of charge carriers in the system. Equations 
(|5,2p and (|5.3p are fundamental. The Ohm law is not. However, in the absence of nonlocal effects, 
Eq. (|5.ip is still valid. In an electric conductor with finite cross section it must be possible a surface 
conditions on the current density J. Ohm law does not permit this and cannot, therefore, be quite 
correct. It has to be supplemented by terms describing a viscous flow. Ohm law is a statement of 
the behavior of many, but not all conducting bodies, and in this sense should be looked upon as 
describing a special property of certain materials and not a general property of all matter. 



5.2 Drude-Lorentz model 

The phenomenological picture described above requires the microscopic justification. We are 
concerned in this paper with the transport of electric charge and heat by the electrons in a 
solid. When our sample is in uniform thermal equilibrium the distribution of electrons over 
the eigenstates available to them in each region of the sample is described by the Fermi-Dirac 
distribution function and the electric and heat current densities both vanish everywhere. Non- 
vanishing macroscopic current densities arise whenever the equilibrium is made non-uniform by 
varying either the electrochemical potential or the temperature from point to point in the sample. 
The electron distribution in each region of the crystal is then perturbed because electrons move 
from filled states to adjacent empty states. 

The electrical conductivity of a material is determined by the mobile carriers and is proportional 
to the number density of charge carriers in the system, denoted by n, and their mobility, fx, 
according to 

a ~ nefM. (5-4) 

Only in metallic systems the number density of charge carriers is large enough to make the 
electrical conductivity sufficiently large. The precise conditions under which one substance has a 
large conductivity and another substance has low ones are determined by the microscopic physical 
properties of the system such as energy band structure, carrier effective mass, carrier mobility, 
lattice properties, and the presence of impurities and imperfections. 

Theoretical considerations of the electric conductivity were started by P. Drude within the classical 
picture about hundred years ago.^1^31 jj e p U ^ forward a free electron model that assumes a 
relaxation of the independent charge carriers due to driving forces (frictional force and the electric 
field). The current density was written as 

j = V^-Et. (5.5) 
m 

Here r is the average time between collisions, E is the electric field, m and e are the mass and 
the charge of the electron. The electric conductivity in the Drude modeP^ is given by 

a = (5.6) 
m 
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The time r is called the mean lifetime or electron relaxation time. Then the Ohm's law can be 
expressed as the linear relation between current density J and electric field E 

J = cjE. (5.7) 

The electrical resistivity R of the material is equal to 

R=j. (5.8) 

The free-electron model of Drude is the limiting case of the total derealization of the outer 
atomic electrons in a metal. The former valence electrons became conduction electrons. They 
move independently through the entire body of the metal; the ion cores are totaly ignored. The 
theory of Drude was refined by Lorentz. Drude-Lorentz theory assumed that the free conduction 
electrons formed an electron gas and were impeded in their motion through crystal by collisions 
with the ions of the lattice. In this approach, the number of free electrons n and the collision time 
r, related to the mean free path 77 = 2tv and the mean velocity v, are still adjustable parameters. 
Contrary to this, in the Bloch model for the electronic structure of a crystal, though each valence 
electron is treated as an independent particle, it is recognized that the presence of the ion cores 
and the other valence electrons modifies the motion of that valence electron. 
In spite of its simplicity, Drude model contains some delicate points. Each electron changes its 
direction of propagation with an average period of 2r. This change of propagation direction is 
mainly due to a collision of an electron with an impurity or defect and the interaction of electron 
with a lattice vibration. In an essence, t is the average time of the electron motion to the first 
collision. Moreover, it is assumed that the electron forgets its history on each collision, etc. To 
clarify these points let us consider the notion of the electron drift velocity. The electrons which 
contribute to the conductivity have large velocities, that is large compared to the drift velocity 
which is due to the electric field, because they are at the top of the Fermi surface and very 
energetic. The drift velocity of the carriers is intimately connected with the collision time r 

V d = OCT, 

where a is a constant acceleration between collision of the charge carriers. In general, the mean 
drift velocity of a particle over iV free path is 

v d ~ \a[r + (At) 2 /r] . 

This expression shows that the drift velocity depends not only on the average value r but also on 
the standard deviation (At) of the distribution of times between collisions. An analysis shows that 
the times between collisions have an exponential probability distribution. For such a distribution, 
At = t and one obtains Vd = olt and J = ne 2 jvnEr. Assuming that the time between collisions 
always has the same value r we find that (At) = and Vd = \olt and J = ne 2 /2mEr. 
The equations (|5.7p and (|5.8p are the most fundamental formulas in the physics of electron con- 
duction. Note, that resistivity is not zero even at absolute zero, but is equal to the so called 
" residual resistivity" . For most typical cases it reasonably to assume that scattering by impurities 
or defects and scattering by lattice vibrations are independent events. As a result, the relation 
(|5.6p will take place. There is a big variety (and irregularity) of the resistivity values for the 
elements not speaking on the huge variety of substances and materials P^^L^ 
In a metal with spherical Fermi surface in the presence of an electric field E, the Fermi surface 
would affect a Ak displacement, Ak = k — ko. The simplest approximation is to suppose a rigid 
displacement of the Fermi sphere with a single relaxation time r, 

, dk . (k — ko) f± „ s 

h— + h- — = eE. (5.9) 

dt t 
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Thus we will have at equilibrium 



er -> 



Ak = —E. (5.10) 
n 



The corresponding current density will take the form 



We get from Eq.([5 



J = } nk <**>>> = W? j Sk e * AkSS *°- < 5J1) 

Let us consider briefly the frequency dependence of a. Consider a gas of noninteracting electrons 
of number density n and collision time r. At low frequencies collisions occur so frequently that 
the charge carriers are moving as if within a viscous medium, whereas at high frequencies the 
charge carriers behave as if they were free. These two frequency regimes are well-known in the 
transverse electromagnetic response of metalsP^^DE^^I] The electromagnetic energy given 
to the electrons is lost in collisions with the lattice, which is the "viscous medium". The relevant 
frequencies in this case satisfy the condition lot <C 1. Thus in a phenomenological description^^ 
one should introduce a conductivity a and viscosity r\ by 

e 2 f c n 1 — ~- _ 

a = , r\ = -rav l nT c . (5.13) 

m 2 

On the other hand, for lot ^> 1 viscous effects are negligible, and the electrons behave as the nearly 
free particles. For optical frequencies they can move quickly enough to screen out the applied field. 
Thus, two different physical mechanisms are suitable in the different regimes defined by lot <C 1 
and lot S> 1. 

In a metal impurity atoms and phonons determine the scattering processes of the conduction 
electrons. The electrical force on the electrons is eE. The "viscous" drag force is given by —ttiv/t. 
Then one can write the equation 

TflV 

mv = eE . (5-14) 

r 

For E ~ exp(— iuit), the oscillating component of the current is given by 

J(lo) = nev(uo) = g{uj)E{uS), (5.15) 

where 

°"M = i : > °0 = • (5.16) 

1 — ilot m 

For low frequencies we may approximate Eq. (|5.14p as v ~ (er/m)E. For high frequencies we may 
neglect the collision term, so v ~ (e/m)E. Thus the behavior of the conductivity as a function of 
frequency can be described on the basis of the formula Eq. (|5.16p . 

Let us remark on a residual resistivity, i.e. the resistivity at absolute zero. Since real crystals 
always contain impurities and defects the resistivity is not equal zero even at absolute zero. If one 
assume that the scattering of a wave caused by impurities (or defects) and by lattice vibrations are 
independent events, then the total probability for scattering will be the sum of the two individual 
probabilities. The scattering probability is proportional to 1/r, where r is the mean lifetime 
or relaxation time of the electron motion. Denoting by 1/n the scattering probability due to 
impurities and defects and by I/V2 the scattering probability due to lattice vibrations we obtain 
for total probability the equality 

1/r = 1/n + l/r 2 ; 1/(7 = 1/aa + l/a 2 . (5.17) 
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This relation is called Matthiessen rule. In practice, this relation is not fulfilled well (see RefsP^SES) 
The main reason for the violation of the Matthiessen rule are the interference effects between 
phonon and impurity contributions to the resistivity. RefsP^M^Q] gj ve a comprehensive review of 
the subject of deviation from Matthiessen rule and detailed critical evaluation of both theory and 
experimental data. 

5.3 The low- and high-temperature dependence of conductivity 

One of the most informative and fundamental properties of a metal is the behavior of its elec- 
trical resistivity as a function of temperature. The temperature dependence of the resistivity is 
a good indicator of important scattering mechanisms for the conduction electrons. It can also 
suggests in a general way what the solid-state electronic structure is like. There are two limiting 
cases, namely, the low temperature dependence of the resistivity for the case when T < 0£>, where 
Od is effective Debye temperature, and the high temperature dependence of the resistivity, when 

t > e D . 

The electrical resistivity of metals is due to two mechanisms, namely, (i) scattering of electrons on 
impurities (static imperfections in the lattice), and (ii) scattering of electrons by phonons. Sim- 
plified treatment assumes that one scattering process is not influenced by the other (Matthiessen 
rule). The first process is usually temperature independent. For a typical metal the electrical 
resistivity R(T), as a function of the absolute temperature T, can be written as 

R(T) = R + Ri(T), (5.18) 

where i?o is the residual electrical resistivity independent of T, and Ri(T) is the temperature- 
dependent intrinsic resistivity. The quantity i?o is due to the scattering of electrons from chemical 
and structural imperfections. The term Ri(T) is assumed to result from the interaction of elec- 
trons with other degrees of freedom of a crystal. In general, for the temperature dependence of the 
resistivity three scattering mechanisms are essential, (i) electron-phonon scattering, (ii) electron- 
magnon scattering and (iii) electron-electron scattering. The first one gives T 5 or T 3 dependence 
at low temperaturesP The second one, the magnon scattering is essential for the transition metals 
because some of them show ferromagnetic and antiferromagnetic properties.^ This mechanism 
can give different temperature dependence due to the complicated (anisotropic) dispersion of the 
magnons in various structures. The third mechanism, the electron-electron scattering is responsi- 
ble for the R ~ T 2 dependence of resistivity. 

Usually, the temperature-dependent electrical resistivity is tried to fit to an expression of the form 

R(T) = R + Ri(T) = R + AT 5 + BT 2 + (CT 3 ) + . . . (5.19) 

This dependence corresponds to Mathiessen rule, where the different terms are produced by dif- 
ferent scattering mechanisms. The early approach for studying of the temperature variation of 
the conductivity^^M^ was carried out by Sommerfeld, Bloch and Houston. Houston explained 
the temperature variation of conductivity applying the wave mechanics and assuming that the 
wave-lengths of the electrons were in most cases long compared with the interatomic distance. 
He then solves the Boltzmann equation, using for the collision term an expression taken from the 
work of Debye and Waller on the thermal scattering of x-rays. He obtained an expression for the 
conductivity as a function of a mean free path, which can be determined in terms of the scattering 
of the electrons by the thermal vibrations of the lattice. Houston found a resistance proportional 
to the temperature at high temperatures and to the square of temperature at low temperatures. 
The model used by Houston for the electrons in a metal was that of Sommerfeld - an ideal gas 
in a structureless potential well. Bloch improved this approach by taking the periodic structure 
of the lattice into account. For the resistance law at low temperatures both Houston and Bloch 
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results were incorrect. Houston realized that the various treatments of the mean free path would 
give different variations of resistance with temperature. In his later workP^ he also realized that 
the Debye theory of scattering was inadequate at low temperatures. He applied the Brillouin 
theory of scattering and arrived at T 5 law for the resistivity at low temperatures and T at high 
temperatures. Later on, it was shown by many authors^ that the distribution function obtained 
in the steady state under the action of an electric field and the phonon collisions does indeed lead 
to R ~ T 5 . The calculations of the electron-phonon scattering contribution to the resistivity by 
BloctpZl anc j GruneiseEpSl i eac i to the following expression 

t 5 r e/T z b 
R ( T >~WL - V-im -«-■) ' <5 ' 20) 

which is known as the Bloch-Gruneisen law. 

A lot of efforts has been devoted to the theory of transport processes in simple met als KllIMIlMil 
such as the alkali metals. The Fermi surface of these metals is nearly spherical, so that band- 
structure effects can be either neglected or treated in some simple approximation. The effect of 
the electron-electron interaction in these systems is not very substantial. Most of the scattering is 
due to impurities and phonons. It is expected that the characteristic T 2 dependence of electron- 
electron interaction effects can only be seen at very low temperatures T, where phonon scattering 
contributes a negligible T 5 term. In the non-simple metallic conductors, and in transition metals, 
the Fermi surfaces are usually far from being isotropic. Moreover, it can be viewed as the two- 
component systems^!] where one carrier is an electron and the other is an inequivalent electron 
(as in s — d scattering) or a hole. It was shown that anisotropy such as that arising from a 
nonspherical Fermi surface or from anisotropic scattering can yield a T 2 term in the resistivity at 
low temperatures, due to the deviations from Mathiessen rule. This term disappears at sufficiently 
high T. The electron-electron Umklapp scattering contributes a T 2 term even at high T. It was 
conjectured (see RefP^Zl) that the effective electron-electron interaction due to the exchange 
of phonons should contribute to the electrical resistivity in exactly the same way as the direct 
Coulomb interaction, namely, giving rise to a T 2 term in the resistivity at low temperatures. 
The estimations of this contribution show^^Sl that it can alter substantially the coefficient of 
the T 2 term in the resistivity of simple and polyvalent metals. The role of electron-electron 
scattering in transition metals was discussed in Refs PUSHED a calculation of the electrical and 
thermal resistivity of Nb and Pd due to electron-phonon scattering was discussed in Refill A 
detailed investigation^ 3 -^ of the temperature dependence of the resistivity of Nb and Pd showed 
that a simple power law fit cannot reconcile the experimentally observed behavior of the transition 
metals. Matthiessen rule breaks down and simple Bloch-Gruneisen theory is inadequate to account 
for the experimental data. In particular, in Ref. 153 it has been shown that the resistivity of Pd 
can be expressed by a T 2 function where, on the other hand, the temperature dependence of the 
resistivity of Nb should be represented by a function of T more complicated than the T 3 . It seems 
to be plausible that the low-temperature behavior of the resistivity of transition metals may be 
described by a rational function of (AT 5 + BT 2 ). This conjecture will be justified in section [TT1 
For real metallic systems the precise measurements show a quite complicated picture in which 
the term Ri(T) will not necessarily be proportional to T 5 for every metal (for detailed review see 
RefsP™^. The purity of the samples and size-effect contributions and other experimental 
limitations can lead to the deviations from the T 5 law. There are a lot of other reasons for such 
a deviation. First, the electronic structures of various pure metals differ very considerably. For 
example, the Fermi surface of sodium is nearly close to the spherical one, but those of transition 
and rare-earth metals are much more complicated, having groups of electrons of very different 
velocities. The phonon spectra are also different for different metals. It is possible to formulate 
that the T 5 law can be justified for a metal of a spherical Fermi surface and for a Debye phonon 
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spectrum. Moreover, the additional assumptions are an assumption that the electron and phonon 
systems are separately in equilibrium so that only one phonon is annihilated or created in an 
electron-phonon collision, that the Umklapp processes can be neglected, and an assumption of a 
constant volume at any temperature. Whenever these conditions are not satisfied in principle, 
deviations from T 5 law can be expected. This takes place, for example, in transition metals as a 
result of the s — d transitions^^ due to the scattering of s electrons by phonons. This process can 
be approximately described as being proportional to T 7 with 7 somewhere between 5 and 3. The 
s — d model of electronic transport in transition metals was developed by MottP 531531 In this 
model, the motion of the electrons is assumed to take place in the nearly-free-electron-like s-band 
conduction states. These electrons are then assumed to be scattered into the localized d states. 
Owing to the large differences in the effective masses of the s and d bands, large resistivity result. 
In RefP^ the temperature of the normal-state electrical resistivity of very pure niobium was 
reported. The measurements were carried out in the temperature range from the superconducting 
transition (T c = 9.25K) to 300-FT in zero magnetic field. The resistance-versus-temperature data 
were analyzed in terms of the possible scattering mechanisms likely to occur in niobium. To fit 
the data a single-band model was assumed. The best fit can be expressed as 

R{T) = (4.98 ± 0.7)1(T 5 + (0.077 ± 3.0)10~ 7 T 2 (5.21) 
+(3.10 ± 0.23)10~ 7 [T 3 J 3 (0 D /T)/7.212] + (1.84 ± 0.26)10~ 10 [T 5 J 5 (0 D /T)/124.4], 

where J3 and J5 are integrals occurring in the Wilson and Bloch theories^ and the best value 
for 0£>, the effective Debye temperature, is (270 ± W)K. Over most of the temperature range 
below 300 K, the T 3 Wilson term dominates. Thus it was concluded that inter-band scattering 
is quite important in niobium. Because of the large magnitude of inter-band scattering, it was 
difficult to determine the precise amount of T 2 dependence in the resistivity. Measurements of 
the electrical resistivity of the high purity specimens of niobium were carried out in RefsP^tU^ 
It was shown that Mott theory is obeyed at high temperature in niobium. In particular, the 
resistivity curve reflects the variation of the density of states at the Fermi surface when the 
temperature is raised, thus demonstrating the predominance of s — d transitions. In addition, 
it was found impossible to fit a Bloch-Gruneisen or Wilson relation to the experimental curve. 
Several arguments were presented to indicate that even a rough approximation of the Debye 
temperature has no physical significance and that it is necessary to take the Umklapp processes into 
account. Measurements of low-temperature electrical and thermal resistivity of tungsteri^"^' and 
vanadiuuP^l showed the effects of the electron-electron scattering between different branches of 
the Fermi surface in tungsten and vanadium, thus concluding that electron-electron scattering 
does contribute measurable to electrical resistivity of these substances at low temperature. 
In transition metal compounds, e.g. MnP the electron-electron scattering is attributeoP^ to be 
dominant at low temperatures, and furthermore the 3d electrons are thought to carry electric 
current. It is remarkable that the coefficient of the T 2 resistivity is very large, about a hundred 
times those of Ni and Pd in which s electrons coexist with d electrons and electric current is mostly 
carried by the s electrons. This fact suggests strongly that in MnP s electrons do not exist at the 
Fermi level and current is carried by the 3d electrons. This is consistent with the picture^^ that 
in transition metal compounds the s electrons are shifted up by the effect of antibonding with the 
valence electrons due to a larger mixing matrix, compared with the 3d electrons, caused by their 
larger orbital extension. 

It should be noted that the temperature coefficients of resistance can be positive and negative 
in different materials. A semiconductor material exhibits the temperature dependence of the 
resistivity quite different than in metal. A qualitative explanation of this different behavior follows 
from considering the number of free charge carriers per unit volume, n, and their mobility, fi. In 
metals n is essentially constant, but /i decreases with increasing temperature, owing to increased 
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lattice vibrations which lead to a reduction in the mean free path of the charge carriers. This 
decrease in mobility also occurs in semiconductors, but the effect is usually masked by a rapid 
increase in n as more charge carriers are set free and made available for conduction. Thus, intrinsic 
semiconductors exhibit a negative temperature coefficient of resistivity. The situation is different 
in the case of extrinsic semiconductors in which ionization of impurities in the crystal lattice is 
responsible for the increase in n. In these, there may exist a range of temperatures over which 
essentially all the impurities are ionized, that is, a range over which n remains approximately 
constant. The change in resistivity is than almost entirely due to the change in fj,, leading to a 
positive temperature coefficient. 

It is believed that the electrical resistivity of a solid at high temperatures is primarily due to the 
scattering of electrons by phonons and by impurities P It is usually assumed, in accordance with 
Matthiessen rule that the effect of these two contributions to the resistance are simply additive. At 
high temperature (not lower than Debye temperature) lattice vibrations can be well represented by 
the Einstein model. In this case, I/V2 ~ T, so that I/02 ~ T. If the properties and concentration 
of the lattice defects are independent of temperature, then \/o\ is also independent of temperature 
and we obtain 

l/a~a + bT, (5.22) 

where a and b are constants. However, this additivity is true only if the effect of both impurity 
and phonon scattering can be represented by means of single relaxation times whose ratio is in- 
dependent of velocity P^ It was showrP^Sl that the addition of impurities will always decrease 
the conductivity. Investigations of the deviations from Matthiessen rule at high temperatures in 
relation to the electron-phonon interaction were carried out in Refs .USSHIMI n was gnown [159] j n 
particular, that changes in the electron-phonon interaction parameter A, due to dilute impurities 
were caused predominantly by interference between electron-phonon and electron-impurity scat- 
tering. 

The electronic band structures of transition metals are extremely complicated and make calcula- 
tions of the electrical resistivity due to structural disorder and phonon scattering very difficult. In 
addition the nature of the electron-phonon matrix elements is not well understood P^ The analy- 
sis of the matrix elements for scattering between states was performed in RefP^l n W as concluded 
that even in those metals where a fairly spherical Fermi surface exists, it is more appropriate to 
think of the electrons as tightly bound in character rather than free electron-like. In addition, the 
"single site" approximations are not likely to be appropriate for the calculation of the transport 
properties of structurally-disordered transition metals. 

5.4 Conductivity of alloys 

The theory of metallic conduction can be applied for explaining the conductivity of alloy g.E2HlZ2| 
According to the Bloch-Gruneisen theory, the contribution of the electron-phonon interaction to 
the dc electrical resistivity of a metal at high temperatures is essentially governed by two factors, 
the absolute square of the electron-phonon coupling constant, and the thermally excited mean 
square lattice displacement. Since the thermally excited mean lattice displacement is proportional 
to the number of phonons, the high temperature resistivity R is linearly proportional to the 
absolute temperature T, and the slope dR/dT reflects the magnitude of the electron-phonon 
coupling constant. However, in many high resistivity metallic alloys, the resistivity variation 
dR/dT is found to be far smaller than that of the constituent materials. In some cases dR/dT 
is not even always positive. There are two types of alloy, one of which the atoms of the different 
metals are distributed at random over the lattice points, another in which the atoms of the 
components are regularly arranged. Anomalous behavior in electrical resistivity was observed in 
many amorphous and disordered substances r 68 ^ 173 ^ At low temperatures, the resistivity increases 
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in T 2 instead of the usual T 5 dependence. Since T 2 dependence is usually observed in alloys which 
include a large fraction of transition metals, it has been considered to be due to spins. In some 
metals, T 2 dependence might be caused by spins. However, it can be caused by disorder itself. The 
calculation of transport coefficients in disordered transition metal alloys become a complicated 
task if the random fluctuations of the potential are too large. It can be shown that strong potential 
fluctuations force the electrons into localized states. Another anomalous behavior occurs in highly 
resistive metallic systems^^^l which is characterized by small temperature coefficient of the 
electrical resistivity, or by even negative temperature coefficient. 

According to Matthiessen rule 1129JH3Q] ^ e electrical resistance of a dilute alloy is separable into a 
temperature-dependent part, which is characteristic of the pure metal, and a residual part due to 
impurities. The variation with temperature of the impurity resistance was calculated by Taylor 
The total resistance is composed of two parts, one due to elastic scattering processes, the other to 
inelastic ones. At the zero of temperature the resistance is entirely due to elastic scattering, and 
is smaller by an amount 70 than the resistance that would be found if the impurity atom were 
infinitely massive. The factor 70 is typically of the order of 10~ 2 . As the temperature is raised the 
amount of inelastic scattering increases, while the amount of elastic scattering decreases. However, 
as this happens the ordinary lattice resistance, which varies as T 5 , starts to become appreciable. 
For a highly impure specimen for which the lattice resistance at room temperature, Rg, is equal 
to the residual resistance, Rq, the total resistance at low temperatures will have the form 

i?(T)«10- 2 (^) 2 + 500(^) 5 + i? - (5.23) 

The first term arises from incoherent scattering and the second from coherent scattering, accord- 
ing to the usual Bloch-Gruneisen theory. It is possible to see from this expression that T 2 term 
would be hidden by the lattice resistance except at temperatures below 6 /AO. This represents a 
resistance change of less than 10 _5 i?o, an d is not generally really observable. 
In disordered metals the Debye- Waller factor in electron scattering by phonons may be an origin 
for negative temperature coefficient of the resistivity. The residual resistivity may decreases as T 2 
with increasing temperature because of the influence of the Debye- Waller factor. But resulting 
resistivity increases as T 2 with increasing temperature at low temperatures even if the Debye- 
Waller factor is taken into account. It is worthy to note that the deviation from Matthiessen 
rule in electrical resistivity is large in the transition metal alloys^SE^ and dilute alloys .EZttiZMI T n 
certain cases the temperature dependence of the electrical resistivity of transition metal alloys at 
high temperatures can be connected with change electronic density of statesP^ The electronic 
density of states for V — Cr, Nb — Mo and Ta — W alloys have been calculated in the coherent 
potential approximation. From these calculated results, temperature dependence of the electrical 
resistivity R at high temperature have been estimated. It was shown that the concentration vari- 
ation of the temperature dependence in R/T is strongly dependent on the shape of the density of 
states near the Fermi level. 

Many amorphous metals and disordered alloys exhibit a constant or negative temperature coef- 
ficient of the electrical resistivity^^-^ in contrast to the positive temperature coefficient of the 
electrical resistivity of normal metals. Any theoretical models of this phenomenon must include 
both the scattering (or collision) caused by the topological or compositional disorder, and also the 
modifications to this collision induced by the temperature or by electron-phonon scattering. If 
one assumes that the contributions to the resistivity from scattering mechanisms other than the 
electron-phonon interaction are either independent of T, like impurity scattering, or are saturated 
at high T, like magnetic scattering, the correlation between the quenched temperature dependence 
and high resistivity leads one to ask whether the electron-phonon coupling constant is affected by 
the collisions of the electrons. 
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The effect of collisions on charge redistributions is the principal contributor to the electron-phonon 
interaction in metals. It is studied as a mechanism which could explain the observed lack of tem- 
perature dependence of the electrical resistivity of many concentrated alloys. The collision time- 
dependent free electron deformation potential can be derived from a self-consistent linearized 
Boltzmann equation. The results indicates that the collision effects are not very important for 
real systems. It can be understood assuming that the charge redistribution produces only a neg- 
ligible correction to the transverse phonon-electron interaction. In addition, although the charge 
shift is the dominant contribution to the longitudinal phonon-electron interaction, this deforma- 
tion potential is not affected by collisions until the root mean square electron diffusion distance 
in a phonon period is less than the Thomas-Fermi screening length. This longitudinal phonon- 
electron interaction reduction requires collision times of the order of 10 -19 sec in typical metals 
before it is effective. Thus, it is highly probable that it is never important in real metals. Hence, 
this collision effect does not account for the observed, quenched temperature-dependence of the 
resistivity of these alloys. However, these circumstances suggest that the validity of the adiabatic 
approximation, i.e., the Born-Oppenheimer approximation, should be relaxed far beyond the pre- 
viously suggested criteria. All these factors make the proper microscopic formulation of the theory 
of the electron-phonon interaction in strongly disordered alloys a very complicated problem. A 
consistent microscopic theory of the electron-phonon interaction in substitutionally disordered 
crystalline transition metal alloys was formulated by Wysokinski and KuzemskjP^Sl within the 
MTBA. This approach combines the Barisic, Labbe and Friedel modeP^ with the more complex 
details of the CPA (coherent potential approximation). 

The low-temperature resistivity of many disordered paramagnetic materials often shows a T 3 / 2 
rather than a T 2 dependence due to spin-fluctuation-scattering resistivity. The coefficient of the 
T 3 / 2 term often correlates with the magnitude of the residual resistivity as the amount of disorder is 
varied. A model calculations that exhibits such behavior were carried our in RefP^E In the absence 
of disorder the spin-fluctuations drag suppresses the spin-fluctuation T 2 term in the resistivity. 
Disorder produces a finite residual resistivity and also produces a finite spin-fluctuation-scattering 
rate. 



5.5 Magnetoresistance and the Hall effect 

The Hall effect and the magnetoresistanc e 1 182 ^ ^ are the manifestations of the Lorentz force on a 
subsystem of charge carries in a conductor constrained to move in a given direction and subjected 
to a transverse magnetic field. Let us consider a confined stream of a carriers, each having a 
charge e and a steady-state velocity v x due to the applied electric field E x . A magnetic field H in 
the z direction produces a force F y which has the following form 



F = e[E + (l/c)v x Hj . (5.24) 
The boundary conditions lead to the equalities 

F y = = E y - (l/c)v x H z . (5.25) 
The transverse field E y is termed the Hall field E y = E H and is given by 

E H = (l/c)v x H z = J x = nev x , (5.26) 

nec 

where J x is the current density and n is the charge carrier concentration. The Hall field can be 
related to the current density by means of the Hall coefficient Rh 

E H = R H J X H Z ; R H = — . (5.27) 

nec 
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The essence of the Hall effect is that Hall constant is inversely proportional to the charge carrier 
density n, and that is negative for electron conduction and positive for hole conduction. A useful 
notion is the so-called Hall angle which is defined by the relation 

9 = tan' 1 (E y /E x ) . (5.28) 

Thus the Hall effect may be regarded as the rotation of the electric field vector in the sample as a 
result of the applied magnetic field. The Hall effect is an effective practical tool for studying the 
electronic characteristics of solids. The above consideration helps to understand how thermomag- 
netic effects^^"^ can arise in the framework of simple free-electron model. The Lorentz force acts 
as a velocity selector. In other words, due to this force the slow electrons will be deflected less 
than the more energetic ones. This effect will lead to a temperature gradient in the transverse 
direction. This temperature difference will result in a transverse potential difference due to the 
Seebeck coefficient of the material. This phenomenon is called the Nernst-Ettingshausen effect 
It should be noted that the simple expression for the Hall coefficient Rh is the starting point only 
for the studies of the Hall effect in metals and alloys JESlLLMI j^. implies Rh is temperature inde- 
pendent and that E H varies linearly with applied field strength. Experimentally the dependence 
Rh = 1/nec do not fit well the situation in any solid metal. Thus there are necessity to explain 
these discrepancies. One way is to consider an effective carrier density n*(n) which depends on 
n, where n is now the mean density of electrons calculated from the valency. This interrelation is 
much more complicated for the alloys where n*(n) is the function of the concentration of solute 
too. It was shown that the high-field Hall effect reflects global properties of the Fermi surface such 
as its connectivity, the volume of occupied phase space, etc. The low-field Hall effect depends 
instead on microscopic details of the dominant scattering process. A quantum-mechanical theory of 
transport of charge for an electron gas in a magnetic field which takes account of the quantization 
of the electron orbits has been given by Ar gyres. 

Magnetoresistance^E^IE^E^ is an important galvanomagnetic effect which is observed in a wide 
range of substances and under a variety of experimental conditions JSSfflsSI The transverse mag- 
netoresistance is defined by 

, , R(H) - R AR . 
Qmr(H) = = — , (5.29) 

where R(H) is the electrical resistivity measured in the direction perpendicular to the magnetic 
field H, and R is the resistivity corresponding to the zero magnetic field. The zero-field resistivity 
R is the inverse of the zero-field conductivity and is given approximately by 

according to the simple kinetic theory applied to a single-carrier system. Here e, m* , n, (v) and I 
are respectively, charge, effective mass, density, average speed, and mean free path of the carrier. 
In this simplified picture the four characteristics, e, m* , n, and (v), are unlikely to change sub- 
stantially when a weak magnetic field is applied. The change in the mean free path I should then 
approximately determine the behavior of the magnetoresistance AR/R at low fields. 
The magnetoresistance practically of all conducting pure single-crystals has been experimentally 
found to be positive and a strong argument for this were given on the basis of nonequilibrium sta- 
tistical mechanics.^ In some substances, e.g. carbon, CdSe, E^CuSis, etc., m agnetoresistance 
is negative while in CdMnSe is positive and much stronger than in CdSe^^221 A qualitative in- 
terpretation of the magnetoresistance suggests that those physical processes which make the mean 
free path larger for greater values of H should contribute to the negative magnetoresistance. Mag- 
netic scattering leads to negative magnetoresistance^^ characteristic for ferro- or paramagnetic 
case, which comes from the suppression of fluctuation of the localized spins by the magnetic field. 
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Figure 3: Schematic form of the thermal conductivity of various materials 



A comprehensive derivation of the quantum transport equation for electric and magnetic fields was 
carried out by MahanPSSl More detailed discussion of the various aspects of theoretical calculation 
of the magnetoresistance in concrete substances are given in RefsJ 188 ^ 198 ^ 200 " 202 ^ 



5.6 Thermal conduction in solids 

Electric and thermal conductivities are intimately connected since the thermal energy also is 
mainly transported by the conduction electrons. The thermal conductivity^^] G f a variety of 
substances, metals and nonmetals, depends on temperature region and varies with temperature 
substantially^^ (see Fig. 3). Despite a rough similarity in the form of the curves for metallic 
and nonmetallic materials there is a fundamental difference in the mechanism whereby heat is 
transported in these two types of materials. In metals^Ql] heat is conducted by electrons; in 
nonmetals3"2SSl it is conducted through coupled vibrations of the atoms. The empirical datcP^ 
show that the better the electrical conduction of a metal, the better its thermal conduction. Let 
us consider a sample with a temperature gradient dT/dx along the x direction. Suppose that the 
electron located at each point x has thermal energy E(T) corresponding to the temperature T at 
the point x. It is possible to estimate the net thermal energy carried by each electron as 

E(T) - E(T + ^-tv cos 9) = -^p^rt> cos 9. (5.31) 
v 1 v dx ' dT dx y 1 

Here we denote by 6 the angle between the propagation direction of an electron and the x direction 
and by v the average speed of the electron. Then the average distance travelled in the x direction 
by an electron until it scatters is tvcosO. The thermal current density J q can be estimated as 

dE(T) dT 9 , , „, ,„ . 

J, = - n ^rte TV (cos e) > (5 - 32) 

where n is the number of electrons per unit volume. If the propagation direction of the electron 
is random, then (cos 2 6) = 1/3 and the thermal current density is given by 

1 dE(T) 2 dT 
J a = — n — — — tv — . (5.33) 
9 3 dT dx ( ' 
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Here ndE(T)/dT is the electronic heat capacity C e per unit volume. We obtain for the thermal 
conductivity n the following expression 

K = -C e rv 2 ; Q = -K-. (5.34) 
3 ax 

The estimation of k for a degenerate Fermi distribution can be given by 

where l/2mv 2 = 3/5Co and ks is the Boltzmann constant. It is possible to eliminate nr with the 
aid of equality r = ma/ne 2 . Thus we obtain^ 

- = —-%-T. 5.36 
a 5 e z 

This relation is called by the Wiedemann-Franz law. The more precise calculation gives a more 
accurate factor value for the quantity 7r 2 /5 as 7r 2 /3. The most essential conclusion to be drawn 
from the Wiedemann-Franz law is that k/u is proportional to T and the proportionality constant 
is independent of the type of metal. In other words, a metal having high electrical conductivity has 
a high thermal conductivity at a given temperature. The coefficient k/ctT is called the Lorentz 
number. At sufficiently high temperatures, where a is proportional to 1/T, k is independent 
of temperature. Qualitatively, the Wiedemann-Franz law is based upon the fact that the heat 
and electrical transport both involve the free electrons in the metal. The thermal conductivity 
increases with the average particle velocity since that increases the forward transport of energy. 
However, the electrical conductivity decreases with particle velocity increases because the collisions 
divert the electrons from forward transport of charge. This means that the ratio of thermal to 
electrical conductivity depends upon the average velocity squared, which is proportional to the 
kinetic temperature. 

Thus, there are relationships between the transport coefficients of a metal in a strong magnetic 
field and a very low temperatures. Examples of such relations are the Wiedemann-Franz law for 
the heat conductivity k, which we rewrite in a more general form 

k = LTa, (5.37) 

and the Mott rule^EJ for the thermopower S 

S = eLT(r- 1 ^. (5.38) 
d/j, 

Here T is the temperature, \i denotes the chemical potential. The Lorentz number L = 1/3(7t/cb /e) 2 , 
where kg is the Boltzmann constant, is universal for all metals. Useful analysis of the Wiedemann- 
Franz law and the Mott rule was carried out by NagaevP^J 



6 Linear Macroscopic Transport Equations 

We give here a brief refresher of the standard formulation of the macroscopic transport equations 
from the most general point of viewPSl One of the main problems of electron transport the- 
ory is the finding of the perturbed electron distribution which determines the magnitudes of the 
macroscopic current densities. Under the standard conditions it is reasonably to assume that the 
gradients of the electrochemical potential and the temperature are both very small. The macro- 
scopic current densities are then linearly related to those gradients and the ultimate objective of 
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Table 4: Fluxes and Generalized Forces 

Table 4. Fluxes and Generalized Forces 



Process 


Flux 


Generalized force 


Tensor character 


Electrical conduction 


Je 


V<p 


vector 


Heat conduction 


■h 


V(l/T) 


vector 


Diffusion 


Diffusion flux J p 


-(l/T)[Vn] 


vector 


Viscous flow 


Pressure tensor P 


-(1/T)W 


Second-rank tensor 


Chemical reaction 


Reaction rate W r 


Affinity a r /T 


Scalar 



the theory of transport processes in solids (see Table 4). Let 77 and T denote respectively the elec- 
trochemical potential and temperature of the electrons. We suppose that both the quantities vary 
from point to point with small gradients V77 and VT. Then, at each point in the crystal, electric 
and heat current densities J e and J q will exist which are linearly related to the electromotive force 
£ = 1/eVr/ and VT by the basic transport equations 

J e = L u £ + L l2 VT, (6.1) 

J q = L 2l £ + L 22 VT. (6.2) 

The coefficients Ln, L12, T21 and L22 in these equations are the transport coefficients which 
describe irreversible processes in linear approximation. We note that in a homogeneous isothermal 
crystal, £ is equal to the applied electric field E. The basic transport equations in the form (16.11) 
and (|6.2p describe responses J e and J q under the influence of £ and VT. The coefficient L\\ = a 
is the electrical conductivity. The other three coefficients, L12, T21 and L22 have no generally 
accepted nomenclature because these quantities are hardly ever measured directly. From the 
experimental point of view it is usually more convenient to fix J e and VT and then measure £ 
and J q . To fit the experimental situation the equations (16. ip and (I6.2p must be rewritten in the 
form 

£ = RJ e + SVT, (6.3) 
J q = UJ e - kVT, (6.4) 

where 

R = a~ 1 , (6.5) 

S = -a~ 1 L 12 , (6.6) 

n = L 2l a-\ (6.7) 

K = L 2 i<T _1 Li2 = L22, (6.8) 

which are known respectively as the resistivity, thermoelectric power, Peltier coefficient and ther- 
mal conductivity. These are the quantities which are measured directly in experiments. 
All the coefficients in the above equations are tensors of rank 2 and they depend on the magnetic 
induction field B applied to the crystal. By considering crystals with full cubic symmetry, when 
B = 0, one reduces to a minimum the geometrical complications associated with the tensor char- 
acter of the coefficients. In this case all the transport coefficients must be invariant under all the 
operations in the point group m3mP^ This high degree of symmetry implies that the coefficients 
must reduce to scalar multiples of the unit tensor and must therefore be replaced by scalars. When 
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B / the general form of transport tensors is complicated even in cubic crystals JSOHS] T n the 
case, when an expansion to second order in B is sufficient the conductivity tensor takes the form 

a a p(B) - a a/3 (0) + ^ 1)B~ U + 2 ^ dB u B u B ^ v + " " " (6 ' 9) 

Here the (a/3)— th element of a referred to the cubic axes (Oxyz). For the case when it is possible 
to confine ourselves by the proper rotations in m3m only, we obtain that 

^(O) = a 6 a/ 3, (6.10) 

where o"o is the scalar conductivity when B = 0. We also find that 

dB^ ~ ?£a/37 ' 

OB^B u 

where C are all scalar and e a ^~ is the three-dimensional alternating symbol.'^' Thus we 

obtain a relation between j and £ (with VT = 0) 

j = cr Q £ |^xfi| £i? 2 £ + rjB(B£) + C$£, (6.11) 



where $ is a diagonal tensor with <3> QQ , = i? 2 (see RefsJ^E^J). The most interesting transport 
phenomena is the electrical conductivity under homogeneous isothermal conditions. In general, 
the calculation of the scalar transport coefficients &o,<;,£ t ,r],( i is complicated task. As was men- 
tioned above, these coefficients are not usually measured directly. In practice, one measures the 
corresponding terms in the expression for £ in terms of j up to terms of second order in B. To 
show this clearly, let us iterate the equation (|6.1ip . We then find that 

£ = R j - R H J xB + Rq {bB 2 j + cB(Bj) + d$j) , (6.12) 

where 

R = ctq" 1 , R h = Gq \, (6.13) 
are respectively the low-field resistivity and Hall const an t(inillS2l[IS3lllSZI anc j 

b = -R (( + R o ey, c=-Ro(r]-R o ey, d = -g c (6.14) 

are the magnetoresistance coefficients.^ These are the quantities which are directly measured. 



7 Statistical Mechanics and Transport Coefficients 

The central problem of nonequilibrium statistical mechanics is to derive a set of equations which 
describe irreversible processes from the reversible equations of motionpZJSSlEZlErJHSS] rj^g consis- 
tent calculation of transport coefficients is of particular interest because one can get information 
on the microscopic structure of the condensed matter. There exist a lot of theoretical methods 
for calculation of transport coefficients as a rule having a fairly restricted range of validity and 
applicability 032130111 The most extensively developed theory of transport processes is that based 
on the Boltzmann equationP ^^^ ^^^ However, this approach has strong restrictions and 
can reasonably be applied to a strongly rarefied gas of point particles. For systems in the state of 
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statistical equilibrium, there is the Gibbs distribution by means of which it is possible to calculate 
an average value of any dynamical quantity. No such universal distribution has been formulated 
for irreversible processes. Thus, to proceed to the solution of problems of statistical mechanics of 
nonequilibrium systems, it is necessary to resort to various approximate methods. During the last 
decades, a number of schemes have been concerned with a more general and consistent approach 
to transport theory H21HSIH5j[2iT}t2i5Ji2i7| ^his field i s very active and there are many aspects to the 
problemP^^Z^SSl 

7.1 Variational principles for transport coefficients 

The variational principles for transport coefficients are the special techniques for bounding trans- 
port coefficients, originally developed by Kohler, Sondheimer, Ziman, etc. (see RefsP^™). This 
approach is equally applicable for both the electronic and thermal transport. It starts from a 
Boltzmann-like transport equation for the space-and-time-dependent distribution function fg or 
the occupation number n q {r,t) of a single quasiparticle state specified by indices q (e.g., wave 
vector for electrons or wave vector and polarization, for phonons). Then it is necessary to find or 
fit a functional F[fff, n q (f, t)] which has a stationary point at the distribution fx, n q (r, t) satisfying 
the transport equation, and whose stationary value is the suitable transport coefficient. By evalu- 
ating F for a distribution only approximately satisfying the transport equation, one then obtains 
an upper or lower bound on the transport coefficients. Let us mention briefly the phonon-limited 
electrical resistivity in metalsP^ With the neglect of phonon drag, the electrical resistivity can 
be written 

1 $ ${^ k -®g?W{k,ki)dkdk> 
~2k B T \J e ^( d fO/de(k))dk\2 

Here W(k, k') is the transition probability from an electron state k to a state k! , is the electron 
velocity and /° is the equilibrium Fermi-Dirac statistical factor. The variational principle 2 tell us 
that the smallest possible value of the right hand side obtained for any function is also the 
actual resistivity. In general we do not know the form of the function that will give the right 
hand side its minimal value. For an isotropic system the correct choice is = uk, where u is 
a unit vector along the direction of the applied field. Because of its simple form, this function is 
in general used also in calculations for real systems. The resistivity will be overestimated but it 
can still be a reasonable approximation. This line of reasoning leads to the Ziman formula for the 
electrical resistivity 

R< 3vr v r rd 2 kd 2 k' (k-k') 2 u u (g)Al(k,k') 

~ e 2 k B TS 2 k 2 F ^ J J v v> (exp(hu u (q)/k B T-l)(l-eM-^u{q)/k B T)' l ' J 

Here kp is the average magnitude of the Fermi wave vector and S is the free area of the Fermi 
surface. It was shown in Ref.^3 that the formula for the electrical resistivity should not contains 
any electron-phonon enhancements in the electron density of states. The electron velocities in 
Eq. (|7.2p are therefore the same as those in Eq. (]7.ip . 

7.2 Transport theory and electrical conductivity 

Let us summarize the results of the preceding sections. It was shown above that in zero magnetic 
field, the quantities of main interest are the conductivity a (or the electrical resistivity, R = 
1/cr) and the thermopower S. When a magnetic field B is applied, the quantity of interest is a 
magnetor esist ance 

R(B) - R(B = 0) 
0MR = R(B^0) (7 - 3) 
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and the Hall coefficient Rh- The third of the (generally) independent transport coefficients is the 
thermal conductivity k. The important relation which relates k to R at low and high temperatures 
is the Wiedemann-Franz lawP^Sl T n s i m pl e metals and similar metallic systems, which have well- 
defined Fermi surface it is possible to interpret all the transport coefficients mentioned above, 
the conductivity (or resistivity), thermopower, magnetoresistance, and Hall coefficient in terms 
of the rate of scattering of conduction electrons from initial to final states on the Fermi surface. 
Useful tool to describe this in an approximate way is the Boltzmann transport equation, which, 
moreover usually simplified further by introducing a concept of a relaxation time. In the approach 
of this kind we are interested in low-rank velocity moments of the distribution function such as 
the current 

vf(v)d 3 v. (7.4) 



In the limit of weak fields one expects to find Ohm law j = oE. The validity of such a formulation of 
the Ohm law was analyzed by Bakshi and Gross) 2 ^ This approach was generalized and developed 
by many authors. The most popular kind of consideration starts from the linearized Boltzmann 
equation which can be derived assuming weak scattering processes. For example, for the scattering 
of electrons by "defect" (substituted atom or vacancy) with the scattering potential V d (r) the 
perturbation theory gives 

3vr 2 m 2 lW r 2kF tll „ , .^^mw. 



R~ 4n3e2fc 2 J Q \(k + q\V d (m)\ 2 q"dq, (7.5) 

where N and Oq are the number and volume of unit cells, kp the magnitude of the Fermi wavevec- 
tor, the integrand (A; + q\V d (r)\k) represents the matrix elements of the total scattering poten- 
tial V d (f), and the integration is over the magnitude of the scattering wavevector q defined by 
q = k — k' . 

Lax have analyzed in detail the general theory of carriers mobility in solidsP^ Luttinger and 
Kohnp^ESl GreenwoodpSS Chester and Thellungj^SH and FujiitsPZI developed approaches to the 
calculation of the electrical conductivity on the basis of the generalized quantum kinetic equations. 
The basic theory of transport for the case of scattering by static impurities has been given in the 
works of Kohn and Luttinger^^^ and GreenwoocP^ (see also Ref! 2 ^) . In these works the usual 
Boltzmann transport equation and its generalizations were used to write down the equations for 
the occupation probability in the case of a weak, uniform and static electric field. It was shown 
that in the case of static impurities the exclusion principle for the electrons has no effect at all 
on the scattering term of the transport equation. In the case of scattering by phonons, where the 
electrons scatter inelastically, the exclusion principle plays a very important role and the transport 
problem is more involved. On the other hand transport coefficients can be calculated by means 
of theory of the linear response such as the Kubo formulae for the electrical conductivity. New 
consideration of the transport processes in solids which involve weak assumptions and easily gen- 
eralizable methods are of interest because they increase our understanding of the validity of the 
equations and approximations used-^^E!' 2 ^ Moreover, it permits one to consider more general 
situations and apply the equations derived to a variety of physical systems. 



8 The Method of Time Correlation Functions 

The method of time correlation functiongSZlBSlBSl^ 1 221 11 223 112 ^ 12 ^ 1 is an attempt to base a linear 
macroscopic transport equation theory directly on the Liouville equation. In this approach one 
starts with complete iV-particle distribution function which contains all the information about the 
system. In the method of time correlation functions it is assumed that the iV-particle distribution 
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function can be written as a local equilibrium iV-particle distribution function plus correction 
terms. The local equilibrium function depends upon the local macroscopic variables, temperature, 
density and mean velocity and upon the position and momenta of the N particles in the system. 
The corrections to this distribution functions is determined on the basis of the Liouville equation. 
The main assumption is that at some initial time the system was in local equilibrium (quasi- 
equilibrium) but at later time is tending towards complete equilibrium. It was shown by many 
authors (for comprehensive review see RefsP^ESKD) that the suitable solutions to the Liouville 
equation can be constructed and an expression for the corrections to local equilibrium in powers 
of the gradients of the local variables can be found as well. The generalized linear macroscopic 
transport equations can be derived by retaining the first term in the gradient expansion only. In 
principle, the expressions obtained in this way should depend upon the dynamics of all N particles 
in the system and apply to any system, regardless of its density. 



8.1 Linear response theory 

The linear response theory was anticipated in many works (see Ref s jm[ Pll^ |2371|238| £ Qr details) 
on the theory of transport phenomena and nonequilibrium statistical mechanics. The important 
contributions have been made by many authors. By solving the Liouville equation to the first order 
in the external electric field, Kubo 239 " 2 ' formulated an expression for the electric conductivity in 
microscopic terms. 

He used linear response theory to give exact expressions for transport coefficients in terms of 
correlation functions for the equilibrium system. To evaluate such correlation functions for any 
particular system, approximations have to be made. 

In this section we shall formulate briefly some general expressions for the conductivity tensor 
within the linear response theory. Consider a many-particle system with the Hamiltonian of a 
system denoted by H. This includes everything in the absence of the field; the interaction of the 
system with the applied electric field is denoted by H ex t. The total Hamiltonian is 

U = H + H ext . (8.1) 

The conductivity tensor for an oscillating electric field will be expressed in the forrrP321 

<J^ U = / Tr Po j u (0)j^t + ihX)e- iujt dtd\, (8.2) 
Jo Jo 

where po is the density matrix representing the equilibrium distribution of the system in absence 
of the electric field 

p = e^ H /[Tre^ H ], (8.3) 

/3 being equal to 1/ksT. Here j^, j v are the current operators of the whole system in the p, v direc- 
tions respectively, and jn represents the evolution of the current as determined by the Hamiltonian 
H 

Ut) = e^V"^. (8.4) 

Kubo derived his expression Eq. (|8.2p by a simple perturbation calculation. He assumed that at 
t = — oo the system was in the equilibrium represented by po- A sinusoidal electric field was 
switched on at t = — oo, which however was assumed to be sufficiently weak. Then he considered 
the equation of motion of the form 

ih^- t p=[H + H ext (t),p]. (8.5) 
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The change of p to the first order of H ex t is given by 

P- P0 = ^ f e^ m 'l ih )[H ext {t%P^ m ' m + 0(H ext ). (8.6) 
in 7_oo 

Therefore the averaged current will be written as 

(j^)) = \[ t Tr[H ext (t'),p )U-t')dt', (8.7) 

where H ex t(t') will be replaced by —edE(t'), being the total dipole moment of the system. 
Using the relation 

[A,e-f> H ] = -e-f> H ^ e XH [A,p]e- XH d\, (8.8) 
i Jo 

the expression for the current can be transformed into Eq. (|8.7|) . The conductivity can be also 
written in terms of the correlation function (j v (P)ju(t))o. The average sign (. . .)q means the average 
over the density matrix p . 

The correlation of the spontaneous currents may be described by the correlation functionP^ 

S^(t) = 0,(0)^(t)) = (jAr)j^(t + r)) . (8.9) 

The conductivity can be also written in terms of these correlation functions. For the symmetric 
(" s" ) part of the conductivity tensor KubcP^l derived a relation of the form 

1 f°° 

Rea s (uj) = — -— Bavit) cos cotdt, (8.10) 

£ ^{ UJ ) Jo 

where is the average energy of an oscillator with the frequency oj at the temperature 

T = l/kB/3. This equation represents the so-called fluctuation-dissipation theorem, a particular 
case of which is the Nyquist theorem for the thermal noise in a resistive circuit. The fluctuation- 
dissipation theorems were established^"^] f or systems in thermal equilibrium. It relates the 
conventionally defined noise power spectrum of the dynamical variables of a system to the corre- 
sponding admittances which describe the linear response of the system to external perturbations. 
The linear response theory is very general and effective tool for the calculation of transport co- 
efficients of the systems which are rather close to a thermal equilibrium. Therefore, the two 
approaches, the linear response theory and the traditional kinetic equation theory share a domain 
in which they give identical results. A general formulation of the linear response theory was given 
by Kubol™I3]for the case of mechanical disturbances of the system with an external source in 
terms of an additional Hamiltonian. 

A mechanical disturbance is represented by a force F(t) acting on the system which may be given 
function of time. The interaction energy of the system may then be written as 

H ext (t) = —AF(t), (8.11) 

where A is the quantity conjugate to the force F. The deviation of the system from equilibrium is 
observed through measurements of certain physical quantities. If AB(t) is the observed deviation 
of a physical quantity B at the time t, we may assume, if only the force F is weak enough, a linear 
relationship between AB(t) and the force F(t), namely 

AB(t)= f <j> BA {t,t')F(t')dt', (8.12) 

J —oo 
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where the assumption that the system was in equilibrium at t = — oo, when the force had been 
switched on, was introduced. This assumption was formulated mathematically by the asymptotic 
condition, 

F(t) ~ e £t as t -> -oo (e > 0). (8.13) 

Eq. (|8.12p assumes the causality and linearity. Within this limitation it is quite general. Kubo 
called the function efts A of response function of B to F, because it represents the effect of a delta- 
type disturbance of F at the time t' shown up in the quantity B at a later time t. Moreover, as it 
was claimed by Kubo, the linear relationship (|8.12p itself not in fact restricted by the assumption 
of small deviations from equilibrium. In principle, it should be true even if the system is far from 
equilibrium as far as only differentials of the forces and responses are considered. For instance, 
a system may be driven by some time-dependent force and superposed on it a small disturbance 
may be exerted; the response function then will depend both on t and t' separately. If, however, 
we confine ourselves only to small deviations from equilibrium, the system is basically stationary 
and so the response functions depend only on the difference of the time of pulse and measurement, 
t and t', namely 

cl>BA(t,t') = <l> BA (t-t'). (8.14) 
In particular, when the force is periodic in time 

F{t) = ReFe iult , (8.15) 

the response of B will have the form 

AB(t) = Re X BA^)Fe iut , (8.16) 

where Xba(w) is the admittance 

Xba{u)= [ <i)BA{t)e- iujt dt. (8.17) 

J — oo 

More precisely the response AB(t) to an external periodic force F(t) = F cos(wi) conjugate 
to a physical quantity A is given by Eq. (|8.16|) . where the admittance Xba(w) is defined as 

/"OO 

Xba(u>) = lim / cjiBA^e-^+^dt. (8.18) 

The response function 4>ba{^) is expressed as 

<t> BA (t) = iTr[A, P ]B(t) = -iTrp[A, B(t)] (8.19) 

rP rP 
= / Tr pA(-i\)B(t)dt = - TrpA(-i\)Bd\, 
Jo Jo 

where p is the canonical density matrix 

p = exp(-/3(H - «)), exp(-/3ft) = Tr exp( -/?#). (8.20) 
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In certain problems it is convenient to use the relaxation function denned by 

$BA(t)= lim / (t>BA{t')e- £t ' dt' (8.21) 
£ ^+° Jt 

/oo 
([B(t>),A]}dt> 

poo pP 

= - dt' d\(A(i\)B(t')) 
Jt Jo 

= f d\((A(-i\)B(t)) - lim (A(-iX)B(t))) 
Jo V > 



P 

d\(A(-i\)B(t)} - p lim (AB(t)) 

o 



t— >oo 



= / d\(A(-i\)B(t)) - (3(A°B°). 
Jo 

It is of use to represent the last term in terms of the matrix elements 

/ d\(A(—i\)B(t)} — (3(A°B ) (8.22) 
Jo 



i/j2 e M-P E d) Y,( n \ A \ m )( m \ B \ n ) e ~ it(En ~ Em) ■ 



E — E 

n,m ^ m Cjn 



Here \m) denotes an eigenstate of the Hamiltonian with an eigenvalue E m and A and B° are the 
diagonal parts of A and B with respect to H. The response function Xba{w) can be rewritten in 
terms of the relaxation function. We have 

Xba(w) = — lim l°° ^ BA {t)e~^ + ^dt (8.23) 

e^+OjQ 

poo 

= <j)BA(0) - ioj lim foA^e-^+^dt 

e^+0J 
poo pP J 

= - lim / dt d\e-^ +£)t ^-(AB(t + i\)) 

e^+0j J Q dt 11 

= i lim f°° dte- {iu+£ ^ ( (AB(t + 0)) - (AB{t))) 
£ ^+°Jo ^ > 

poo 

= i lim / dte-^ +e ^{[B(t),A]) 
_ A mn B nm i c -PE n _ e -PE m \ 

!,) -A- M__ -I- ?'P 



where 



w + co mn + ie 

n,m 



(m\A\n) 
(E„e-^)V2 
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In particular, the static response xba(0) is given by 

Xba(0) = 0ba(O) (8.24) 

= / dx((AB(iX))- lim (AB(t + ixW 
Jo v > 

= i lim / e- £t dt((AB(t + i/3)) - (AB(t)) 

/•oo 

= i lim / e - et dt([sm,i4]>. 
This expression can be compared with the isothermal response defined by 

X T ba = J ((AB(iX)) - (A)(B))dX. (8.25) 
The difference of the two response functions is given by 

Xba-Xba(0)= lim / dX(AB(t + iX))-(3(A)(B) (8.26) 



= /3 l Km(AB(t))-{A)(B) . 

Vt— i>oo / 

The last expression suggests that it is possible to think that the two response functions are 
equivalent for the systems which satisfy the condition 

lim(AB(i)) = (A)(B). (8.27) 

t— >oo 

It is possible to speak about these systems in terms of ergodic (or quasi-ergodic) behavior, however, 
with a certain reservation (for a recent analysis of the ergodic behavior of many-body systems see 

It may be of use to remind a few useful properties of the relaxation function. 
If A and B are both Hermitian, then 



28) 



/>oo 

&BA(t) = real, / <f> AA (t) > 0. 
J o 

The matrix-element representation of the relaxation function have the form 

$BA(t) = y^( R ™n cos(w mn ) - R 2 mn sin(uw)) (8.29) 

m,n 

+i yZ\ R mn sin(w mn ) - R 2 mn cos(a; m) 



where 



R mn — l ^{A nm B mn -\- A mn B nm )R mn , 

R mn — ~z~\A nm B mn A mn B nm )R mn , 
1% 

(yn\A\n) 



(En e-^)V2 : 
c — e~ 

R mn — » Ldmn — B m E n . 
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This matrix-element representation is very useful and informative. It can be shown that the 
relaxation function has the property 

W flA (i) = 0, (8.30) 

which follows from the odd symmetry of the matrix-element representation. The time integral of 
the relaxation function is given by 

/'OO 77"/? / \ 

/ $BA(t)dt = -Y^[A nrn B mn + B nm A mn )e-P En 5(uj mn ). (8.31) 
In particular, for A = B, we obtain 

/•oo 

/ $AA(t)dt = TT/3Y,\A nm \ 2 e-^5(u; mn ) > 0. (8.32) 

J° m,n 

It can be shown alscP^SI that if A and B are both bounded, then we obtain 

/■OO /"OO /'OO 

/ & 6A (t)dt = i / dt ([B(t'),A])dt' (8.33) 
Jo '" Jo Jt 



dt([B(t),A])=i([A,B}). 

Jo 

For A = B, we have 

/•oo 

<f> M (t)dt = 0, / dt I dX(A(-i\),A(t)) =0 (8.34) 



/•oo 

L 



10 

and 

lim (AB(t)) = 0, (8.35) 

t— >oo 

if A and .B are both bounded. 

Application of this analysis may not be limited to admittance functions! 2 ^ p or example, if one 
write a frequency dependent mobility function n(uS) as 

ft(u) = [iLU + ^uj)}- 1 , (8.36) 

the frequency-dependent friction 7(0;) is also related to a function <j)(t), which is in fact the corre- 
lation function of a random forceJ^^ An advanced analysis and generalization of the Kubo linear 
response theory was carried out in series of papers by Van Vliet and co-authorspSl Fluctuations 
and response in nonequilibrium steady state were considered within the nonlinear Langevin equa- 
tion approach by Ohta and OhkumaP^ It was shown that the steady probability current plays 
an important role for the response and time-correlation relation and violation of the time reversal 
symmetry. 



8.2 Green functions in the theory of irreversible processes 

Green functions are not only applied to the case of statistical equilibrium HZE Pol l 22 1 II259H261J They 
are a convenient means of studying processes where the deviation from the state of statistical equi- 
librium is small. The use of the Green functions permits one to evaluate the transport coefficients 
of these processes . Moreover, the transport coefficients are written in terms of Green functions 
evaluated for the unperturbed equilibrium state without explicitly having recourse to setting up 
a transport equation. The linear response theory can be reformulated in terms of double-time 
temperature-dependent (retarded and advanced) Green functions! 259 " 261 ^ We shall give a brief 
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account of this reformulatiorpSESi] an( j n s simplest applications to the theory of irreversible pro- 
cesses. 

The retarded two-time thermal Green functions arise naturally within the linear response for- 
malism, as it was shown by ZubarevP^E^ To show this we consider the reaction of a quantum- 
mechanical system with a time-independent Hamiltonian H when an external perturbation 

H ext (t) = —AF(t), (8.37) 

is switched on. The total Hamiltonian is equal to 

H = H + H ext , (8.38) 

where we assume that there is no external perturbation at limt — > — oo 

H ex t{t)\u m t^-oo = 0. (8.39) 

The last condition means that 

lim p(t) = po = e~ pH /[Fce- pH ], (8.40) 

t— >— oo 

where p(t) is a statistical operator which satisfies the equation of motion 

ih^- t p(t) = [H + H ext (t),p}, (8.41) 

This equation of motion together with the initial condition (|8,40p suggests to look for a solution 
of Ea. (|gH]) of the form 

p{t) = p + Ap(t). (8.42) 
Let us rewrite Eq. (|8.42p . taking into account that [H,p] = 0, in the following form 

d d 
ih-(p + Ap(t)) = ih-Ap(t) = (8.43) 

[H + H ext (t),p + Ap(t)] = [H,Ap(t)] + [H ext (t),p] + [H ext {t),Ap{t)}. 

Neglecting terms H ex t(t)Ap, since we have assumed that the system is only little removed from a 
state of statistical equilibrium, we get then 

ih^Apit) = [H, Ap(t)] + [H ext (t),p], (8.44) 

where 

Ap(t)\ limt ^oo = 0. (8.45) 

Processes for which we can restrict ourselves in Eq. (|8.44p to terms linear in the perturbation 
are called linear dissipative processes. For a discussion of higher-order terms it is convenient to 
introduce a transformation 

Ap{t) = e~ iHt/h Q(t)e iHt/h , (8.46) 

Then we have 

ih^Apit) = [H, Ap(t)} + e-* Ht / h (ih^Qit)) e im / h . (8.47) 
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This equation can be transformed to the following form 

ih^g(t) = [e iHt / h H ext (t)e- iHt /\p] + [e iHt / R H ext (t) e - iHt /\ <?(*)], (8.48) 

where 

e(t)|limt->-oo = 0. (8.49) 

In the equivalent integral form the the above equation reads 

rt 



Q(t) = \f d\[e tHX / h H ext (\)e- iHX/ \p] (8.50) 
+ Jh f dX l eiHX/RH ext(\)e- iHX /\ e (\)]. 
This integral form is convenient for the iteration procedure which can be written as 

e (t) = ±-f d\[e mx ' h H ext We- lHX/ \p] (8.51) 
+ (^) 2 f dX t d\'\e iHX / h H ext (\)e- iHX /\[e iHX '/ h H e ^^ 

\ltl/ j_oo J —oo L 



+ ••• 



In the theory of the linear reaction of the system on the external perturbation usually the only 
first term is retained 

M*) = ir dre- iH ^l n \H ext {r\pY H ^l n . (8.52) 

in y_ 00 

The average value of observable A is 

(A) t = Tr(Ap(t)) = Tr(Ap ) + Tr(AAp(t)) = (A) + A(A) t . (8.53) 

From this we find 

Tr [e iH ^/ h Ae- iH ^/ n H ext {T)p - H ext {r)e iH ^l R Ae~ iH ^/ n p) +... 

= 1 f dr([A(t - T),H ext ( T )}-) + ... (8.54) 
Introducing under the integral the sign function 

«*-r)-\\ * T< J: (^) 
10 if r > t, 

and extending the limit of integration to — oo < r < +oo we finally find 

A(A) t = J°° dr^e(r)([A(r), H ext (t - r)]_) + . . . (8.56) 

Let us consider an adiabatic switching on a periodic perturbation of the form 

H ext (t) = Bexp —(E + ie)t. 
in 
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The presence in the exponential function of the infinitesimal factor e > 0, e — > make for the 
adiabatic switching of the perturbation. Then we obtain 

1 f 00 -1 1 

A(A) t = exp— (E + ie)t / dr exp — (E + ie)r — 6(t){\A{t), B]\. (8.57) 
in J -oo ^n. 

It is clear that the last expression can be rewritten as 

1 f 00 —1 

A(A) t = exp—( y E + ie)t dr exp — (£ + ie)r G ret (A, B; r) = (8.58) 

i/l J —oo 

exp4(£ + ^)tG ret (A5;£) = expl(£ + * £ )i((A| J B)) i5+ie . 

Here E = fkj and ((yl|-B)) s+«e is the Fourier component of the retarded Green function ((A(t); B(r)}}. 
The change in the average value of an operator when a periodic perturbation is switched on adia- 
batically can thus be expressed in terms of the Fourier components of the retarded Green functions 
which connect the perturbation operator and the observed quantity. 
In the case of an instantaneous switching on of the interaction 

fO if t < t , 

H ext (t) = \ fo./.fcV -f (8 - 59) 

I Ln ex P( Ut/ihjVn if t>t Q , 



where Vq is an operator which does not explicitly depend on the time, we get 

A(A) t = J2 dT((A(t);Vn(T)))exp-(n + te)r, (8.60) 

i.e., the reaction of the system can also be expressed in terms of the retarded Green functions. 
Now we can define the generalized susceptibility of a system on a perturbation H ex t(t) as 

X (A, B; E) = X (A, zB; E) = lim -A(A) t exp ^-{E + ie)t = ((A\B)) E+i£ . (8.61) 

z->o z in 

In the time representation the above expression reads 

1 f°° -1 

X(A,B;E) = -J dtex V —(E + ie)t6(t)([A(t),B}_). (8.62) 

This expression is an alternative form of the fluctuation-dissipation theorem, which show explic- 
itly the connection of the relaxation processes in the system with the dispersion of the physical 
quantities. 

The particular case where the external perturbation is periodic in time and contains only one 
harmonic frequency to is of interest. Putting in that case Vt = ±hw in Eq. (l8.60p . since 

H ext {t) = -ho cos ute £t B, (8.63) 

where ho, the amplitude of the periodic force, is a c— number and where B is the operator part of 
the perturbation, we get 

A(A) t = -ho exp (J^ojt + etj {{A\B)) E=hw (8.64) 
-/i exp (-^-ut + etj ((A\B))e=-Hu ■ 
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Taking into account that (A)t is a real quantity we can write it as follows 

A(A) t = Re (x{E)h e^ Et+£t ) . (8.65) 

Here x(E) is the complex admittance, equal to 

X(E) = -27r((A\B)) E=hw . (8.66) 

These equations elucidate the physical meaning of the Fourier components of the Green func- 
tion {(A\B)) E=hu> as being the complex admittance that describes the influence of the periodic 
perturbation on the average value of the quantity A. 



8.3 The electrical conductivity tensor 

When a uniform electric field of strength £ is switched on then the perturbation acting upon the 
system of charged particles assumes the form H ext = —£ ■ d(t), where d(t) is the total dipole 
moment of the system. In this case the average operator A(t) is the current density operator j 
and the function \ is t ne complex electrical conductivity tensor denoted by cr a p(uj). If the volume 
of the system is taken to be equal to unity, then we have 

j t d a {t)= j a {t). (8.67) 

The Kubo formula (|8.65p relates the linear response of a system to its equilibrium correlation 
functions. Here we consider the connection between the electrical conductivity tensor and Green 
functions E2126QI L e ^ us start with a simplified treatment when there be switched on adiabatically 
an electrical field £{t), uniform in space and changing periodically in time with a frequency u 

£(t)=£cosut. (8.68) 

The corresponding perturbation operator is equal to 

H ext (t) = -e ^2(£rj) cos ute et . (8.69) 
j 

Here e is the charge of an electron, and the summation is over all particle coordinates fj. Under 
the influence of the perturbation there arises in the system an electrical current 



j*(t)= dr((j a (t);H ext {r))), (8.70) 

J — oo 

where 

H ext (T) = Hl(r) cos wre £T , (8.71) 



Here j a is the current density operator, if the volume of the system is taken to be unity. The 
equation (|8.70p can be transformed to the following form 



j a (t) = -Re{ / dr(0«(t); Hl( T )))— — - + (8.72) 



(L 7a (o), J ff T 1 (o)]_) e ^+ £t -i— }. 



U) — l£ ' 
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Noting that 

Hlir) = -(Sj(r)), [r?,rj] = r^8 aP 6 ij , (8.73) 

we get from this equation 

j a (t) = Re{a a ^(u)£p exp(iujt + et)}, (8.74) 

where 

(Ta0(p) = <W + / dr (( Ja );# r ) . , - 8.75 

raw J-oo ia; + e 

is the conductivity tensor, and n the number of electrons per unit volume. The first term in 
Eq. (|8.75p corresponds to the electrical conductivity of a system of free charges and is not connected 
with the inter-particle interaction. As oj — > oo the second term decreases more strongly than the 
first one 

2 

en 

lim Imwc7 Q/3 (w) = 5 a /3, 

and the system behaves as a collection of free charges. 

Let us discuss the derivation of the conductivity tensor in general form. Consider a system of 
charged particles in electrical field £, which is directed along the axis (3 (J3 = x,y,z). The 
corresponding electrostatic potential ipw 

has the following form 



(r,t) =J2yJ2 V W (?> n ) ex P (i nt + i£t)/ih - iq^j • (8.76) 



y 

n q 



In the momentum representation the above expression reads 

SjfXq, SI) = % -q a ^\q, SI) = 5 af3 £p(q, SI). (8.77) 
Consider the case when the perturbation H ext has the form 

H fr = ^*P m & exp((Ot + iet)/ih) . (8.78) 

Here 

Et t _ 

a p+q a Pi Vq — V—q- 

V 

is the particle density operator. 

A reaction of the system on the perturbation is given by 

lW(q,Sl) = A(ej a (q)) = ^((ej a (q)\e^\q,Sl)4))n +i s (8.79) 

i r 00 e^l® — i 

= e 2 ^(q,n)- dcoJ(r)l,j a (q);co) 



V J_oo y Sl — UJ + lE 

where 1^ (q, SI) denotes the component of the density of the current in the a-direction when 
external electric field directed along the /3-axis. 
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8.4 Linear response theory: pro et contra 

It was shown in the previous sections that the formulation of the linear response theory can be 
generalized so as to be applied to a rather wide class of the problems. It is worth to note that the 
"exact" linear expression for electrical conductivity for an arbitrary system was derived originally 
by KubcP^in a slightly different form than Eq. (j8.75p 



1 



oo 



<V = lim - ^(0) + / dte- £t <f)^{t) , (8.80) 

e— ^ e V J — oo / 

where ^ 

<?V(i) = ^Tr (n, ^ e^x^ ^ ei± iu (t) (8.81) 

i i 

is the current response in the \i— direction when a pulse of electric field £{t) is applied in the 
v— direction at t = 0; e» is the charge of the ith particle with position vector fj and n is a density 
operator. In the one-electron approximation Eq. (l8.80p reduces to (cf. RefP^) 

a l3 = 27re 2 h^(m\v l \n)(n\v J \m) . (8.82) 

nm ^ ' n 

Here Vi is the velocity operator and / is the Fermi function. 

To clarify the general consideration of the above sections it is of interest to consider here a 
simplified derivation of this formula, using only the lowest order of time-dependent perturbation 
theoryPSl This approach is rooted in the method of derivation of Callen and WeltonpSS Dl rt 
takes explicitly into account degeneracy of the states. The linear response theory formulated by 
KubcP^HI was motivated in parlP^ES] D y the prior work of Callen and WeltonpSSl wriQ proposed 
a quantum mechanical perturbative calculation with the external forces exerted on a dissipative 
system. They pointed out a general relationship between the power dissipation induced by the 
perturbation and the average of a squared fluctuation of the current of the system in thermal 
equilibrium. 

Let us discuss first the role of dissipation. A system may be called to be dissipative if it absorbs 
energy when subjected to a time-periodic perturbation, and linear if the dissipation (rate of 
absorption of energy) is quadratic in the perturbation. For a linear system an impedance may be 
defined and the proportionality constant between the power and the square of the perturbation 
amplitude is simply related to the impedance. In the case of electrical current in a material one 
can write down that 



W = ~ - ^ = V 



2 



2 (Impedance) 2 

where W is the average power and V is the voltage. If we calculate the power microscopically in 
some way and find it quadratic in the applied force (voltage), then comparison with this equation 
will give the conductivity of the substance. 

Consider a situation when an electron of charge e is situated a distance x from the end of a resistor 
of length L and then a voltage V = Vosinwt is applied in the x direction! The perturbation 
term in the Hamiltonian will be of the form 



H* xt = V ey sin ut. (8.83) 
Li 

The Hamiltonian of a system in the absence of the field (but including all other interactions) is 
denoted by Hq with corresponding wave function ip n such that 

H ^ n = E n ^ n . (8.84) 
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The total wave function may be expanded, in terms of the ip n 

*n = ^ a "(#™. ( 8 - 85 ) 

n 

where the coefficients a n (t) may be approximately determined by first-order perturbation theory 
The rate of transition is then given by 

= E K"»|x|n)| 2 [5(£? m - (£ n + fej)) + 5(E m - (E n - hto))]. (8.86) 

mn 

The first term corresponds to a transition to a state 2£ m = E n + hui in which energy foj is absorbed, 
whereas the second term corresponds to a transition to a state E m = E n — Huj in which energy is 
emitted. Hence, the net rate of absorption of energy is given by 

dK = l^^l^y, \( m \x\n)\ 2 [5(E m - E n - frw) + 5(E m - E n + fiw)], (8.87) 

ran 

which is quadratic in Vo- 

This equation gives the absorption rate for a single, isolated electrons, but in a real system we 
are dealing with an ensemble of these, which we shall represent by the Fermi function. One must 
therefore fincP^ the average absorption by averaging over all initial states | n) and taking the Pauli 
exclusion principle as well as the two spin directions into account. The result is^^l 

(^)= 7 re 2 y 2 -ElHx|n)| 2 (8.88) 

mn 

{f(E n )5(E m - E n -hw)(l- f{E m )) 
-f(E n )5(E m -E n + MU " f(E m ))} 
= TreVfM^ \(n + hw\x\n)\ 2 [f(E n ) - f(E n )f(E n + fiw)] - 

n 

£|(n'|z|n' + hw)\ 2 [f(E n , + huj)- f(E n , + hu)f(E n ,)]}, 

n> 

where we have put n' = n — fkj. The above formula can be transformed to the form 

(^) = vre 2 y 2 W ^|(mNn)| 2 . (8.89) 

mn 

[f (En) - f(E m )]5(E m -E n - M- 
This expression may be simplified by introducing the matrix element of the velocity operator 

^ = {[H,x]. (8.90) 

In principle, on the right-hand side of Eq. (|8.90p the total Hamiltonian 

H = H + H ext (8.91) 

should be written. In this case, however, the terms of higher than quadratic order in Vo appear. 
For a linear system one can neglect these and use Eq. (l8.90p . Thus we have 

% i 

\(m\x\n)\ = —(m\[H, x]\n) = —(E m (m\x\n) — E n {m\x\n)) (8.92) 
n n 
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and 



Therefore 



\(m\x\n)\ 2 = -T^-{E m - E n ) 2 \(m\x\n)\ 2 = uj 2 \(m\x\n)\ 2 . (8.93) 

fir 

(~~TT~) = ~—V 2 V \(m\x\n)\ 2 [f(E n ) - f(E n + fiw)]5(E m - E n - fiu). (8.94) 
dt 'jJ 

van 

If we now assume the current to be in phase with the applied voltage, the average energy dissipation 
becomes 

1 V 2 

Power = f- (8.95) 

2 R(u) 

as the resistance R(oj) is now equal to the impedance Z(lj). Referring everything to unit volume, 
and noting that the resistance per unit volume is the resistivity, we get for the conductivity a 

&xx{u) = \(m\v \n)(n\v\m)\ (8.96) 

id — J 
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If (En) - f(E n + fiw)]8(E m -En-hu). 
A straightforward generalization of this procedure, using a perturbation 

H ext = Voexi/L sm{ujt) , (8.97) 



leads to the definition of an impedance matrix and a conductivity tensor 

• |Hwil™}(nM™>l[/( s n) - f( E n + tuMEm -E n -Hu), (8. 
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which is the Kubo-Greenwood equationP^ESl The derivatiorPSSl presented here confirms that 
fact that the linear response theory is based on the fluctuation-dissipation theorem, i.e. that 
the responses to an external perturbation are essentially determined by fluctuations of relevant 
physical quantities realized in the absence of the perturbation. Thus the linear response theory 
has a special appeal since it deals directly with the quantum mechanical motion of a process. 
The linear response theory (or its equivalent) become soon a very popular tool of the transport 
theory ES32S1] ^ s was expressed by Langerj^^the Kubo formula " probably provides the most rig- 
orous possible point of departure for transport theory. Despite its extremely formal appearance, it 
has in fact proved amenable to direct evaluation for some simple models" . Edwards 26 ^ and Chester 
and Thellung^H have used the Kubo formula to calculate the impurity resistance of a system of 
independent electrons, and have recovered the usual solution of the linearized Boltzmann equa- 
tion. VerboverP^Sl hag extended this work to higher orders in the concentration of impurities and 
has found corrections to the conductivity not originally derived via Boltzmann techniques. It was 
concluded that the Kubo formula might be most fruitfully applied in the full many-body problem, 
where it is not clear that any Boltzmann formulation is valid. However Izuyama^^ casted doubt 
on the Kubo formula for electrical conductivity. He claimed that it is not, in fact, an exact formula 
for electrical conductivity, but is rather a coefficient relating current to an "external" field, which 
coefficient is equal to the conductivity only in special case. A correlation-function formula for 
electrical conductivity was derived by MagarP^Sl by a formalism which gives prominence to the 
total electric field, including fields which may arise from the charged particles which are part of 
the sy stem being studied. 

Langer^^ have evaluated the impurity resistance of an interacting electron gas on the basis of 
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the Kubo formula at low but finite temperatures. The calculations are exact to all orders in 
the electron-electron interactions and to lowest order in the concentration of impurities. In the 
previous papers the impurity resistance of this gas was computed at absolute zero tem- 

perature. It was 

showiPS that the zero-temperature limit of this calculation yields the previous 
result. In RefJ^^SI Kubo formula for thermal conductivity was evaluated for the case of an inter- 
acting electron gas and random, fixed, impurities. The heat flux was examined in some detail and 
it was shown that in a normal system where the many-body correlations are sufficiently weak, 
the Wiedemann-Franz law remains valid. The relationships between the transport coefficients 
of a metal in a strong magnetic field and at very low temperatures were discussed by Smrcka 
and StredaJ2221 Formulae describing the electron coefficients as functions of the conductivity were 
derived on the basis of the linear response theory. As was mentioned earlier, examples of such 
relations are the Wiedemann-Franz law for the heat conductivity k and the Mott rule^^ for the 
thermopower S. It was shown that that the Wiedemann-Franz law and the Mott rule are obeyed 
even in the presence of a quantized magnetic field ujt > 1 if the scattering of electrons is elastic 
and if Hui 3> kT. 

A theoretical analysis, based on Kubo formalism, was made for the ferromagnetic Hall effect by 
Leribaux^ 4 in the case of transport limited by electron-phonon scattering. The antisymmetric, 
off-diagonal conductivity was, to first order in magnetization, found to be of order zero in the 
electron-phonon interaction (assumed to be weak) and, to this order, was equivalent to Karplus 
and Luttinger resultsP^ Tanaka, Moorjani and MoritsP^ expressed the nonlinear transport coef- 
ficients in terms of many-time Green functions and made an attempt to calculate the higher-order 
transport coefficients. They applied their theory to the calculation of the nonlinear susceptibil- 
ity of a Heisenberg ferromagnet and nonlinear polarizability problem. Schotte^^teconsidered the 
linear response theory to show the closeness of it and kinetic equations. 

At the same time, several authors^^SQI ra i sec [ an important question as to the general validity 
of the correlation formulae for transport coefficients. The relation for the electrical conductivity 
was, in principle, not questioned as in this case (as was shown above) one may obtain it by a 
straightforward application of basic statistical mechanical principles and perturbation techniques. 
One simply calculates the response of the system to an electric field. What has been questioned, 
however, was the validity of the correlation relations for the transport coefficients of a fluid - the 
diffusion constant, the thermal conductivity, and the viscosity - where no such straightforward 
procedure as was used for the electrical conductivity was availablepSD 

However, Jackson and Mazur^^ presented a derivation of the correlation formula for the viscosity 
which is similar in spirit, and as free of additional assumptions, as that for the electrical conduc- 
tivity. The correlation formula for the viscosity was obtained by calculating statistically the first 
order response of a fluid, initially in equilibrium, to an external shearing force. It was shown that 
on the basis of this derivation, the correlation formula for the viscosity, the exactness of which 
had been questioned, was placed on as firm a theoretical basis as the Kubo relation for the electri- 
cal conductivity. In addition, Resibois^^" 2 ^ demonstrated the complete equivalence between the 
kinetic approach developed by Prigogine and coworkers 2 ^! and the correlation function formal- 
ism for the calculation of linear thermal transport coefficients. It was shown that in both cases 
these transport coefficients are determined by the solution of an inhomogeneous integral equation 
for a one-particle distribution function which is the generalization to strongly coupled system of 
the Chapman-Enskog first approximation of the Boltzmann equation! 45 ^ 286 ! Interesting remarks 
concerning the comparison of the linear response theory and Boltzmann equation approach were 
formulated by R. PeierlsP^E^ SchofielcP^^ elaborated a general derivation of the transport 
coefficients and thermal equilibrium correlation functions for a classical system having arbitrary 
number of microscopic conservation laws. This derivation gives both the structure of the correla- 
tion functions in the hydrodynamic (long wavelength) region and a generalized definition of the 
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transport coefficients for all wave-lengths and frequencies. 

A number of authors have given the formulation of nonlinear responses j25M29QH295| wag s ] lown 
that since in a nonlinear system fluctuation sources and transport coefficients may considerably 
depend on a nonequilibrium state of the system, nonlinear nonequilibrium thermodynamics should 
be a stochastic theory. From the other side the linearity of the theory itself was a source of many 
doubts. The most serious criticism of the Kubo linear response theory was formulated by van 
KampenP2Sl He argued strenuously that the standard derivation of the response functions are 
incorrect. In his own words, " the basic linearity assumption of linear response theory ... is 
completely unrealistic and incompatible with basic ideas of statistical mechanics of irreversible 
processes" . The main question raised by van Kampen concerned the logic of the linear response 
theory, not the results. As van KampenP^ expressed it, "the task of statistical mechanics is not 
only to provide an expression for this (transport) coefficient in term of molecular quantities, but 
also to understand how such a linear dependence arises from the microscopic equations of motion" . 
The van Kampen's objections'^ to Kubo linear response theory can be reduced to the following 
points. In the linear response theory one solves the Liouville equation to first order in the ex- 
ternal field E (electric field). This is practically equivalent to following a perturbed trajectory 
in phase space in a vicinity of the order of \E\ of unperturbed trajectory. In the classical case 
trajectories are exponentially unstable and corresponding field E should be very small. Kubo's 
derivation supposed non-explicitly that macroscopic linearity (Ohm law, etc.) is the consequence 
of microscopic linearity, but these two notions are not identical. Macroscopic linearity is the re- 
sult of averaging over many trajectories and is not the same as linear deviations from any one 
trajectory. In other words, van Kampen's argument was based on the observation that due to the 
Lyapunov instability of phase-space trajectories, even a very small external field will rapidly drive 
any trajectory far away from the corresponding trajectory without field. Hence linear response 
theory which is based on the proportionality of the trajectory separation with the external field 
could only be expected to hold for extremely short times, of no physical interest. 
Responses to van Kampen's objections were given by many 

authorsPEHESZISIIl The mam argu- 
ments of these responses, were based on the deep analysis of the statistical mechanical behavior of 
the many-body system under the external perturbation. It was shown that in statistical mechani- 
cal calculations one deals with the probability distributions for the behavior of the many particles 
rather than to the behavior of an individual particle. An analysis of the structural stability of 
hyperbolic dynamics-averaging and some other aspects of the dynamical behavior shows that the 
linear separation of trajectories goes on long enough for Green-Kubo integrals to decay. Moreover, 
Naudts, Pule and Verbeure^^l analyzed the long-time behavior of correlations between extensive 
variables for spin-lattice systems and showed that the Kubo formula, expressing the relaxation 
function in terms of of the linear response function, is exact in the thermodynamic limit. 
It was mentioned above that in the 1950s and 1960s the fluctuation relations, the so-called Green- 
Kubo relations SSES1I3Q3H305] were derived for the causal transport coefficients that are defined by 
causal linear constitutive relations such as Fourier law of heat flow or Newton law of viscosity. 
Later it was shown also that it was possible to derive an exact expression for linear transport 
coefficients which is valid for systems of arbitrary temperature, T, and density. The Green-Kubo 
relations give exact mathematical expression for transport coefficients in terms of integrals of time 
correlation functions HSHSEQBHHQS] More precisely, it was shown that linear transport coefficients 
are exactly related to the time dependence of equilibrium fluctuations in the conjugate flux. For 
a more detailed discussion of these questions see Ref s HSHHDI2Q2H2Q5] 

To summarize, close to equilibrium, linear response theory and linear irreversible thermodynam- 
ics provide a relatively complete treatment. However, in systems where local thermodynamic 
equilibrium has broken down, and thermodynamic properties are not the same local functions 
of thermodynamic state variables such that they are at equilibrium, serious problems may ap- 
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9 The Nonequilibrium Statistical Operator Method and Kinetic 
Equations 

The method of the nonequilibrium statistical operatoJ^ ^ 211 ^^^^ was reviewed already in de- 
tails 4 - by us. In this section, we remind very briefly the main ideas of the nonequilibrium statistical 
operator (NSO) ^proactfMmmMMmWBmi f or the sake of a self-contained formulation. The 
precise definition of the nonequilibrium state is quite difficult and complicated, and is not uniquely 
specified. Thus the method of reducing the number of relevant variables was proposed. A large 
and important class of transport processes can reasonably be modelled in terms of a reduced 
number of macroscopic relevant variables JSBS5] j^. was SU pp 0sec l that the equations of motion 
for the "relevant" variables (the space- and time-dependent thermodynamic "coordinates" of a 
many-body nonequilibrium system), can be derived directly from the Liouville equation. This 
can be done by defining a generalized canonical density operator depending only upon present 
values of the thermodynamic "coordinates". According to D. N. ZubarevJ^ED the NSO method 
permits one to generalize the Gibbs ensemble method to the nonequilibrium case and to con- 
struct a nonequilibrium statistical operator which enables one to obtain the transport equations 
and calculate the kinetic coefficients in terms of correlation functions, and which, in the case of 
equilibrium, goes over to the Gibbs distribution. The basic hypothesis is that after small time- 
interval t the nonequilibrium distribution is established. Moreover, it is supposed that it is weakly 
time-dependent by means of its parameter only. Then the statistical operator p for t > r can be 
considered as an " integral of motion" of the quantum Liouville equation 

f + i[p,ff] = 0. (9.1) 

Here dp/dt denotes time differentiation with respect to the time variable on which the relevant 
parameters F m depend. It is important to note once again that p depends on t by means of F m (t) 
only. We may consider that the system is in thermal, material, and mechanical contact with a 
combination of thermal baths and reservoirs maintaining the given distribution of parameters F m . 
For example, it can be the densities of energy, momentum, and particle number for the system 
which is macroscopically defined by given fields of temperature, chemical potential and velocity. 
It is assumed that the chosen set of parameters is sufficient to characterize macroscopically the 
state of the system. The set of the relevant parameters are dictated by the external conditions 
for the system under consideration and, therefore, the term dp/dt appears as the result of the 
external influence upon the system; this influence causes that the system is non-stationary. 
In order to describe the nonequilibrium process, it is also necessary to choose the reduced set of 
relevant operators P m , where m is the index (continuous or discrete). In the quantum case, all 
operators are considered to be in the Heisenberg representation 

fiHt\ f-iHt\ , , 

P m (t) = expf— JP m exp( — — J, (9.2) 

where H does not depend on the time. The relevant operators may be scalars or vectors. The 
equations of motions for P m will lead to the suitable "evolution equations" H2E51 j n the quantum 

case 

2^-I[P ra (t),fl] = 0. (9.3) 

The time argument of the operator P m (t) denotes the Heisenberg representation with the Hamil- 
tonian H independent of time. Then we suppose that the state of the ensemble is described by a 
nonequilibrium statistical operator which is a functional of P m (t) 



56 



p(t) = p{...P m (t)...}. (9.4) 

Then p{t) satisfies the Liouville equation (19, lh . Hence the quasi-equilibrium (local-equilibrium) 
Gibbs-type distribution will have the form 

p q = Q q l exp ^- Fm{t)P^)j , (9.5) 

where the parameters F m (t) have the sense of time-dependent thermodynamic parameters, e.g., 
of temperature, chemical potential, and velocity (for the hydrodynamic stage), or the occupation 
numbers of one-particle states (for the kinetic stage). The statistical functional Q q is defined by 
demanding that the operator p q be normalized and equal to 



Q q = Tr exp ^- ^ F m (i)P m j . 



(9.6) 



The kinetic equations are of great interest in the theory of transport processes. In the NSO 
approachpKSEEl main quantities involved are the following thermodynamically conjugate 
values: 

(P m )=--^L-; Fm ( t ) = -?jL. (9.7) 
oF m (t) o{P m ) 

The generalized transport equations which describe the time evolution of variables (P m ) an d 
F m follow from the equation of motion for the P m , averaged with the nonequilibrium statistical 
operator (|9.4p . It reads 



10 Generalized Kinetic Equations and Electroconductivity 
10.1 Basic formulas 

Let us consider a many-particle system in the quasi-equilibrium state. It is determined completely 
by the quasi-integrals of motion which are the internal parameters of the system. In this and 
following sections we will use the notation Aj for the relevant observables to distinguish it from 
the momentum operator P. Here, for the sake of simplicity, we shall mainly treat the simplest case 
of mechanical perturbations acting on the system. The total Hamiltonian of the system under the 
influence of homogeneous external perturbation, depending on time as ~ exp(iujt) is written in 
the following form 

H(t) = H + H F (t), H F (t) = - AjFj exp(iut). (10.1) 

i 

In the standard approach the statistical operator p can be considered as an "integral of motion" 
of the quantum Liouville equation 

^ + ^[p(t),n(t)] = o. (io.2) 

Using the ideas of the method of the nonequilibrium statistical operator, as it was described above, 
we can write 

p(t) = e f dhe'^-^Uit^^pihp^t^i). (10.3) 
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The time-evolution operator U(t,t\) satisfy the conditions 

d U (t,t 1 ) = l-H(t)U(t,t 1 ), 
ot in 

^u(t, h) = -±n{t)u(t, h), u(t, t) = 1. 

oti in 

If we consider the special case in which p{t\) — > po, where po is an equilibrium solution of the 
quantum Liouville equation (110. 2D . then we can find the nonequilibrium statistical operator from 
the following equation: 

^fr + h^ m] = - £{p * ~ po) - (ia4) 

In the lim e _ > .Q+ the nonequilibrium statistical operator will corresponds to the Kubo density matrix 
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pf {t) = Po = ~^ I dhe-^-^U^h^HFit^pH^h). (10.5) 



]imTr(p?(t)B j ) = (B j 

E-J-0+ 



An average of the observable Bj are defined as 

The linear response approximation for the statistical operator is given by 

-^^AjFjexpiiut) J dt ie ( £+ ^exp(^i) [p , Aj] exp ( Z! f^) • (10-6) 

Using the obtained expression for the statistical operator the mean values of the relevant observ- 
ables Bi can be calculated. To simplify notation only observables will be considered for which the 
mean value in the thermal equilibrium vanishes. In other words, in general the Bi will be replaced 
by Bi - (B,) , where (BJq = Tr{p G Bi). We find 

(Bi) t = Tr(pK(t)Bi) =J2 L ij(") F j expM), (10.7) 

j 

where the linear response coefficients (linear admittances) are given by 

Lij(u) = (ij-s Si}*-*, (e + ). (10.8) 

This expression vanish for operators Aj commuting with the Hamiltonian H of the system, i.e. 
the Kubo expressions for Lij(uo) vanish for all uj where for ui = the correct result is given by 

(B^ =PJ2 FjTt( Po A j B 1 ) = UM F i expM). (10.9) 

3 j 

The scheme based on the NSO approach starts with the generalized Liouville equation 

^ + l[p £ (t),n(t)} = -e{Pe(t) - Pq {t))- (io.io) 
For the set of the relevant operators P m it is follow that 

^-Tv( Pe (t)P m ) = 0. (10.11) 
at 
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Here notation are 

p q (t) = Q~ l exp[-P(H -fi.N-^2 F rn (t)P rn )}, (10.12) 

m 

F m (t) = F m exp(wi), P m ^ P m - {P m )o- 

To find the approximate evolution equations, the p e (t) can be linearized with respect to the 
external fields Fj and parameters F m 

dtie (e+iw) * 1 / dr (10.13) 

-oo JO 

[- ^2 FjMt ~ iT )+Yl F m p m (t - ir) +i^Y. F m p m(t - ir)]) . 

j m m 

As a result we find 

Y,Fn ((P n ;P m )u-ie + iu(Pn;Pm)u-ie) = ^F j {A j - P m )^ ie . (10.14) 

n j 

In other notation we obtain 

Pm) ~\~ {Pfi'i Pm) co— is 

+ iuj[(P n \Pm) + {P n ;P m ) uj— is ] (10.15) 

n 

Y. 'l {i A j\ P m) + (Af,Pn 
j 

where 

/oo 
dte^-^\A(t)\B); (10.16) 
-oo 

(A\B) = f dXTr(p A(-i\)B), p = -e^; 
Jo V 

Xij = {Ai\Aj), ^ XikU^kj = 

k 

Thus our generalized transport equation can be written in the following abbreviated form 

F nPnm = ^ F 3 K 3™- ( 10 - 17 ) 

n j 

10.2 Electrical conductivity 

The general formalism of the nonequilibrium statistical operator has been the starting point of 
many calculations of transport coefficients in concrete physical systems. In the present section 
we consider some selected aspects of the theory of electron conductivity in transition metals and 
disordered alloy sGSHMloi yy e p U ^ /j = i f or simplicity of notation. Let us consider the dc electrical 
conductivity 

a= 3^ {P]P) = 3^f dte£t l PdXTI{pP{t ~ iX)P), (10 ' 18) 
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where P is the total momentum of the electrons. Representing P as a sum of operators Ai (relevant 
observables) which can be chosen properly to describe the system considered (see below), the 
corresponding correlation functions (Ai;Aj) can be calculated by the set of equations^^ 

J2(ieSik - iJ2( A l\ A l)(x~ X )lk + i^2u u ( yX ' 1 )ik){A k ;A j ) =i Xij , (10.19) 



where 



Xij = {Ai\Aj), X l k{x~ 1 )kj = Sij] (10.20) 



k 

rP i 

(Ai\Aj)= / dXTr(pAi(—i\)Aj), p=-e^ H ; (10.21) 
Jo v 

Uij = {AiClCAjf, A t = i[H,A t ]. (10.22) 

The operator C = 1 — P is a projection operator^® with P = and (. . -) c 

denotes that in the time evolution of this correlation function, operator L = i[H,...] is to be 
replaced by CLC. By solving the set of equations (|10.19|) the correlation functions (A^Aj) which 
we started from are replaced by the correlation functions Uij f|10.22[) . and with a proper choice of 
the relevant observables these correlation functions can be calculated in a fairly simple approxi- 
mationP^ n however, difficult to go beyond this first approximation and, in particular, to take 
into account the projection operators. Furthermore this method is restricted to the calculation 
of transport coefficients where the exact linear response expressions are known; generalization to 
thermal transport coefficients is not trivial. In RefP^ we described a kind of general formalism 
for the calculation of transport coefficients which includes the approaches discussed above and 
which can be adapted to the problem investigated. 

For simplicity of notation we restrict our consideration here on the influence of a stationary ex- 
ternal electrical field. In the linear response theory the density matrix of the system becomes 

p LR = — [ dt'e^'-V [ d\pP(t' -t-i\), e^0, (10.23) 
m J-oo Jo 

where the time dependence in P(t) is given by the total Hamiltonian of the system without the 
interaction term with the electrical field. From another point of view we can say that there is 
a reaction of the system on the external field which can be described by relevant observables 
such as shift of the Fermi body or a redistribution of the single particle occupation numbers, etc. 
Hence, for small external fields the system can be described in a fairly good approximation by the 
quasi-equilibrium statistical operator of the form 

p q = _L e -^ H -£* aiA (10.24) 

Qq 

where the Ai are the observables relevant for the reaction of the system and the «j are parameters 
proportional to the external field. Of course, the statistical operator (|10.24|) is not a solution 
of the Liouville equation, but an exact solution can be found starting from (|10.24p as an initial 
condition: 

/o 
dt'e^* -t) exp(iH a tf - t))p q exp(-iH s (t' - t)), (10.25) 

where H s = H — eE^-r^. In order to determine the parameters «j we demand that the mean 
values of the relevant observables A^ are equal in the quasi-equilibrium state p q and in the real 
state p s , i.e. 

Tr(p s Ai) = Tr(p q Ai). (10.26) 
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This condition is equivalent to the stationarity condition 



^(MO = o. 

For a sufficiently complete set of operators Ai the condition (|10.26j) ensures that p q describes the 
system with a sufficient accuracy. Linearizing Eq. (ll0.26p with respect to the parameters ai and 
the external field E, one obtains the set of equations 

Y^Oj^-iTrWAjtAi]) + = ^((Mi) + <^>)- (10-27) 

j 

This set of equations can be shown to be equivalent to Eq. (|10.19p ; whereas in the higher orders 
of interaction the equations (|10.27|) are more convenient to handle because the time dependence 
is given here in an explicit form without any projection. 
With the parameters at the current density is given by 

j 

Supposing that the total momentum P of the electrons can be built up by the operators Ai, it can 
be shown that the Eqs. (|10.27p and (|10.28[) include the Kubo expression for the conductivity (see 
below) . In order to solve the system of equations (I10.27P a generalized variational principle can be 
formulated, but the reduction of the number of parameters aj by a variational ansatz corresponds 
to a new restricted choice of the relevant observable A4. 



11 Resistivity of Transition Metal with Non-spherical Fermi Sur- 
face 

The applicability of the transport equations (110. 27j) and (I10.28P derived above to a given problem 
depends strongly on the choice of the relevant operators Aj. The first condition to be fulfilled 
is that the mean values of the occupation numbers of all quasiparticles involved in the transport 
process should be time independent, i.e. 

Ti(p s (ai)n k ) = Tr(p q (ai)n k ), —Tr(p s (ai)n k ) = 0. (11.1) 

Of course, this condition is fulfilled trivially for Ai — > n k , but in this case the equations (|10.27p 
cannot be solved in practice. In the most cases, however, there can be found a reduced set 
of operators sufficiently complete to describe the reaction of the system on the external field. 
It can be shown that under certain conditions the scattering process can be described by one 
relaxation time only and then the set of relevant operators reduces to the total momentum of 
the electrons describing a homogeneous shift of the Fermi body in the Bloch (momentum) space. 
These conditions are fulfilled for a spherical Fermi body at temperatures small in comparison 
to the degeneration temperature and for an isotropic scattering mechanism where the scatterers 
remain in the thermal equilibrium. In this simple case the Eqs. (ll0.27p are reduced to the so-called 
resistivity formulaP^ For non-spherical Fermi bodies the set of relevant observables has to be 
extended in order to take into account not only its shift in the /c-space but also its deformation. 
Our aim is to develop a theory of electron conductivity for the one-band model of a transition 
metal. The model considered is the modified Hubbard model, which include the electron-electron 
interaction as well as the electron-lattice interaction within MTBA. The non-spherical Fermi 
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surface is taken into account. 

The studies of the electrical resistivity of many transition metals revealed some peculiarities. It 
is believed that these specific features are caused by the fairly complicated dispersion law of the 
carriers and the existence of the two subsystems of the electrons, namely the broad s — p band and 
relatively narrow d-band. The Fermi surface of transition metals is far from the spherical form. In 
addition, the lattice dynamics and dispersion relations of the phonons is much more complicated 
than in simple metals. As a result, it is difficult to attribute the observed temperature dependence 
of the resistivity of transition metals to definite scattering mechanisms. 

Here we investigate the influence of the electron dispersion on the electrical resistivity within a 
simplified but workable model. We consider an effective single-band model of transition metal with 
tight-binding dispersion relation of the electrons. We take into account the electron-electron and 
electron-lattice interactions within the extended Hubbard model. The electron-lattice interaction 
is described within the MTBA. For the calculation of the electrical conductivity the generalized 
kinetic equations are used which were derived by the nonequilibrium statistical operator method. 
In these kinetic equations the shift of the non-spherical Fermi surface and its deformation by the 
external electrical field are taken into account explicitly. By using the weak scattering limit the 
explicit expressions for the electrical resistivity are obtained and its temperature dependence is 
estimated. 

We consider a transition metal model with one non-spherical Fermi body shifted in the fc-space 
and deformed by the external electrical field. Hence, the Fermi surface equation E(k) = Ep is 
transformed into E{k) = Ep, where 

E{k) = E{k) + mi/i^ + mV + . . . (11.2) 

dk i=Xa dk 

The term proportional to v\ describes a homogeneous shift of the Fermi surface in the A;-space 
and the last terms allow for deformations of the Fermi body. The polynomials have to be 

chosen corresponding to the symmetry of the cry st a pHDl£L!ML2l anc i j n consequence of the equality 
E{k + G) = E{k) (which is fulfilled in our tight-binding model) they should satisfy the relation 

&(k + G) = 

where G is a reciprocal lattice vector. Thus the relevant operators are given by 

dE 

Ai->mJ2 = 1. (H.3) 

For our tight-binding model we are restricting by the assumption that <J> 1 — > for i > 3. 
11.1 Generalized kinetic equations 

The quantum-mechanical many-body system is described by the statistical operator p s obeying 
the modified Liouville equation of motion Eq. tllO.l'Ol) 

^± + ^[ Ps ,H s ] = -e(p s -p q ), (11.4) 

where H s = H + Hp is the total Hamiltonian of the system including the interaction with an 
external electrical field Hp = —eE^i^i and p q is the quasi-equilibrium statistical operator. 
According to the NSO formalism the relevant operators P m should be selected. These operators 
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include all the relevant observables which describe the reaction of the system on the external 
electrical field. These relevant operators satisfy the condition 

Tr( Ps (t, 0)P m ) = (P m ) = (P m ) q ; Tip s = 1. (11.5) 

This condition is equivalent to the stationarity condition 

£(P m )=0^{[H s ,P m ]) + {[H E ,P m ]) = 0. (11.6) 
at 

In the framework of the linear response theory the operators p s and p q in Eqs. (|11.5|) and (|11.6|) 
should be expanded to the first order in the external electrical field E and in the parameters F m . 
Thus Eq. (fTTBI) becomes® 3 

^F n (-iTr(p[P n ,P m ]) + (P n ;P m )) = ^-(Tv(pP(-iX)P m ) + (P;P m >), (11.7) 

n 

where 

P m = i[H s ,P m ], (11.8) 
{A;B)= dte et d\Tr(pA(t - iX)B), (11.9) 

J-oo JO 

A(t) = exp(iHt)Aexp(-iHt); p = Q' 1 exp(-/3H), 
and P is the total momentum of the electrons. 

The equations (I11.7P are generalized kinetic equations in which the relaxation times and particle 
numbers are expressed via the correlation functions. It will be shown below that these equations 
can be reduced to the Kubo formula for the electrical conductivity provided a relation P e = ^ otiPi 
can be found. 

The generalized kinetic equations (jll.7[) can be solved and the parameters F m can be determined 
by using a variational principle. The current density is given by 

3 = -^Tv( Ps P e ) = -^Tr(p q P e ), (11.10) 
mil mil 

where the conditions (|ll,5p have been used and £1 is the volume of the system. Linearizing 
Eq. (|ll,10p in the parameter F m we find 

j = — V F m / dXTv(pP m (t - iX)P e ) = -E, (11.11) 

m JU 

where the proportionality of the F m to the external electrical field has been taken into account. 
For the tight-binding model Hamiltonian (|4.65p . (|4.66p the proper set of operators P m is given by 

dE 

P e = P 1=m ^2—ala ka , (11.12) 

ka 

B F 

P i = m ^&0)—al a a ka , (i = 2...n). (11.13) 

kcr 

The parameters F m are replaced by the generalized drift velocities. Then the quasi-equilibrium 
statistical operator p q take the form 

Pq = ^^v{-P[H + rnv 1 ^—a\ a a ka + mv^ l {k)— a\ a a kcr y) . (11.14) 



ka i=2... n 
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It should be mentioned here that in general in p q the redistribution of the scatterers by collisions 
with electrons should be taken into consideration. For the electron-phonon problem, e.g., the 
phonon drag can be described by an additional term v p hP p h in Eq. (|11.2j) where v p h is the mean 
drift velocity and P p h is the total momentum of the phonons. Here it will be supposed for 
simplicity, that due to phonon-phonon Umklapp processes, etc., the phonon subsystem remains 
near thermal equilibrium. In the same way the above consideration can be generalized to many- 
band case. In the last case, the additional terms in (|11.2p . describing shift and deformation of 
other Fermi bodies should be taken into account. 

For the tight-binding model Hamiltonian (|4.65p . (|4.66p the time derivatives of the generalized 
momenta (generalized forces) in Eq. (111.71) are given by 

/''„ • /' ; Pf • Pf, (11.15) 

See 



P J ee = i[H ee ,P j }= (11.16) 
iUm v-^ v-^ (*i„-'\dE xi „->>9fi ^in^^^E\ 

feiK2 k^k^G 

a I 1 t afc 2t a fc 3 4. a fc4l^(^i -k 2 + k 3 -k A + G), 
/ B F 1 B F 1 \ 

p f = m EEEi4 $3 ^fe - ^^dkd ■ (1L17) 

fcifc 2 qG uv 1 1 

% 2 a a kiA b lv + b- qu )S(k 2 -h + q + G). 

We confine ourselves by the weak scattering limit. For this case the total Hamiltonian H in 
evolution A(t) = exp(iHt)Aexp(—iHt) we replace by the Hamiltonian of the free quasiparticles 
H® + Hf. With this approximation the correlation functions in Eq. (|11.7j) can be calculated 
straightforwardly. We find that 

(Pj-Pi) « <Pf;if e ) + (Pf;Pf). (11.18) 

Restricting ourselves for simplicity to a cubic system, the correlation functions of the generalized 
forces are given by 

(P J ee ;Pn= (11-19) 

U 2 m 2 pTT 



/l Aj(k 1 ,k2,k 3 ,ki)Ai{ki,k 2 ,k3,ki) 



where 



N 2 

k\ki k:>,k^G 

" /fe 2 )/fc 3 (l " fkME(h) - E(k 2 ) + E(k 3 ) - E(k 4 )) 

5(k x -k 2 + k 3 -k4 + G), 

{Pf;Pf) = (11.20) 

2vrm 2 /3 £ 5>jU 2 ) 2 ^(*i> fe 2 )^(fci, fe) • 
fcifc 2 g^G 

/fa(l " f kl )N(qv)5(E(k 2 ) - E(kx)+uj(qu))d(k 2 - h - q + G), 
/ BF BP BW BT* 1 \ 

+ ("-2D 

/ <9 F BW \ 

BiihM) = (^^)qJ- " ^ W (1L22) 
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and 

f(E(k)) = f k = [exp /?(£(*;) - E F ) + 
N{uj(qij)) = N(qu) = [exp Pu(qv) - 

The correlation functions (P\\P{) vanish in the weak scattering limit. The generalized electron 
numbers in Eq. (jll.l3j) become 

ty = -Tr(pi\(i\);i\) =mPY i $\%)- F f k (l- f k ). (11.23) 
11.2 Temperature dependence of R 

Let us consider the low-temperature dependence of the electrical resistivity obtained above. For 
this region we have 

lim (3f k (l-f k )^S(E(k)-E F ). 

T->0 

Thus the generalized electron numbers Ni in Eq. (|11.23p do not depend on temperature, and the 
temperature dependence of R in Eq. ljll.lip is given by the correlation functions (|11.19|) and 
(jll.2ip . For the term arising from the electron-electron scattering we find 



(P? e ;Pt)=P ' 



/ / / dE(k 1 )dE(k 2 )dE(k 3 )F} l (E(ki),E(k 2 ),E(k 3 ))- 

" fk 2 )fk 3 [l ~ f(E(h) - E(k 2 ) + E(k 3 ))\, (11.24) 



where 



3 ' 



F^EihlEik^Eik,)) = ^^^E / ^ J ^ J ^ 

G 

A j (k 1 ,k 2 ,k 3 , fci - fc 2 + fc 3 + G)Ai(k 1 ,k 2 , h, k 1 -k 2 + k 3 + G) 

I dE 1 1 8E 1 1 dE I 
I dk\ 1 1 dk2 1 1 <9A,'3 I 

5(E(ki) - E(k 2 ) + E{k 3 ) - E{k x -k 2 + k 3 + G)). (11.25) 

With the substitution 

x = P{E{k x ) - E F ), y = f}(E(k 2 ) — E F ), z = (3{E{k 3 ) — E F ), 
the expression (lll.24j) reads 

<"pE max —Ep 111 



(Jf ;i? e )=/r 2 / 

J —I 



PE F J J 1 + exp(x) 1 + exp(-y) 1 + exp(z) 

dxdydz „i ,x _ y „ z , 

F 1 (- + J E F ,^ + S F ,- + J E F ) 



r 2 ^f- (ii.26) 



1 + exp(-x + y - z) v /3 (3 p 

It is reasonably to conclude from this expression, that in the limits 

X\vn.fiE F -> oo, lim/3(£ ma x - £ F ) oo , 

the electron-electron correlation function for low temperatures becomes proportional to T 2 for 
any polynomial § 3 {k). 

For the electron-phonon contributions to the resistivity the temperature dependence is given by 
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the Bose distribution function of phonons N(qv). Because of the quasi-momentum conservation 
law q = — k\ — G the contribution of the electron-phonon Umklapp processes freezes out at 
low temperatures as [exp f3uj(q m i n ) — 1] , where q m in is the minimal distance between the closed 
Fermi surfaces in the extended zone scheme. For electron-phonon normal processes as well as for 
electron-phonon Umklapp processes in metals with an open Fermi surface the quasi-momentum 
conservation law can be fulfilled for phonons with q — > being excited at low temperatures solely. 
To proceed further we use the relation oo(qv) = v^q/q) and the periodicity of the quasiparticle 
dispersion relation with the reciprocal lattice vector G. Taking into account the relation + 
G) = & l (k), we find to the first non- vanishing order in q 

— * Q 

9k,k+q+G ~ 9(-g^)9k,k'\k'=k, (11.27) 

— * r\ 

Bj (k, k + q + G)^ q{l—)Bj{k, h>)\ w=k , (11.28) 

5(E(k + q + G)- E(k) + u(qv)) « -5{-^r + v v Q ). (11.29) 

q q OK 

Hence, we have 

(pei.peij^ (11.30) 

q2 rqmax r r 

m2(3 J^f E I J M0 q )d9 q J ^Jexp(^g) - l]- l Fl (6 q , <p q ), 



where 



(cf ^{k, k% l=k )f k (i - /*wf H + (ii.3i) 

In the k- integral in Eq. (|ll,3U|) the integration limits have to be chosen differently for normal and 
Umklapp processes. With the substitution x = ftv^q we find 

(P?;Ph = (n.32) 

777Q 2 -^-^ 1 rPv^qmax r 5 r r 

^J^f E JffiB I dx —i I J d<P q F^9 q , Vq ) = Af t T . 

Thus we can conclude that the electron-phonon correlation function is proportional to T 5 for any 
polynomial & l k . It is worthy to note that for an open Fermi surface this proportionality follows for 
normal and Umklapp processes either. For a closed Fermi surface the Umklapp processes freeze 
out at sufficiently low temperatures, and the electron-phonon normal processes contribute to the 
electrical resistivity only. With the aid of the Eqs. (jll.26p and fjll .32 j) the generalized kinetic 
equations (jll.7p becomes 

n 

Vi [a^T 2 + A^T^ = eENj. (11.33) 

i=l 

For simplicity we restrict our consideration to two parameters x>\ and vi describing the homo- 
geneous shift and one type of deformation of the Fermi body. Taking into consideration more 
parameters is straightforward but does not modify very much qualitative results. Finally, the 
expression for the electrical resistivity b ecomed^SSEnS] 

R _ g (AftT 2 + Afir 5 )(Af 2 r 2 + Af 2 T*) - {a\%t 2 + ^f 2 r 5 ) 2 

3e 2 N 2 (A e 2 %T 2 + Af 2 T 5 ) + N^A^T 2 + A^T 5 ) - 2N 1 N 2 {A\IT 2 + Af 2 T^>)' K ' ' 
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In general, a simple dependence R ~ T 2 or P ~ T 5 can be expected only if one of the scattering 
mechanisms dominate. For example, when A?J ~ 0, we find 

p_ " (Af^f 2 -ylf 2 )T 5 

3e 2 iV^ll + iV|^4fi - 2iViiV 2 ^fi ' ' 

If, on the other hand, the deformation of the Fermi body is negligible (t>2 = 0), then from 
Eq. (|11.34jl it follows that 

R=^^(AttT* + A<*T 5 ). (11.36) 

It is interesting to note that a somewhat similar in structure to expression (|11.34p have been 
used to describe the resistivity of the so-called strong-scattering metals. In order to improve our 
formula for the resistivity derived above the few bands and the interband scattering (e.g. s — d 
scattering) as well as the phonon drag effects should be taken into account. 

11.3 Equivalence of NSO approach and Kubo formalism 

Equivalence of the generalized kinetic equations to the Kubo formula for the electrical resistivity 
can be outlined as following. Let us consider the generalized kinetic equations 

^P n (-iTr(p[P n ,P m ]) + (P n ;P m >) = ^ (Tr(pP e (-iX)P m ) + (P e ;P m >). (11.37) 

n 

To establish the correspondence of these equations with the Kubo expression for the electrical 
resistivity, it is necessary to express the operators of the total electron momentum P e and the 
current density j in terms of the operators P m . In the other words, we suppose that there exists 
a suitable set of coefficients with the properties 

P e =J2®i P i> 3 = ^Pe- (11.38) 

i 

We get, by integrating Eq. (|11.37p by parts, the following relation 

TfJ f dXTr(pP n (-iX)P m ) - s{P n ;P m )) = — (P e ;P m ). (11.39) 

Supposing the correlation function (P n ; P m ) to be finite and using Eq. (jll.38p we find in the limit 

e ~* ° p 

EE^"/ d\Tr{pP n (-i\)P i ) = —(P e -P e ). (11.40) 

From the equality (jll.40p it follows that 

V F n f d\Tr(pP n (-i\)P e ) = — (P e ; P e >. (11.41) 
n Jo 

Let us emphasize again that the condition 

lim e(P n ;P m ) =0 (11.42) 

e^O 

is an additional one for a suitable choice of the operators P m . It is, in the essence, a certain 
boundary condition for the kinetic equations (|11.37p . Since the Kubo expression for the electrical 
conductivity 

a- (P e ;P e ) ~ (P n ;P m > (11.43) 
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should be a finite quantity, the condition (jll.42p seems quite reasonable. To make a following 
step, we must take into account Eq. (111.10p and (lll.lip . It is easy to see that the r.h.s. of the 
Eq. (jll.4ip is proportional to the current density. Finally, we reproduce the Kubo expression 
(110.18P (for cubic systems) 

" = 1=3^<^>' < 1L44) 
where the proportionality of the F m to the external electrical field has been taken into account. 
It will be instructive to consider a concrete problem to clarify some points discussed above. There 
are some cases when the calculation of electrical resistivity is more convenient to be performed 
within the approach of the generalized transport equations (|1Q.27|) and ()10.28p than within the 
Kubo formalism for the conductivity. We can use these two approaches as two complementary 
calculating scheme, depending on its convenience to treat the problem considered P^l To clarify 
this, let us start from the condition 

Tr(p LR B i )=Tr(p a B i ), (11.45) 

where the density matrices plr and p s were given by equations fllQ.23j) and (110.25P - For operators 
Bi which can be represented by linear combinations of the relevant observables Ai the equation 
(I11.45P is fulfilled exactly where for other operators the equations (I11.45P seems to be plausible 
if the relevant observables have been chosen properly. The conditions (jl 1.45|) make it possible to 
determine a set of parameters which can be used in approximate expressions for the correlation 
functions. In simple cases the conditions (I11.45P even allow to calculate the correlation functions in 
(|10.27p without resorting to another technique. As an example we consider the one-band Hubbard 
model (14"T4"1) 

H = H e + H ee (11.46) 

in the strongly-correlated limrtP^ES \t\/U <C 1. It well known that in this limit the band splits into 
two sub-bands separated by the correlation energy U. In order to take into account the band-split 
we have to project the one-electron operators onto the sub-bands. The relevant operators will 
have the form 

ka 

where 

<i = Nl Yl em ^ ] 4^ta^_ a (11.48) 
ij 

with the projection operators 

nt„ = h-° ^ a = + ' (11.49) 

[(1 - ni-o) if a = + . 

Here P aa is the operator of the total momentum of the electrons in the sub-band a, and P a @ , (a 7^ 
/3) describes kinematical transitions between the sub-bands. It can be shown that correlation func- 
tions (paP-pib) and (p* a -pPP) vanish for a ^ /3 in the limit \t\/U < 1. To make an estimation 
it is necessary to decouple the higher correlation functions in \t\/U and take into account near- 
est neighbor hopping terms only. Then the correlation functions (p QC *;P Qa ) ( a = i) can De 
calculated directly by means of Eq. (|11.45p . where Bi —> p aa . The conductivity becomes 
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where 

k 

and 

(W) = Tr(pW), P = ±exp(-PH). 

Here z is the number of nearest neighbours and t the nearest neighbour hopping matrix ele- 
ment. With the well-known expressions for the mean values {n c ^n a _ a ) we find the conductivity in 
dependence on the electron number in the formPSSl 

a = WT- 1 - , ' , (0<n<l), (11.52) 

a = WT- l -=exp(-U/2kT), (n = l), (11.53) 
v2 

a = WT- W 1 - n/2) ^- n \ (2>n>l). (11.54) 



n 



Thus, it was shown in this and in the preceding sections that the formalism of the generalized 
kinetic equations has certain convenient features and its own specific in comparison with Kubo 
formalism. The derived expressions are compact and easy to handle with. 



11.4 High-temperature resistivity and MTBA 

At high temperatures the temperature dependence of the electrical resistivity R of some transition 
metals and highly resistive metallic systems such as A15 compounds may deviate substantially 
from the linear dependence, which follows from the Bloch-Gruneisen law. These strong deviations 
from the expected behavior with a tendency to flatten to a constant resistivity value was termed 
by resistivity saturatiotW^ and have been studied both experimentally and theoretically by many 
authors ! 315 " 32 ' The phenomenon of resistivity saturation describes a less-than- linear rise in dc 
electrical resistivity R when temperature T increases. It was found that this effect is common in 
transition metal compounds (with pronounced d-band structure) when R exceeds ~ 80 pflcm, and 
that R seems bounded above by a value R m ax ~ 150 fiilcm which varies somewhat with material. 
In Ref.j 3 -^ in particular, it was formulated that the electrical resistivity, R, of a metal is usually 
interpreted in terms of the mean free path (the average distance, I, an electron travels before it 
is scattered). As the temperature is raised, the resistivity increases and the apparent mean free 
path is correspondingly reduced. In this semi-classical picture, the mean free path cannot be much 
shorter than the distance, a, between two atoms. This has been confirmed for many systems and 
was considered to be a universal behavior. Recently, some apparent exceptions were found, includ- 
ing alkali-doped fullerenes and high-temperature superconductors P^Sl However, there remains the 
possibility that these systems are in exotic states, with only a small fraction of the conduction elec- 
trons contributing to the conductivity; the mean free path would then have to be correspondingly 
larger to explain the observed resistivity. The authors of RefF^Sl performed a model calculation of 
electron conduction in alkali-doped fullerenes, in which the electrons are scattered by intramolec- 
ular vibrations. The resistivity at large temperatures implies I ~ a, demonstrating that there is 
no fundamental principle requiring I > a. At high temperatures, the semi-classical picture breaks 
down, and the electrons cannot be described as quasiparticles. Recent review of theoretical and 
experimental investigations in this field was given in RefsP2BH230] ^£ or clis CUSS ion of the electronic 
thermal conductivity at high temperatures see RefP^j). 

The nature of saturation phenomenon of electrical resistivity is not fully understood. Resistivity 
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of a metallic system as a function of temperature reflects an overall electron-phonon interaction 
effects as well as certain contribution effects of disorder J222H3MI There have been some attempts 
to explain the saturation phenomenon in the framework of of the Boltzmann transport theory 
using special assumptions concerning the band structure, etc. The influence of electron-phonon 
scattering on electrical resistivity at high temperatures was investigated in R,ef s E2H1I3231 i n the 
framework of the Frohlich Hamiltonian for the electron-phonon interaction. In Ref. 323 authors 
calculated a temperature-dependent self-energy to the lowest non- vanishing order of the electron- 
phonon interaction. 

However, as it was shown above, for transition metals and their disordered alloys the modified 
tight-binding approximation is more adequate. Moreover the anisotropic effects are described bet- 
ter within MTBA. Here we consider a single-band model of transition metal with the Hamiltonian 

H = H e + Hi + H ei . (11.55) 

The electron subsystem is described by the Hubbard model (I4.14p in the Hartree-Fock approxi- 
mation 

H e = E(ka)al a a ka , E(ka) = E{k) + ^ ^(n fc _ CT ). (11.56) 

ka p 

For the tight-binding electrons in crystals we use E(k) = 2 J2 a ^ C^«) c os(kR K ), where t°(R K ) is 
the hopping integral between nearest neighbours, and R k (k = x, y.z) denotes the lattice vectors in 
a simple lattice in an inversion centre. For the electron-phonon interaction we use the Hamiltonian 
(14661) 

H ei = Y,12 V $^+q)Qqat +q aak*, Qq= J— (b q + bl q ), (11.57) 

a kg V 2u} {q) 

where 

V(k, k + q)= {n 2 )1/2 Y^t°(R K )^^[ S mRj - smR K (k + q)]. (11.58) 

The one-electron hopping t°(R K ) is the overlap integral between a given site R m and one of the 
two nearby sites lying on the lattice axis R K . Operators b\ and b q are creation and annihilation 
phonon operators and ui(q) is the acoustical phonon frequency. N is the number of unit cells in 
the crystal and M is the ion mass. The e u (q) are the polarization vectors of the phonon modes. 
For the ion subsystem we have 

H i = J2^Q)(4 b i + 1 / 2 )- ( 1L59 ) 

For the resistivity calculation we use the following formula^ 3 -' 

R = ro ^ F,F ^ ^ ^ . (11.60) 

3eW 2 l + (l/3mM)(P;F) 

Here M is the effective number of electrons in the band considered 

1 [V 

N = / dXTr(pP(-ihX)P) (11.61) 

3m J 

and P is the total momentum operator 

m sr ^, I dEika) \ + . . 

P = -r >J ~ — n ka, n ka = al a a ka . (11.62) 

n \ dk ) 
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The total force F acting on the electrons is giving by 

F = -AH, P) = — r- ^ V(k, k + q){v k+qi a - v k>a )Q q al +qrT a ka , (11.63) 

with the velocity defined as v kcr = dE(ka)/hdk. It is convenient to introduce a notation 

V(k, k + q) _ iAF q 

Correlation functions in Eq. (jll.60|) can be expressed in terms of the double-time thermodynamic 
Green functions 

-, -> 27T? -> -. 

(F;F) = —((F\P)). m , (11.64) 



h 

(P;F) = ™«P|F»_, te , (11.65) 



«pV»-^ = 

l f°° 

— / dte et #(-t)Tr(pL4,PW]), (11.66) 

where j4 represents either the momentum operator P or the position operator R with P = 
im[H,R]/h. 

We find the following relation 

(F;F) = 

lixm V(k, k + q) _ + ,+ 

2^ — ~ n,a){{ak +qa ak*V>q + V_ q )\P))- m . (11.67) 

Thus we obtain 

~ { X, — k~T\ ^ k+q ' a ~ ^' <J - ) (^ a fc+ga a ^ 6 </l - (11.68) 

kqa 

(( a k+qa a k<T b -q)\A))-ihe) ■ 

Calculation of the higher-order Green functions gives 

(E(k + qa) - E(ka) - % - ihs)((al +qa a ka B q \A)}„ m = 

Tkq[A) + ?(^Sr ((a ^ afe - 9 ' CT(v ~ - 

q' 

V(k + q,k + q + q') + , + \ 



Ey(F, k' — q) 4 + 
/^-t^t ( K+^afc-rOfc'-^afc'o- 1 A) ) ( 1 1 .69) 
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with notation 



We find also 



tt q = huj(q) -> B q = b q ; E q = -fuj(q) ->- B q = b ] 



im 



T kq(P) = ^~{Vk+q,a ~ V k ,a) {a\ +qa a ka B q ) , 



2tt 

T kq (R) = -— Q^(a k+q+qfa ak+q'<r B <i) 6 q',o, 
M= ^J2 5 q'flVk,a-^{al +q , a a k+ql(T ). 

kq'cr 



{{ a l +q a a k-q'ob q b ql \A))_ ih£ = (11.70) 

t V(k,k-q') ((al +qa a ka b q \A)). m 



V/0- 



y^tyy (E(k + qa) - E(k - q'a) - e q > - tutf) - ihe) ' 



({al +q+q/a a ka b q b q/ \A))_ ihe = (U-71) 

_ , t , ) V{k + g + q',k + q) ({a k+qa akgb q \ A))_ ite 

( q ' q ' )0 {E(k + q + q'a) - E(ka) - E q , - hoj(q') - ihe) ' 

Here the symmetry relations 

V(k - q', k) = V*(k, k - q>); V(k + q + q', k + q) = V*(k + q,k + q + q') 

were taken into account. 

Now the Green function of interest can be determined by introducing the self-energ; 



((ai +qa ak-q>*b q B q \A))^ ih£ = 

_/u\ u \ T kq (A) 

U' i'° ( E (k + qa) - E{k - q'a) - E q - M kqa (E q , -ihe) - ihe) ' { > 

The self-energy is given by 

M kqa (E q ,-ihe) = 

V — (V b ,)\ W$j-?)\ 2 + 

2co(q') V 1' q /l E{k + qa) - E{k - q'a) - E q > - hu(q') - ihe 



\V(k + q,k + q + q- 



hi2 



E(k + q + q'a) - E(ka) - E q < - hco(q') - ihe 



I (b b )\ \ v \ n i n y j\ 1 

[ q ' q ' )[ E{k + qa) - E(k - q'a) - E q > + huj(q') - ihe 



,tu \V(k,k-q')\ 2 

\V{k + q,k + q + q')\ 2 . 

E(k + q + q'a) - E(ka) - E q > + frw(q') - ihe 



(11.73) 



In equation ()11.73|) the energy difference (E(k + qa) — E(k — qa) — E q )) is that for the scattering 
process of electrons on phonons, while emission or absorption of one phonon is possible, corre- 
sponding to E q . These scattering processes are contained in the usual Boltzmann transport theory 



72 



leading to the Bloch-Gruneisen law. The self-energy M kqcr describes multiple scattering correc- 
tions to the Bloch-Gruneisen behavior to second order in V, which depends on the temperature 
via the phonon occupation numbers. 

Furthermore, it is assumed that the averages of occupation numbers for phonons in the self-energy 
and for electrons in the effective particle number, are replaced by the Bose and Fermi distribution 
functions, respectively 

(b\b q ) = N q , N q = [exp(phtj(q)) - l] -1 , (11.74) 
(a{ a a ka ) = f k , f k = [exp(/3E(ka) - E F ) + l]" 1 . (11.75) 

This corresponds to neglecting the influence of multiple scattering corrections on the phonon and 
electron distribution functions. 

In order to calculate the expectation values in the inhomogeneities, Eqs. (|11.72p and (jll.73|) . the 
spectral theoremP^E^ should be used. In the lowest non-vanishing order of the electron-phonon 
interaction parameter V we obtain 

(al +qa a ka B q ) = 

V(k + g,k) f [exp(/?(£(fc + go) - Ejka) - E q )] - 1 
-jm~ h+qi " h) q{ q) E(k + q *)-E(kc)-E q ' (1L76) 

with 

-• 1-expi-pE,) \N, if B, = (.*_,; E q — -Hw(q). 

Applying the approximation scheme discussed above we have found the following expressions for 
the Green function with A = P 

{{a\ +qa a k(T B q \P))_ ih£ = 

im V(k + q,k) _ P k q a( E q ) . , 

o rm~\ w^-^,«r h — TW\ iT7 ^7' 11.78 

2vr y/2uj(q) \l kqa {E q ) - M kqa - ihe 

and for the Green function with A = R 



((al +qa a ka B q \R)). ihe 
hV(k + q,k).^ _ p l q(T { E q) ~ p l q M 



.-> {v k+qa -v ka )- — — — — . (11.79) 

2vr y/2uj(q) ^kqa{E q ) - M kqa - the 

We have introduced in the above equations the following notation: 

P$A E l) = fk+qO- ~ fk>q(Eqhl(n kq(T (E q )), (11.80) 
P^(Eq) = (V k+ q„ ~ V ka )f k+q (l - f k )Vq(E q )( l2 (n kqa (E q )) - ^^^^) > (U-81) 

P k q l( E i) = W 1 ~ fk)vq(E q )f3 7l (Q kq(7 (E q ))lf k v k(7 - (1 - f k+q )v k+q<T ], (11.82) 

with 

fn PeMfMkqa(Eq))-l 

ln(^(E q) ) = { ^ {Eq))n , (11.83) 
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and 

Sl kqa (E q ) = E{k + qa) - E{ka) - E q . (11.84) 
For the effective particle number we find 

N=^mf3 y £(v k ) 2 f k (l-f k ). (11.85) 

k 

Before starting of calculation of the resistivity it is instructive to split the self-energy into real 
and imaginary part (e — > 0) 

lim M kqa {huj{q) ± ihe) = ReM kqa (hio(q)) T HmM kq(T (frw(q)) (11.86) 

and perform an interchange of variables k + q — )• k;q — > —q. Now for the relevant correlation 
functions we obtain the expressions 

= E f f * - t^ 2 P<£(Mg)) WMg)), (n.87) 

(P;F)= (11.88) 

where 

l m -^fcgo-(^ ii; (9)) 

5^(^(9)) = _ ReM fc?(T (^;(g)))2 + (ImM fe9(T (M?))) 2 ' (U ' 89) 

In order to obtain Eq. ()11.89p the following symmetry relation for the self-energy was used 

M kqa {hw{q) - ihe) = -M kqa {huj{q) + ihe). (11.90) 

The inspection of both the correlation functions (F; F) and (P; F) shows that its include two 
dominant parts. The first one is the scattering part S kqa (huj(q)), which contains all the information 
about the scattering processes. The second part describes the occupation possibilities before 
and after the scattering processes {P^ a {hw{q)),P^ a {hLo{q)),P^ a {huj{q))) : and includes both the 
Fermi and Bose distribution functions. The approximation procedure described above neglects 
the multiple scattering corrections in these factors for the occupation possibilities. 
For further estimation of the correlation functions the quasi-elastic approximation can be used. 
In this case, in the energy difference Q kqr7 (hu(q)), the phonon energy huj(q) can be neglected 
n kq(7 {hw{q)) ~ O, kqa (0). The phonon wave number q only is taken into account via the electron 
dispersion relation. Furthermore, for the Bose distribution function it was assumed that 

(blb q ) = (b q b\) ~ (^(g))" 1 . (11.91) 

This approximation is reasonable at temperatures which are high in comparison to the Debye 
temperature @d- 

It is well knowrPI^H f r0 m Bloch-Gruneisen theory that the quasi-elastic approximation does not 
disturb the temperature dependence of the electrical resistivity at high or low temperatures. The 
absolute value of the resistivity is changed, but the qualitative picture of the power law of the 
temperature dependence of the resistivity is not influenced. 



74 



In the framework of the quasi-elastic approximation the scattering contribution can be represented 
in the form 

Skqa = m - ^mj 2 + QnMz,.)) 2 ' (1L92) 

Here the real and imaginary parts of the self-energy have the following form 

r> ii#e _kT^ ( \V(k,k+p)\* 1 

kqa h uj( P ) 2 K E{k + qa) - E{k - pa)' 

\V(k + q,k + q+p)\ 2 1 N 

w(p) 2 \E(Jfe + g+|X7) -E{ka)')' K ' 

ImM,V = ^E( |F( !; ( ^ Pl|2 ^( fc + «") " E(k-pa)) + 

(n) 

The occupation possibilities given by P^ qu can be represented in the quasi-elastic approximation 
asPSl 



P ^^f^qj 6{L: ' - t:{krr)]: (IL<)r,) 



= = 0- ( 1L96 ) 

In this approximation the momentum-force correlation function disappears (P; F) ~ 0. Thus we 
have 

R -3^ { ^' (1L97) 

,-, 2m 2 kT ^\V{k,k + q)\ 2 , NX 

(F; F) = —f— £ J,( )2 " v **?&{Ef ~ E{ka)) ■ 

kqa 

ImMf 

(11.98) 



The explicit expression for the electrical resistivity was calculated in Ref. 3 — The additional 
simplifying assumptions have been made to achieve it. For the electrons and phonons the following 
simple dispersion relations were taken 

^) = !4; u(q) = v \q\, V(k,k + q)~ v^f. 
2m* 

(l)e 

It was shown that the estimation of P^J a is given by 

Then, for the electrical resistivity the following result was found 

„ ZV^qp T ["■> f 1 4 



-1 



PV™*)(^F + q/2) - ReM^J 2 + [ImM^ 2 )]2 " 



(11.100) 
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This result shows that the usual Bloch-Gruneisen theory of the electrical resistivity can be cor- 
rected by including the self-energy in the final expression for the resistivity. The Bloch-Gruneisen 
theory can be reproduced in the weak scattering limit using the relation 

lim kFqz 



Re(lm)A/-,o [(H 2 q/m*)(zk F + q/2) - ReM fe e F(? J 2 + [ImM^J] 2 

nm* q . 

= ——5(z + ---). (11.101) 

Inserting Eq. (|11.10ip in the resistivity expression Eq. (|ll,100p gives the electrical resistivity R w in 
the weak scattering limit, showing a linear temperature dependence 

3irV 2 (m*) 2 q 5 n T 

^t^M^w (1L102) 

In this form, the resistivity formula contains two main parameters that influence substantially. 
The first is the Debye temperature characterizing the phonon system, and, the second, the 
parameter a = (Vm*) 2 describing the influence of both the electron system and the strength of 
the electron-phonon interaction. The numerical estimations^^ were carried out for Nb and gave 
the magnitude of the saturation resistivity 207 ^0,cm. 

In the theory described above the deviation from linearity in the high-temperature region of the 
resistivity may be caused by multiple scattering corrections. The multiple scattering processes 
which describe the scattering processes of electrons on the phonon system by emission or absorp- 
tion of more than one phonon in terms of self-energy corrections become more and more important 
with increasing temperature. As was shown above, even for simple dispersion relation of electrons 
and phonons within one-band model the thermally induced saturation phenomenon occurs. For 
the anisotropic model within MTBA the extensive numerical calculations are necessary. 
In subsequent paper P^H Christoph and Schiller have considered the problem of the microscopic 
foundation of the empirical formula-^! (parallel resistor model) 

111 , 

+ - j=- (11.103) 



R(T) RSBT(T) Rmax(T) 

within the framework of the transport theory of Christoph and KuzemskyP^SI The parallel resis- 
tor formula describing the saturation phenomenon of electrical resistivity in systems with strong 
electron-phonon interaction was derived. In Eq. (|11.103p Rsbt(T) is the resistivity given by the 
semiclassical Boltzmann transport theory Rsbt(T) ~ T and the saturation resistivity R m ax cor- 
responds to the maximum metallic resistivity The higher-order terms in the electron-phonon 
interaction were described by a self-energy which was determined self-consistently. They found 
for the saturation resistivity the formula-^' 

_ 3^ 3 h q% 1 
Umax ~ 32 e 2 4 \P(q D /2k F )\ ' 1 U > 

Within the frame of this approach, the saturation behavior of the electrical resistivity was ex- 
plained by the influence of multiple scattering processes described by a temperature-dependent 
damping term of one-electron energies. In the standard picture, the conventional linear temper- 
ature dependence of the resistivity R ~ T(T 3> &£>) is explained by taking into account that 
the number of phonons is proportional to the temperature and, moreover, assuming that the 
electron momentum is dissipated in single-phonon scattering processes only. For an increasing 
number of phonons, however, the multiple scattering processes become more important and the 
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single scattering event becomes less effective. This argument coincides in some sense with the 
Yoffe-Regel criteriorP^HSU stating that an increase in the number of scatterers does not result 
in a corresponding increase of the resistivity if the the mean free path of the electrons becomes 
comparable with the lattice distance. Indeed, the saturation resistivity (jll. 104H coincides roughly 
with the inverse minimal metallic conductivity which can be derived using this criterion. 

12 Resistivity of Disordered Alloys 

In the present section a theory of electroconductivity in disordered transition metal alloys with 
the proper microscopic treatment of the non-local electron-phonon interaction is considered. It 
was established long ago that any deviation from perfect periodicity will lead to a resistivity con- 
tribution which will depend upon the spatial extent and lifetime of the disturbance measured in 
relation to the conduction electron mean free path and relaxation time. It is especially important 
to develop a theory for the resistivity of concentrated alloys, because of its practical significance. 
The electrical resistivity of disordered metal alloys and its temperature coefficient is of a consid- 
erable practical and theoretical i n t e rest.^^'IS3{lZ2]|233|325l[336H3Ml The work in this field has been 
considerably stimulated by Mooij paper^^ where it has been shown that the temperature coeffi- 
cient of the resistivity of disordered alloys becomes negative if their residual resistivity exceeds a 
given critical value. To explain this phenomenon, one has to go beyond the weak-scattering limit 
and to take into account the interference effects between the static disorder scattering and the 
electron-phonon scattering I23SJESH 

In the weak-scattering limfP^ the contributions of impurity and phonon scattering add to the 
total resistivity without any interference terms (Matthiessen rule). For disordered systems many 
physical properties can be related to the configuration-averaged Green functions.'^' There were 
formulated a few methods for calculation of these averaged Green functions. It vas found that 
the single-site CPA (coherent potential approximation)P32H352] p rov id es a convenient and accurate 
approximation for it ES&EES The CPA is a self-consistent met ho cP^H26Qj th a t predicts alloy elec- 
tronic properties, interpolating between those of the pure constituents over the entire range of 
concentrations and scattering strengths. The self-consistency condition is introduced by requiring 
that the coherent potential, when placed at each lattice site of the ordered lattice, reproduces 
all the average properties of the actual crystal. The coherent-potential approximation has been 
developed within the framework of the multiple-scattering description of disordered systems!*^ 
A given scatterer in the alloy can be viewed as being embedded in an effective medium with a 
complex energy-dependent potential whose choice is open and can be made self-consistently such 
that the average forward scattering from the real scatterer is the same as free propagation in the 
effective medium. The strong scattering has been first considered by Velicky^H in the framework 
of the single-site CPA using the Kubo-Greenwood formula. These results have been extended 
later to a more general models F^^^ 

The first attempt to include the electron-phonon scattering in the CPA calculations of the re- 
sistivity was given by Chen et aZP^ A model was introduced in which phonons were treated 
phenomenologically while electrons were described in CPA. The electron-phonon interaction was 
described by a local operator. Chen, Weisz and Sher^HI (CWS) have performed a model calculation 
on the temperature dependence of the electronic density of states and the electrical conductiv- 
ity of disordered binary alloys, based on CPA solutions by introducing thermal disorder in the 
single-band model. They found that the effect of thermal disorder is to broaden and smear the 
static-alloy density of states. The electrical conductivity in weak scattering alloys always decreases 
with temperature. However, in the strong-scattering case, the temperature coefficient of conduc- 
tivity can be negative, zero, or positive, depending on the location of the Fermi energy. Brouers 
and Brauwers^MJ have extended the calculation to an s — d two-band model that accounts for 
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the general behavior of the temperature dependence of the electrical resistivity in concentrated 
transition metal alloys. In RefP^ a generalization of CWS theory^Sl was made by including 
the effect of uniaxial strain on the temperature variation of the electronic density of states and 
the electrical conductivity of disordered concentrated binary alloys. The validity of the adiabatic 
approximation in strong-scattering alloys was analyzed by Chen, Sher and WeiszF^l n was shown 
that the electron screening process in the moving lattice may be modified by lattice motion in 
disordered alloys. Were this modification significant, not only the effective Hamiltonian but also 
the whole adiabatic approximation would need to be reconsidered. 

A consistent theory of the electroconductivity in disordered transition metal alloys with the 
proper microscopic treatment of the electron-phonon interaction was carried out by Christoph 
and KuzemskyF^ They used the approach of paper where a self-consistent microscopic theory 
for the calculation of one-particle Green functions for the electron-phonon problem in disordered 
transition metal alloys was developed. However, this approach cannot be simply generalized to 
the calculation of two-particle Green functions needed for the calculation of the conductivity by 
the Kubo formula. Therefore, for the sake of simplicity, in their study Christoph and Kuzemsky^^ 
neglected the influence of disorder on the phonons. Thus, in the model investigated here, in con- 
trast to the GWS approach j^l t ne dynamics of the phonons is taken into account microscopically, 
but they are treated as in a virtual reference crystal. 

For a given configuration of atoms the total Hamiltonian of the electron-ion system in the substi- 
tutionally disordered alloy can be written in the f orrD [!2QJ|3io] 

H = H e + Hi + H ei , (12.1) 

where 

H e = ti a \a a ia + E Uj4a a jT ( 12 - 2 ) 
icr ija 

is the one-particle Hamiltonian of the electrons. For our main interest is the description of the 

electron-phonon interaction, we can suppose that the electron-electron correlation in the Hubbard 

form has been taken here in the Hartree-Fock approximation in analogy with Eq. (lll.56p . 

For simplicity, in this paper the vibrating ion system will be described by the usual phonon 

Hamiltonian 

H, t =J2"{&)(bU«» + V2). (12.3) 

qv 

The electron-phonon interaction term is taken in the following formpQl 

H « = E E T ^ u ? - u >l a ^ ( 12 - 4 ) 

ij OLG 

where uf(a = x,y,z) is the ion displacement from the equilibrium position R{. 
In terms of phonon operators this expression can be rewritten in the form 

Hei = E E + ^kW, (12.5) 

i^zj qua 

where 

A (ij) = 1° t o R J- R : e u (q)(e^ - e^Y (12.6) 

q KJJ ^2{M)Nu{qv) l3 \Rj-Ri\ ^ ' 

Here u{qv) are the acoustic phonon frequencies, (M) is the average ion mass, e u (q) are the polar- 
ization vectors of the phonons, and go is the Slater coefficient originated in the exponential radial 
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decrease of the tight-binding electron wave function. It is convenient to rewrite this expression in 
the form 

H « = Y,Y. A Sj){\ + bl q )a\a j , (12.7) 
where the spin and phonon polarization indices are omitted for brevity. 

The electrical conductivity will be calculated starting with the Kubo expression for the dc con- 
ductivity: 

a(ie) = -((J\P)) i£ , (e->0+), (12.8) 

where P = e^^i?joJoj and R4 is the position vector; m/eJ = m/eP is the current operator of 
the electrons. It has the form 

J = -ie ^(Ri - Rj)tija\aj. (12.9) 
u 

Then the normalized conductivity becomes 

s 2 
Q 



ij I 

where Q, is the volume of the system. It should be emphasized here that a 01 " depends on the 
configuration of the alloy. A realistic treatment of disordered alloys must involve a formalism to 
deal with one-electron Hamiltonian that include both diagonal and off-diagonal randomness PiMSBD] 
In the present study, for the sake of simplicity, we restrict ourselves to a diagonal disorder. Hence 
we can rewrite hopping integral ty as 

Uj = ^Y1 exp\ik(Ri - Rj)] . (12.11) 

k 

Thus to proceed it is necessary to find the Green function Gy^m = {{a\aj\a\a m )) . It can be calcu- 
lated by the equation of motion method. Using the Hamiltonian ()12.ip we find by a differentiation 
with respect to the left hand side 

Hij^nG n r t i m (ui) = {a\a m )5ij — {a\aj)8 m i + 

nr 

- n)e^ {{a\a n {b q + bl q )\a\a m )) - 

qn 

A q (n - i)e iA {(ai aj (b q + bl q )\a}a m ))) , (12.12) 

where 

Hij,m — (^ ~i~ ^rj^ni^rj ^jr^ni ~t~ ^ni^rj- (12.13) 

We define now the zeroth-order Green functions G®- lm that obey the following equations of motion 

^2 H ij,mG^ rM = {a\a m )8ij - (aJaj)S mi , (12.14) 
nr 

^ H rn ,lm G{ ij,nr = { a \ a m)^lj ~ (al a j)$mi, (12.15) 
nr 
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where Eq. (|12.15p has been obtained by a differentiation with respect to the right hand side of 
G®j lm . Using these definitions it can be shown that 

^((aja n )(5 rt - (ala t )5 sn )G nrt i m (uj) = 

nr 

5E((afa m )% - (o}o J -)5 rn i)G ! S tiii (w) + 
u 

EE(^' " n)e^((a\a n {b q + bl q )\a\a m )) - 



ijn q 

Aq{n _ i) e ^((aia 3 (b q + bl q )\a\a m )))G stjji (u;) . (12.16) 

The r.h.s. higher order Green functions can be calculated in a similar way. To proceed we 
approximate the electron-phonon Green function as 

{{aia r b\b q \B)) ~ N(q)((ala r \B)). (12.17) 

Here N(q) denotes the Bose distribution function of the phonons. 
As a result we find 

^2{(ala n )5 rt - {ala t )8sn){{an a rbq\a\a m )) = 



w (?) 5E ( ( a J a i h i I a I a ™) ) G °st , j» ( w ) 

£(1 + JV(g))(A_,0' " n)e-^G in , /m M 
A_ g (n - i)e-^G nhlm ^))G° st ^) 



Y J Y, A -^ n -P>^ Rn (M)Gn j ,i m ~ (ala 3 )G ipM )G° stJi (u;) (12.18) 



ij np 



and a similar equation for {{ana r b- q \a\a m )) . 

In the above equations the Green functions G and G° as well as the mean values (ajou) which 
can be expressed by one-particle Green functions depend on the atomic configuration. For the 
configuration averaging (which we will denote by G) we use the simplest approximation 



G-G^G-G, (12.19) 

i.e. in all products the configurational-dependent quantities will be averaged separately. Taking 
into account Eqs. (|12.14p and (|12.19D . the averaged zeroth-order Green function G^. ^ is given 
by the well-known CPA solution for two-particle Green function in disordered metallic alloy^^ 

^^M = ^T, e ^^~* ) ^~* ) nfaM (12.20) 

where F 2 {k\,k 2 ) is given by 

F 2 (k 1 ,k 2 )^i(E(k 2 )-E(k 1 )) [ cic^M 1 )] 2 , (12.21) 

J OLO U — Zj{bj) — K{ki) 
for \E{k l )-E{k 2 )\^\nE(ki)% 

F 2 (4fe),^^5M) ( (12.22) 
for \E(h) - E(k 2 )\ y> WEih))]. 
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Here denotes the coherent potential and /(w) is the Fermi distribution function. The con- 

figurational averaged terms {a\a n ) are given by 



(ata n ) = e^-^Fxik), (12.23) 

k 

i r / i 



Fi(jfe) = ~J dujf(oj)lm 



oj - E(w) - 



After the configurational averaging equations (112. 16p and (I12.18P can be solved by Fourier trans- 
formation and we find 

G ijM (oj) = ^ E E e-^e^e-^e^Gih, k 2 ; k 3 ,k 4 ), (12.24) 

where 

F 2 (k u k 2 ) 



G(h, k 2 ; fe 3 , fc 4 ) = G(/ci, fe 2 ) = F 2 (fci, k 2 )5(k 4L , ki)5(k 3 ,k 2 

E 



X(q, k 2 )G(h,k 2 ) + F(£ h,k 2 )G(h - q, k 2 - q) 



and 



„ [FxCfci) - F x {k 2 -q)- u{q)F 2 {k x ,k 2 - qj^ih, k 2 - q))^ 

| X 1 (q, fci, k 2 )G{ki -q,k 2 -q)- Yi(q, ki)G(k 1 ,k 2 ) 
[FSl -0)- Fi(h) ~ Cj(q)F 2 (h - q, k 2 )](F 2 (h - q, &))-! 
-2 terms with u){q) -»■ -u(q), N(q) -»• (-1 - JV(g))Y (12.25) 

A (s> fc ) = ^ E e-^-^^Ci " i)- (12-26) 



Here the following notation were introduced 

X(q, k 2 ) = A(q, k 2 - q)A(-q, fc 2 )(^l(^2 - q) ~ 1 - (12.27) 

Y(<f, fc 2 ) = A(£ k 2 - q)A(-q, + N(q)), (12.28) 

*l(g, fci, fc 2 ) = A(q, h)A(-q, k 2 - q)(l + N(q) - F 1 (k 2 )), (12.29) 

Y 1 (q, h) = A(q, h)A(-q, h - - q) + N{q)). (12.30) 

Equation (|12.25f) is an integral equation for the Green function G(ki, k 2 ) to be determined. 
The structural averaged conductivity can be obtained, in principle, by using Eq. (|12.10p . where the 
Green function ((a\aj\ajai)) is to be replaced by Gijji(uj) as given by Eq. (ll2.24|) . It is, however, 
more convenient to start with the Kubo formula in the following formP^l 



ie 2 sr ^, 1 / dE(k) \ .. + . 

IT^oE^ {-dTP) <K«*+pl^»i., (12-31) 



where r/_ p = a^a^-p is the electron density operator. To find the Green function ((at.a,k+ p \'r]-p)) , 
the integral equation ()12.25j) has to be solved. In general, this can be done only numerically, but 
we can discuss here two limiting cases explicitly. At first we consider the weak-scattering limit 
being realized for a weak disorder in the alloy, and second, we investigate the temperature coeffi- 
cient of the conductivity for a strong potential scattering. 
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In the weak-scattering limit the CPA Green function is given by the expression 

F 2 (k u k 2 ) - i(E(k) - ECh))^ ■ (12.32) 
for - EftjK 

Corresponding to this limit the following solution ansatz for the Green function G(k,k + p) can 
be used 

e(*,* + p) = « !<W,- P » te * p)^j • p^pjg, , (12-33) 

where 7 describes the contribution of the electron-phonon scattering to the coherent potential. 
Taking into account that in the weak-scattering limit |S| <C ^(q), the terms F 2 (k,k — q) in the 
right hand side denominators of Eq. (|12.25|) can be replaced by the expression (|12.32p . and then 
the integral equation (|12.25p becomes for limp -4 

. df(E(k)) f dE{k) \ 1 . df(E{k)) f dEjk) \ 1 _ 

* dE(k) \ dk P J E(E(k)) + j(E(k)) ~ 1 dE(k) \ dk P ) E(£(fc)) 

1 1 V A(a k - a)A(-a k)\ Z^q) + Z 2 (k,q) 
X(E(k))N^ A[q ' k q)A[ q ' k)[ E {k)-E{k-q)-u J {q)+ie + 

Z 3 (k,q) - Z 4 (k,q) 



E(k)-E(k-q)+uj(q) + ie 

2 terms with uj(q) ^ -u(q), N{q) -4 (-1 - N(q))\. (12.34) 

Here the following notation were introduced 

KM = (f(.E(k -,))-!- W W)5 Xj . (12 . 35 ) 

(JEM , 

*<*, «) = (mkn + jvmjs^ • E(£(t -X.'W-»)) - (12 ' 36) 

= (i - W)) + ivta))^^ ' ^r^Fi ' (12 ' 37) 

*<*, ,) = (f(E(k - + ^ ■ (12 . 38 ) 

Approximating the self-energy terms 'E(E(k)) and j(E(k)) by Y>(Ef) = S and j(Ep) = 7, 
respectively, the terms proportional to E cancel and 7 can be calculated by 

df (dE(k) \ n 



7 



(1 - f(E(k - q)) + N(q)) dE{ f_ q) • (®|j£>p)]*(JE7(fc) - - g) - w( ff ))) 



^^(g,£-g>4(-g,£) 



q 



dE(k) 
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Using the approximations 



dE(k) 1 

-k, A(q, k - q)A(-q, k) ~ A 2 q, q 0, 



m 

where effective mass m* = m*(Ep), we find 

■> ~ ^ 2( 2 uX)^ I **^*> W <«> + "<*>) (12 ' 40) 

and 

7~<^ ' (12.41 

|T if T»fl fl . 

For a binary alloy A X B\_ X with concentrations of the constituents ca and cb and the corresponding 
atomic energies €a and e#, in the weak-scattering limit the coherent potential is given by 



Y> = c A c B {e A -e B fD(E F ). (12.42) 

Then the conductivity becomes 

" = 3w/*(^r)«' r - (12 - 43) 

where 

r _1 = E + 7, (12.44) 

in correspondence with the Matthiessen, Nordheim and Bloch-Gruneisen rulesP 
Now we estimate temperature coefficient of the conductivity for a strong potential scattering. 
For a strongly disordered alloy the electron-phonon interaction can be considered as a small 
perturbation and the Green functions G(k, k') on the right hand side of equation (|12.25p can 
be replaced by CPA Green functions F(k,k'). For simplicity, on the right hand side of equation 
(|12.25p we take into consideration only terms proportional to the Bose distribution function giving 
the main contribution to the temperature dependence of the conductivity. Then ((atafal^fc— fci)) 
becomes {k\ = k + p ~ k) 

{{a\a kl |r/ fe _ fcl » = F 2 {k, h) (l - ^yj^TT^ E A & 1 ~ ^ A (~Q, k)N(q) ■ 



{ 



F 2 (k, k - q)[F 2 {k -q,h-q)- F 2 (k, fci)](Fi(fc) - F x {k - q)) 



} ■ (12-45) 



[F x {k) - - g)] 2 - ^ 2 (g)F 2 (fc, fc - <?) 
F 2 (k- q,k)[F 2 (k - g , - g) - F 2 (fc, fci)](Fi(fc - g) - Fi( fc)). 
[Fx (A; - q) - - u; 2 (g)F 2 2 (A; - g, fe) 

Neglecting at low temperatures the terms w 2 (g)F|(/ci — g, &i) ~ g 4 as compared to [F\{k\) — 
Fx(k\ — q)\ 2 ~ g 2 and using equation ()12.2ip for uj(q) <C we find for small q and p — > 

((ata^l^)) ~ F 2 (k, k+p)[l+ (^|y) ' ^(A(g, k - q) - A(q, k))). (12.46) 
Here A(g, fc) is the temperature-dependent correction terms to the CPA Green function are given 

by 

A( S ,*)-2AVrto)(/ ^ H.^., w y)' (12.47) 
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For temperatures kpT <C Ep we can write 
and the conductivity becomes 

a = a CP A + Aa{T), (12.49) 

where 

is the standard CPA expression for the conductivity and 



Aa(r) ^ E( w E5J%) . 



n ^\ dk 

k q 

|Im U - m?) - m - J ' 4 ~ " m U - db - sw ) 14 )- (12 - 51) 

Introducing the effective mass of the electrons with E{k) ~ Ep, the temperature-dependent 
correction to the conductivity becomes 

Here the quantity Aa(T) is positive definite and increasing with increasing temperature. Hence, 
in strongly disordered alloys where the electron-phonon scattering is weak as compared with the 
disorder scattering the temperature coefficient of the resistivity is negative. It should be men- 
tioned, however, that the concrete temperature dependence of the correction term ()12.52p is a 
crude estimation only because in the derivation of (112. 52[) the influence of the disorder on the 
lattice vibrations has been neglected. 

One more remark is appropriate for the above consideration. For the calculation of transport 
coefficients in disordered 3d systems the classical approaches as the Boltzmann equation become 
useless if the random fluctuations of the potential are too large The strong potential fluc- 
tuations force the electrons into localized states. In order to investigate the resistivity of metallic 
alloys near the metal-insulator transitiorP^SK^H the corresponding formula for the resistivity can 
be deduced along the line described above. For a binary transition metal alloy the corresponding 
Hamiltonian is given by 

H = y~] tia\ai + ^ tija\aj (12.53) 

t ij 

(with €i = eA,ep depending on the occupation of the lattice site i). A corresponding integral 
equation for the Green function ((a_j|af ))u can be written down. Using a simple ensemble averaging 
procedure and approximating the averaged Green function by the expression 

((a>J)) w n±Y'exp[ik(R i -R j )]—^exp(-a(e k )\R i -R j \), (12.54) 

k u ~ e k 

the integral equation transforms into an equation for the parameter (a(efc)) -1 which is proportional 
to the averaged mean free path of the electrons. It can be shown then, by solving this equation for 
electrons at the Fermi surface Ep, that (a(Ep))~ 1 and the conductivity a drop in a discontinued 
way from (a)~^ n and a m in, respectively, to zero as the potential fluctuations exceed a critical 
value. Note that (a)~^ n is of the order 1/d, where d is the lattice parameter. 
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13 Discussion 



In the foregoing sections we have discussed some selected statistical mechanics approaches to the 
calculation of the electrical conductivity in metallic systems like transition metals and their dis- 
ordered alloys within a model approach. 

Electrons in metals are scattered by impurities and phonons. The theory of transport processes 
for ordinary metals was based on the consideration of various types of scattering mechanisms and, 
as a rule, has used the Boltzmann equation approach. The aim of the present review was to de- 
scribe an alternative approach to the calculation of electroconductivity, which can be suitable for 
transition metals and their disordered alloys. There is an important aspect of this consideration. 
The approximations used here are the tight-binding and modified tight-binding, which are admit- 
tedly not ideally precise but does give (at least as the first approximation) reasonable qualitative 
results for paramagnetic transition metals and their disordered alloys. We studied the electronic 
conduction in a model of transition metals and their disordered alloys utilizing the method of gen- 
eralized kinetic equations. The reasonable and workable expressions for the electrical conductivity 
were established and analyzed. We discussed briefly various approaches for computing electrical 
conductivity as well. 

We hope that these considerations have been done with sufficient details to bring out their scope 
and workability. In this paper, we have considered the idealized Hubbard model which is the 
simplest (in the sense of formulation, but not solution) and most popular model of correlated 
lattice fermionsP^^ We believe that this technique can be applied to other model systems (e.g. 
multi-band Hubbard model, periodic Anderson model, etc.). As it is seen, this treatment has 
some advantages in comparison with the standard methods of computing electrical conductivity 
within the Boltzmann equation approach, namely, the very compact form. The physical picture 
of electron-electron and electron-phonon scattering processes in the interacting many-particle sys- 
tems is clearly seen at every stage of calculations, which is not the case with the standard methods. 
This picture of interacting many-particle system on a lattice is far richer and gives more possi- 
bilities for the analysis of phenomena which can actually take place in real metallic systems. We 
believe that our approach offers a convenient way for approximate considerations of the resistivity 
of the correlated electron systems on a lattice. We believe that in view of the great difficulty 
of developing a first-principles microscopic theory of transport processes in solids, the present 
approach is a useful alternative for description the influence of electron-electron, electron-phonon 
and disorder scattering effects on the transport properties of transition metals and their disordered 
alloys. 

In the confines of a review of this nature many topics of great practical and theoretical interest 
have necessarily to be omitted (see e.g. RefsP^M^iJ). i n recent years the field of mesoscopic 
physics is developed rapidly G33l2Z2lEZS] j^. deals w jth systems under experimental conditions where 
several quantum length scales for electrons are comparable. The physics of transport processes 
in such systems is rich of quantum effects, which is typically characterized by interplay of quan- 
tum interference and many-body interactions. It would be of interest to generalize the present 
approach to quantum transport phenomena. 

In conclusion, the foregoing analysis suggests that the method of the generalized kinetic equa- 
tions is an efficient and useful formalism for the studying of some selected transport processes in 
metallic systems. 
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